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Abstract 


This  thesis  addresses  the  origin  and  effect  of  noise  and  fluctuations  in  quantum  devices 
that  find  applications  in  quantum  information  processing.  In  the  first  half  of  the  thesis, 
we  consider  a  phenomenological  model  of  the  noise  and  study  its  effect  on  the  fidelity  of 
quantum  measurements  and  operations  that  are  essential  for  quantum  computing.  We  fo¬ 
cus  primarily  on  quantum  gate  operations  in  phase  qubits,  detection  of  microwave  photons 
and  the  measurement  of  the  Berry  curvature.  Specifically,  we  examine  the  effect  of  the 
Ohmic  noise  on  optimally  controlled  flux-biased  phase  qubits  for  one-  and  two-quadrature 
microwave  pulses  and  demonstrate  that  two-quadrature  pulses  with  fixed  driving  frequency 
are  as  robust  as  variable  driving  frequency,  in  the  presence  of  environment.  Next,  we 
present  a  model  to  analyze  the  quantum  efficiency  of  a  microwave  photon  detector  based  on 
a  current-biased  Josephson  junction  and  study  the  effect  of  decoherence  on  the  detection  ef¬ 
ficiency  of  the  detector.  We  also  present  alternative  set-ups  for  microwave  photon  detection 
and  provide  a  systematic  method  to  compute  the  power  absorbed  by  the  detector.  We  then 
consider  the  effect  of  decoherence  on  the  Berry  curvature  measurement,  which  employs  a 
novel  non-adiabatic  protocol,  and  show  that  the  curvature  is  immune  to  decoherence. 

In  the  second  half  of  the  thesis,  we  go  beyond  the  phenomenological  models  of  the 
noises  and  perform  a  detailed  study  of  the  microscopic  description  of  the  Johnson  noise. 
Here  we  focus  primarily  on  quantum  dot  devices.  We  present  a  formalism  to  compute 
relaxation  rates  in  charge  and  spin  qubits  due  to  evanescent  wave  Johnson  noise  (EWJN) 
from  the  metallic  gates  that  are  in  proximity  to  quantum  dots.  The  EWJN  is  analyzed 
for  the  metallic  gates  that  are  characterized  by  both  local  and  nonlocal  dielectric  response 
functions.  Additionally,  we  extend  our  treatment  beyond  the  dipole  approximation  by 
taking  into  account  the  finite  size  of  the  quantum  dots.  We  also  derive  an  enhancement 
of  the  electric  EWJN  that  occurs  outside  a  thin  metallic  film,  relative  to  the  electric  field 
surrounding  a  conducting  half-space. 
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Chapter  1 

Introduction 


Physicists  have  long  envisioned  the  idea  of  harnessing  the  power  of  quantum  mechanics  to 
perform  computational  tasks  mm-  However,  the  actual  work  toward  practical  implemen¬ 
tation  of  such  a  ground-breaking  idea  has  only  begun  recently.  The  concept  of  quantum 
computation  has  gained  considerable  momentum  in  recent  years  due  to  a  discovery  that 
certain  computational  tasks,  such  as  number  factorization,  period  finding  [3],  and  database 
search  [3],  could  be  speeded  up  if  these  problems  were  encoded  in  the  states  of  quantum 
systems  called  quantum  bits  (or  qubits)  and  computed  employing  operations  described  by 
quantum  mechanics  [5].  Due  to  its  quantum  nature,  a  quantum  bit  can  be  in  a  superpo¬ 
sition  of  both  its  logical  states  “0”  and  “1”.  Hence,  while  a  classical  memory  consisting 
of  n  bits  allows  to  store  one  out  of  2n  numbers,  n  quantum  bits  can  store  all  2n  numbers 
simultaneously.  Logical  operations  on  a  quantum  memory  work  in  parallel  to  all  these 
numbers,  whereas  a  classical  computer  would  need  to  repeat  the  calculation  2n  times.  This 
“built-in”  parallelism  facilitates  a  quantum  computer  to  work  exponentially  faster  allowing 
to  solve  computationally  intensive  problems  in  seconds  which  would  otherwise  take  years  on 
today’s  fastest  supercomputers.  Furthermore,  in  a  quantum  processor,  qubit  states  can  be 
superposed  and  entangled  by  sequences  of  externally  controlled  manipulations  that  result 
in  unitary  transformations,  also  known  as  quantum  gates.  Such  quantum  operations  allow 
quantum  computers  to  push  boundaries  of  parallel  computations,  unconditionally  secure 


2 


communications,  etc. 

Numerous  systems  have  now  been  proposed  as  promising  candidates  for  quantum  infor¬ 
mation  processing.  Among  these  are  trapped  ions  EE,  optical  lattices  EE,  photons  [TO], 
nuclear  spins  of  molecules  (nuclear  magnetic  resonance,  Nitrogen-vacancy  centers  in  dia¬ 
mond,  etc.  )  mm ,  quantum  dots  m  and  superconducting  circuits  [14] .  While  there  exists 
different  platforms  on  which  to  build  the  hardware  of  a  quantum  computer,  semiconduct¬ 
ing  and  superconducting  qubits  based  on  solid-state  technologies  are  the  most  promising 
candidates  that  hold  unique  prospects  of  scalability  and  manageability  owing  to  their  ease 
of  integration  with  highly  developed  semiconductor  technology  and  classical  electronics. 

For  practical  realization  of  a  quantum  computer,  advances  must  be  made  to  preserve 
the  quantum  coherence  much  longer  than  gate  operation  times.  The  loss  of  the  quan¬ 
tum  coherence  of  a  qubit  stems  from  the  interaction  of  the  qubit  with  its  local  environ¬ 
ment.  This  phenomenon  is  often  known  as  decoherence  in  the  literatures.  In  the  past 
decade,  tremendous  progress  has  been  made  in  identifying  various  sources  of  decoherence 
that  plague  practical  implementation  of  solid-state  qubits  m •  In  a  semiconducting  qubit, 
some  sources  of  decoherence  are  intrinsic  in  nature,  such  as  charge  coupling  to  phonons, 
while  others  are  extrinsic,  such  as  Johnson  noise  due  to  fluctuating  electromagnetic  fields 
arising  from  the  metallic  gates  or  transmission  lines  used  to  control  the  qubit.  Similarly,  in 
superconducting  qubits,  the  most  notorious  of  all  sources  of  decoherence  is  the  presence  of 
defects  in  dielectric  materials  US]  used  to  construct  Josephson  junctions,  which  are  building 
blocks  of  superconducting  qubits.  Indeed,  decoherence  has  been  a  constant  worry  since  the 
beginning  of  the  field  of  quantum  computing,  which  has  led  to  emergence  of  new  ideas  like 
decoherence  free  subspace  HZ!,  topological  quantum  computing  [ISi  19],  error-correction 
codes  ED],  dynamical  decoupling  m,  which  have  all  galvanized  the  quest  for  harnessing 
the  power  of  quantum  mechanics  to  perform  complex  computational  tasks. 

Since  it  is  not  practically  possible  to  completely  isolate  a  qubit  from  its  environment, 
controlling  the  qubit  in  the  presence  of  environment  is  another  equally  important  aspect 
of  quantum  computation.  Furthermore,  qubits  are  often  approximated  by  a  two-level  sys- 
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tem,  which  is  not  necessarily  accurate  particularly  in  case  of  solid-state  qubits  based  on 
macroscopic  degrees  of  freedom  like  charge  or  current  as  in  charge  or  phase  qubit  in  super¬ 
conducting  qubits.  This  poses  another  hindrance  to  accurate  control  of  qubits  during  gate 
operations.  Inaccurate  control  often  leads  to  additional  error,  which  can  surpass  the  thresh¬ 
old  of  a  fault-tolerant  quantum  computation  [22|  [23] ,  rendering  the  output  from  quantum 
operations  unreliable.  One  way  to  overcome  this  problem  is  by  engineering  shaped  pulses 
that  are  optimized  for  particular  gate  operations  and  are  easy  to  implement  in  the  labora¬ 
tory.  Therefore,  it  is  vitally  important  to  perform  a  detailed  study  of  the  performance  of 
the  optimally  shaped  pulses  in  the  presence  of  environment  that  qubits  are  most  likely  to 
encounter. 

While  decoherence  is  the  most  pressing  issue  and  deserves  undivided  attention  of  the 
scientific  community,  researchers  have  also  managed  to  explore  several  other  avenues  within 
mesoscopic  physics  which  will  ultimately  benefit  the  field  of  quantum  computing.  One  of 
such  avenues  is  circuit  quantum  electrodynamics  (cQED),  which  also  finds  applications  in 
quantum  information  processing.  Indeed,  with  the  generation  of  non-classical  propagating 
fields  p31  [25] .  and  the  demonstrated  potential  for  microwave  radiation  to  mediate  coupling 
between  solid-state  qubits  [26] ,  the  field  of  cQED  has  become  a  rapidly  evolving  test-bed 
for  mesoscopic  physics  and  a  promising  paradigm  for  quantum  computation.  However, 
the  lack  of  a  microwave  photon  detector  is  one  major  deficiency  within  this  held.  High 
efficiency  microwave  photon  detection  is  an  interesting  problem  for  many  areas  of  physics, 
including  low  temperature  measurement,  quantum  information  science,  particle  physics 
and  astrophysics.  For  optical  photons,  there  are  many  examples  of  devices  capable  of 
detecting  single  photons  with  high  efficiency  m-  However,  reliable  single  photon  detection 
of  microwaves  is  very  difficult,  principally  due  to  their  low  energy.  It  is  therefore  imperative 
to  explore  different  models  of  on-chip  microwave  photon  detector  and  analyze  their  efficiency 
in  the  presence  of  coupling  to  the  environment. 

Aware  of  the  difficulty  in  overcoming  decoherence  in  a  practical  quantum  computer 
architecture,  researchers  have  also  started  exploring  alternative  proposals  for  qubit  imple- 
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mentation  that  are  fundamentally  different  from  the  previous  proposals.  In  these  proposals, 
information  is  encoded  in  those  degrees  of  freedom  that  are  insensitive  to  local  perturbations 
and  hence  are  protected  against  errors  caused  by  interaction  with  the  local  environment. 
This  novel  computing  paradigm,  known  as  topological  quantum  computing,  has  gained 
considerable  attention  in  recent  years  due  to  its  intrinsic  fault-tolerant  nature.  It  has  been 
shown  that  universal  quantum  computation  can  be  performed  using  only  geometric  effects 
appearing  in  quantum  phases  [28] .  The  Berry  phase  is  one  of  such  geometric  phases  that 
can  be  used  to  implement  elementary  quantum  gates.  Although  geometric  evolutions  are 
easy  to  control  and  resistant  to  error,  it  is  not  quite  straightforward  to  measure  the  Berry 
phase  in  practice.  Traditional  experiments  for  the  Berry  phase  measurement  require  strict 
adiabatic  conditions  and  rely  on  interference  effects  [29].  The  measurement  process  would 
be  greatly  simplified  if  such  stringent  conditions  could  be  avoided  altogether.  A  more  recent 
proposal  [30]  that  does  not  require  strict  adiabatic  condition  is  a  promising  step  in  that 
direction.  Although  non-adiabatic  measurement  protocols  are  very  attractive,  they  may 
not  be  immune  to  noises  from  the  environment.  It  is,  therefore,  important  to  understand 
how  the  coupling  of  a  qubit  to  its  environment  affects  the  Berry  curvature,  and  in  turn  the 
Berry  phase,  in  non-adiabatic  measurement  protocols. 

To  that  end,  this  thesis  consists  of  two  major  parts.  The  first  part,  comprising  of 
the  first  four  chapters,  discusses  issues  related  to  superconducting  qubits.  In  particular,  we 
examine  the  optimally  controlled  phase  qubit  in  the  presence  of  the  Ohmic  environment  and 
propose  alternative  pulse  sequences  that  are  robust  against  the  presence  of  the  third  level 
in  the  qubit  and  easier  to  implement  in  the  laboratory.  We  also  analyze  the  efficiency  of  a 
microwave  photon  detector  based  on  a  current-biased  Josephson  junction.  Next,  we  consider 
the  effect  of  environment  on  the  Berry  curvature  and  discuss  the  prospects  of  measurement 
of  the  Berry  curvature  in  superconducting  qubits.  The  second  part,  consisting  of  the  next 
two  chapters,  discusses  issues  related  to  semiconducting  qubits.  Here  we  perform  a  detailed 
study  of  the  physics  of  Johnson  noise  and  provide  formula  for  relaxation  rates  in  charge  and 
spin  qubits  for  various  cases.  The  majority  of  the  work  done  in  this  part  has  an  emphasis  on 
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quantum  dot  devices,  but  it  applies  to  other  qubit  architectures  that  incorporate  metallic 
gates,  traps  or  electrodes  in  the  vicinity  of  the  qubit. 

1.1  Thesis  outline 

This  thesis  describes  several  aspects  of  the  theoretical  issues  related  to  semiconducting  and 
superconducting  qubits.  A  common  thread  that  ties  together  different  chapters  in  this 
thesis  is  the  study  of  the  effect  of  decoherence  on  measurements  and  operations  in  quantum 
devices  that  find  applications  in  quantum  information  processing. 

In  chapter  [2j  we  provide  a  brief  overview  of  different  noise  sources  that  are  responsible 
for  dissipation  and  dephasing  in  superconducting  and  semiconducting  qubits.  Some  of  these 
noise  sources  are  due  to  external  circuitry  that  is  used  to  control  and  manipulate  qubits, 
while  the  other  sources  are  intrinsic  to  qubit  devices.  We  then  proceed  to  discuss  different 
methods  that  can  be  used  to  study  the  coupling  of  a  qubit  to  its  environment.  Here  we 
derive  the  master  equation  of  an  open  quantum  system  using  a  phenomenological  model  of 
the  environment. 

In  chapter  [3j  we  study  the  effect  of  environment  on  the  gate  operation  of  flux-biased 
phase  qubits.  We  employ  the  master  equation  for  a  reduced  density  matrix  of  the  qubit 
system  coupled  to  the  Ohmic  environment,  described  by  the  Caldeira-Leggett  model.  Nu¬ 
merically  solving  this  equation,  we  evaluate  the  gate  error  as  a  function  of  gate  time, 
temperature,  and  environmental  coupling  strength  for  experimentally  determined  qubit  pa¬ 
rameters.  The  analysis  is  presented  for  single-quadrature  microwave  (control)  pulses  as  well 
as  for  two-quadrature  pulses,  which  lower  the  gate  error  significantly  for  idealized  systems 
in  the  absence  of  environment.  Our  results  indicate  that  two-quadrature  pulses  with  fixed 
and  variable  driving  frequency  have  similar  performance,  which  outweighs  the  performance 
of  single-quadrature  pulses,  in  the  presence  of  environment. 

In  chapter  [4j  we  analyze  the  quantum  efficiency  of  a  microwave  photon  detector  based 
on  a  current-biased  Josephson  junction.  We  consider  the  Jaynes-Cummings  Hamiltonian 
to  describe  coupling  between  the  photon  field  and  the  junction.  We  then  take  into  ac- 
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count  the  coupling  of  the  junction  and  the  resonator  to  the  environment.  We  solve  the 
equation  of  motion  of  the  density  matrix  of  the  resonator-j  unction  system  to  compute  the 
quantum  efficiency  of  the  detector  as  a  function  of  detection  time,  bias  current,  and  energy 
relaxation  time.  Our  results  indicate  that  junctions  with  modest  coherence  properties  can 
provide  efficient  detection  of  single  microwave  photons,  with  quantum  efficiency  in  excess 
of  80%.  Furthermore,  we  discuss  other  alternative  set-ups  for  microwave  photon  detection 
and  provide  a  systematic  way  of  calculating  the  power  absorbed  by  the  detector  in  those 
set-ups. 

Next,  we  discuss  the  effect  of  decoherence  on  geometric  phases  in  quantum  mechanics. 
In  chapter  [5j  we  perform  a  detailed  analysis  of  the  effect  of  relaxation  and  pure  dephasing 
on  the  Berry  curvature,  which  is  computed  by  employing  a  novel  non-adiabatic  protocol. 
We  then  provide  an  experimentally  convenient  method  to  measure  the  Berry  curvature. 
Our  results  indicate  that  the  Berry  curvature  is  immune  to  decoherence  in  non-adiabatic 
measurement  protocols.  In  addition,  we  also  consider  the  effect  of  the  third  level  on  the 
measurement  of  the  Berry  curvature  in  superconducting  qubits. 

In  the  remaining  chapters,  we  turn  our  attention  to  microscopic  models  of  decoherence 
in  qubits,  focusing  primarily  on  semiconducting  qubit  architectures.  In  chapter  [6j  we  focus 
on  the  origin  and  effect  of  evanescent-wave  Johnson  noise  in  semiconducting  qubits.  In 
many  quantum  computer  architectures,  the  qubits  are  in  close  proximity  to  metallic  device 
elements.  Metals  have  a  high  density  of  photon  modes,  and  the  fields  spill  out  of  the  bulk 
metal  because  of  the  evanescent-wave  component.  Thus  thermal  and  quantum  electromag¬ 
netic  Johnson-type  noise  from  metallic  device  elements  can  decohere  nearby  qubits.  Here 
we  use  quantum  electrodynamics  to  compute  the  strength  of  this  evanescent- wave  John¬ 
son  noise  as  a  function  of  distance  from  a  metallic  half  space.  Previous  treatments  have 
shown  unphysical  divergences  on  the  surface  of  metallic  gate.  We  remedy  this  by  using  a 
proper  nonlocal  dielectric  function.  Decoherence  rates  of  local  qubits  are  proportional  to 
the  magnitude  of  electric  or  magnetic  correlation  functions  evaluated  at  the  qubit  position. 
These  rates  serve  as  an  important  constraint  on  future  device  architectures.  Compari- 
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son  with  single-electron  spin  relaxation  measurements  shows  that  evanescent-wave  Johnson 
noise  may  constitute  the  dominant  relaxation  mechanism  in  experiments  performed  at  low 
magnetic  field. 

In  chapter  [7J  we  continue  our  study  of  decoherence  in  charge  (spin)  qubits  due  to 
evanescent-wave  Johnson  noise  (EWJN)  in  a  laterally  coupled  double  quantum  dot  (sin¬ 
gle  quantum  dot).  We  derive  expressions  for  the  energy  relaxation  rates  of  charge  and 
spin  qubits  in  a  variety  of  dot  geometries,  and  EWJN  is  shown  to  be  a  dominant  source 
of  decoherence  for  spin  qubits  held  at  low  magnetic  fields.  Previous  studies  in  this  field 
approximated  the  charge  or  spin  qubit  as  a  point  dipole.  Ignoring  the  finite  size  of  the 
quantum  dot  in  this  way  leads  to  a  spurious  divergence  in  the  relaxation  rate  as  the  qubit 
approaches  the  metal.  Our  approach  goes  beyond  the  dipole  approximation  and  remedies 
this  unphysical  divergence  by  taking  into  account  the  finite  size  of  the  quantum  dot.  Ad¬ 
ditionally,  we  derive  an  enhancement  of  EWJN  that  occurs  outside  a  thin  metallic  him, 
relative  to  the  held  surrounding  a  conducting  half-space. 

Finally,  in  chapter  [8]  we  offer  concluding  remarks. 

1.2  Publications 

Some  of  the  main  chapters  in  this  thesis  are  based  on  published  results,  each  of  which 
represents  the  work  of  many  individuals.  Here,  I  document  my  contribution  to  each  work, 
as  well  as  my  collaborators. 

Chapter  [3]  is  based  on  Ref.  [31],  titled  Effect  of  the  Ohmic  environment  on  an  optimally 
controlled  flux-biased  phase  qubit ,  which  I  completed  with  Prof.  Maxim  G.  Vavilov.  I  carried 
out  the  analytical  and  numerical  calculations,  under  the  supervision  of  Prof.  Vavilov. 

Chapter  [4]  is  based  on  Ref.  [32j,  titled  Quantum  efficiency  of  a  microwave  photon  de¬ 
tector  based  on  a  current-biased  Josephson  junction ,  which  I  completed  with  Prof.  Robert 
McDermott,  and  Prof.  Maxim  G.  Vavilov.  I  carried  out  the  main  numerical  calculation 
based  on  a  theoretical  model  suggested  by  Prof.  Vavilov,  which  was  based  on  an  experi¬ 
mental  device  first  proposed  and  developed  by  Prof.  McDermott. 
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Chapter  [6]  is  based  on  Ref.  [33],  titled  Qubit  relaxation  from  evanescent-wave  John¬ 
son  noise,  which  I  completed  with  Dr.  Luke  Langsjoen,  Prof.  Maxim  G.  Vavilov,  and 
Prof.  Robert  Joynt.  I  principally  carried  out  the  numerical  work  with  assistance  from  Dr. 
Langsjoen,  based  on  theoretical  work  laid  out  by  Prof.  Vavilov,  and  Prof.  Joynt. 

Chapter  [7]  is  based  on  Ref.  pH] ,  titled  Relaxation  in  quantum  dots  due  to  evanescent- 
wave  Johnson  noise,  which  I  completed  with  Dr.  Luke  Langsjoen,  Prof.  Maxim  G.  Vavilov, 
and  Prof.  Robert  Joynt.  I  carried  out  the  analytical  and  numerical  work  with  assistance 
from  Dr.  Langsjoen,  Prof.  Vavilov,  and  Prof.  Joynt. 
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Chapter  2 

Qubit  devices  and  decoherence 

2.1  Introduction 

Most  proposals  for  implementing  a  quantum  computer  are  based  on  qubits  constructed  from 
microscopic  degrees  of  freedom.  They  consider  spin  of  either  electrons  or  nuclei,  or  transition 
dipoles  of  either  atoms  or  ions  in  vacuum.  These  degrees  of  freedom  are  naturally  very  well 
isolated  from  their  environment,  and  hence  decohere  slowly.  However,  in  other  solid-state 
qubit  architectures  like  integrated  superconducting  circuits  or  semiconducting  heterostruc¬ 
tures,  quantum  coherence  inevitably  suffers  from  fluctuations  due  to  their  environment.  In 
this  chapter,  we  first  briefly  discuss  different  sources  of  noises  that  plague  semiconducting 
and  superconducting  qubit  architectures.  We  then  discuss  theoretical  methods  or  models 
that  are  used  in  this  thesis  to  study  coupling  of  the  qubit  devices  to  environment,  within 
the  framework  of  open  quantum  systems. 

2.2  Decoherence  in  superconducting  qubits 

Superconducting  qubits  are  solid-state  electrical  circuits  based  on  the  Josephson  junction 
(JJ),  which  is  non-dissipative  and  strongly  non-linear  circuit  element  [35 j .  In  contrast  to 
microscopic  entities  such  as  spins  or  atoms,  they  tend  to  be  well  coupled  to  other  circuits, 
which  make  them  very  attractive  candidates  for  readouts  and  gate  implementations.  How- 
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ever,  the  same  biasing  wires  that  make  superconducting  circuits  easy  to  manipulate  and 
measure  also  couple  the  qubit  to  the  electromagnetic  environment  which  is  responsible  for 
dissipation  and  dephasing  of  the  qubit.  Any  external  circuit  that  is  connected  to  the  JJ 
provides  complex  impedance  Z{u)  to  the  Josephson  junction.  This  external  impedance  is 
defined  as  the  ratio  of  the  complex  amplitudes  of  the  voltage  response  of  the  environment 
to  the  the  oscillating  current  at  frequency  oj.  The  real  part  of  this  complex  impedance 
characterizes  the  dissipative  contribution  of  the  external  electromagnetic  environment.  In 
this  and  subsequent  chapters,  we  primarily  focus  on  the  noise  due  to  external  circuits. 

In  addition  to  the  noise  due  to  external  circuits,  various  other  mechanisms  are  responsi¬ 
ble  for  noises  in  superconducting  qubits.  Most  of  these  mechanisms  are  intrinsic  in  nature. 
For  instance,  the  motion  of  charges  in  the  defects  in  the  oxide  tunnel  barrier  and  amorphous 
dielectric  of  the  circuit  have  been  proposed  to  be  the  source  of  critical  current  noise  |36l  137] 
and  charge  noise  [38,  [39] ,  respectively.  The  1//  spectrum  (where  /  is  the  frequency)  of 
the  materials  noises  suggests  that  they  all  originate  from  two-level  fluctuators  in  the  amor¬ 
phous  tunnel  layer  of  the  junction  m-  Similarly,  the  relaxation  of  paramagnetic  spins 
located  at  the  superconductor-insulator  interfaces  can  be  the  source  of  low  frequency  1  // 
flux  noise  mm,  although  the  exact  mechanism  responsible  for  the  origin  of  such  spins  is 
not  quite  clear  and  is  an  area  of  active  research. 

2.3  Decoherence  in  semiconducting  qubits 

The  building  block  of  a  semiconducting  qubit  is  a  quantum  dot,  which  can  be  thought  of 
as  a  quantum  box  that  can  be  filled  with  electrons  (or  holes)  which  occupy  the  available 
discretized  states  of  the  system.  The  electrons  can  tunnel  on  and  off  the  dot,  which  is 
coupled  to  a  large  reservoir  via  tunnel  barriers.  The  gate  electrodes  are  used  to  control  the 
height  of  the  barriers  and  consequently  the  rates  for  tunneling  through  the  barriers  on  and 
off  the  dots  [13].  Although  many  kinds  of  quantum  dots  have  been  realized  so  far,  in  this 
thesis  we  focus  exclusively  on  GaAs  or  Si  based  quantum  dots. 

An  important  class  of  proposals  consists  of  using  the  charge  degree  of  freedom  in  semi- 
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conducting  double  quantum  dots  to  realize  a  charge  qubit  Hang.  A  significant  advantage 
of  the  charge  qubit  is  that  it  can  be  controlled  directly  via  external  voltage  sources.  While 
charge  qubits  are  easy  to  control  and  manipulate,  they  also  couple  easily  to  the  outside  en¬ 
vironment.  It  is  widely  believed  that  the  main  contribution  to  decoherence  in  a  charge  qubit 
comes  from  its  coupling  to  fluctuating  background  charges  (FBC)  that  create  dynamical 
electric  field,  which  affects  the  charge  qubit  state  [J6,  07J  08] .  The  sources  of  FBCs  are  elec¬ 
trons  trapped  in  the  insulating  layer  close  to  the  surface  of  metallic  gates/electrodes.  These 
electrons  can  also  interact  with  the  conduction  electrons  in  the  gates  leading  to  significant 
decoherence  of  charge  qubits  |  i!)  .  Other  probable  alternative  microscopic  mechanisms  of 
decoherence  in  charge  qubits  are  via  electron  coupling  to  the  phonon  bath  [50]  and  due  to 
gate  voltage  fluctuations  EH- 

Other  proposals  utilize  the  electron  spin  in  semiconductor  quantum  dots,  which  can 
be  isolated  and  controlled  with  a  high  accuracy,  to  realize  a  spin  qubit  [52l  153].  However, 
like  a  charge  qubit,  spin  qubit  also  suffers  from  decoherence  due  to  its  coupling  with  the 
surrounding  environment.  Two  kinds  of  environment  have  been  identified  as  the  main 
source  of  decoherence  of  an  electron  spin  in  a  quantum  dot.  They  are  the  phonons  in 
the  lattice  [53,  (55j  and  the  spins  of  atomic  nuclei  in  the  quantum  dot  [56].  In  GaAs 
based  quantum  dots,  strong  spin-orbit  interaction  is  present  even  in  the  absence  of  external 
electric  fields.  This  interaction  couples  spin  with  the  orbital  degrees  of  freedom,  which  is 
then  coupled  to  phonons.  Consequently,  it  leads  to  indirect  coupling  between  the  electron 
spin  and  phonons  m-  The  phonons  constitute  a  large  dissipative  bosonic  reservoir  and 
provide  a  source  of  dephasing  and  relaxation.  Short  time  correlations  in  the  phonon  bath 
induce  a  Markovian  dynamics  of  the  electron  spin. 

The  electron  spin  and  nuclear  spins  in  a  quantum  dot  couple  via  the  hyperfine  inter¬ 
action,  which  creates  entanglement  between  them  and  strongly  affects  the  electron  spin 
dynamics  [58j  [59] .  It  turns  out  that  long  time  correlations  in  the  nuclear  spin  system  in¬ 
duce  a  non-Markovian  dynamics  of  the  electron  spin,  with  non-exponential  decay  in  time 
of  the  expectation  values  of  the  electron  spin  components. 
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In  addition  to  the  sources  mentioned  above,  semiconducting  qubits  can  also  decohere 
due  to  Johnson  noise  from  the  metallic  gates  that  are  in  the  vicinity  of  quantum  dots.  The 
high  density  of  evanescent  modes  (both  electric  and  magnetic)  nearby  the  metallic  gates 
easily  couple  to  quantum  dots  causing  the  loss  of  their  phase  coherence  and  consequently 
affecting  both  charge  and  spin  qubits. 

2.4  Modeling  the  dynamics  of  open  quantum  systems 

Historically,  description  of  open  system  dynamics  has  largely  been  based  on  the  reduced  den¬ 
sity  matrix  formalism  within  which  both  intrinsic  quantum  mechanical  fluctuations  of  the 
system  and  external  noise  of  its  environment  can  be  conveniently  incorporated  in  a  unified 
manner.  During  the  past  few  decades,  efforts  have  been  made  to  devise  various  phenomeno¬ 
logical  models  and  more  recently,  microscopic  Hamiltonians  have  been  used  to  describe  an 
open  quantum  system.  In  this  section,  we  outline  two  theoretical  models/methods  that  are 
used  to  study  the  dynamics  of  a  qubit  device  coupled  to  its  environment  in  this  thesis. 

Lindblad  master  equation 

In  may  cases,  it  is  useful  to  model  the  dynamics  of  an  open  quantum  system  by  means  of 
an  appropriate  equation  of  motion  for  its  density  matrix,  which  is  known  as  the  quantum 
master  equation.  Here  we  consider  the  evolution  of  the  density  matrix  of  an  open  quantum 
system,  without  providing  any  connections  to  underlying  microscopic  physical  model  of  the 
environment  with  which  the  system  is  coupled.  This  leads  to  a  phenomenological  descrip¬ 
tion  of  dephasing  and  relaxation  in  the  open  quantum  system,  often  known  as  Lindblad 
dissipator. 

Consider  a  system-reservoir  coupling  which  is  uncorrelated  at  t  =  0,  namely: 

W(0)  =  p{0)®WR(0)t  (2.1) 

where  W,p ,  and  Wr  are  density  matrices  of  the  entire  system  (qubit  and  its  environment), 
a  qubit  device  and  a  reservoir,  respectively.  Next,  define  a  dynamical  map  M(t)  describing 
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a  transformation  of  the  reduced  system  at  t  =  0  to  some  time  t  >  0  as: 


p(t)  =  M(t)p(  0) 


Trfl 


U(t,  0)W(0)U*(t,  0) 


(2.2) 


where  U(t,0)  is  a  unitary  operator  that  defines  the  evolution  of  the  qubit  and  its  envi¬ 
ronment  combined.  The  map  A4  represents  a  completely  positive  and  trace  preserving 
quantum  operation.  If  t  is  allowed  to  vary  in  Ai(t),  then  it  leads  to  a  family  of  one  param¬ 
eter  dynamical  maps  with  Ad(0)  being  the  identity  map.  If  the  characteristic  time  scales 
over  which  the  reservoir  correlation  function  decays  are  much  smaller  than  the  characteris¬ 
tic  time  scales  of  the  qubit  device,  memory  effects  in  the  dynamics  of  the  reduced  density 
matrix  equation  can  be  neglected.  This  leads  to  a  Markovian  type  behavior  and  may  be 
formalized  with  the  help  of  a  semigroup  property: 


M(t\)  M(t2)  =  M(t\  +  f2)  ,h  ,t2  >  0  . 


(2.3) 


Thus  a  quantum  dynamical  semigroup  is  a  continuous,  one-parameter  family  of  dynamical 
maps  satisfying  the  semigroup  property  (Markovian).  If  the  quantum  dynamical  semigroup 
Ai  (f)  is  contracting,  then  there  exists  a  linear  map  L  called  a  generator  of  the  semigroup  [60| : 


Ai(t)  =  exp(Ct) ,  which  leads  to 

di  =  cm-  (2.4) 

The  construction  of  the  most  general  form  of  the  generator  L  leads  to  the  Lindblad  master 
equation.  In  the  following,  we  present  a  derivation  of  the  Lindblad  form  of  equation  of 
motion  for  the  reduced  density  matrix  of  the  qubit  system,  following  Ref.  [60j. 

The  Hamiltonian  of  a  global  system  is  given  by: 


H  =  Hq  +  HR  +  V, 


(2.5) 


where  Hq  and  Hr  are  the  qubit  and  reservoir  Hamiltonians,  respectively  and  V  is  the 
interaction  Hamiltonian  between  the  qubit  and  the  reservoir.  In  the  interaction  picture, 
equation  of  motion  of  a  full  density  matrix  of  the  entire  system  W  is  given  by  the  Von- 
Neumann  equation: 


dWr(t) 

dt 


Vl(t),W!(t) 


(2.6) 
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where  the  subscript  /  denotes  the  interaction  picture.  The  formal  solution  of  the  above 
equation  is  given  by: 

l  rl 

WI(t)  =  WI(0)  +  -J  ds[yKs),W7(*)].  (2.7) 

Inserting  the  integral  equation  into  the  Von-Neumann  equation  and  tracing  the  bath  vari¬ 
ables  out  we  obtain: 

^  Jo  ^  [^(SW(S)]]  ds’  (2-8) 

where  we  assumed  Tr#  Vj(t),Wj(0)  =  0.  The  right  hand  side  of  the  above  equation  still 
depends  on  the  full  density  matrix,  and  this  is  where  we  make  the  first  approximation. 
Within  the  Born  approximation,  it  is  assumed  that  the  coupling  between  the  system  and 
the  reservoir  is  negligibly  weak  so  that  the  reservoir  density  matrix  is  negligibly  affected  by 
the  interaction,  and  the  total  system  after  some  time  t  is  Wj{t )  =  pi(t)  <8>  WR,  where  WR 
is  the  density  matrix  of  the  reservoir.  This  gives: 

d^  =  ~p{ri 'rR[Vi(t),[Vi(s),pI(s)®Wlt]\  ds.  (2.9) 

A  further  simplification,  pi(s)  — >•  pi(t),  leads  to  the  evolution  of  the  qubit  density 
matrix  at  time  t  dependent  only  on  the  present  state.  This  assumption  leads  to  the  Redfield 
equation,  which  is  local  in  time,  but  depends  on  the  choice  of  the  initial  preparation  at  t  =  0. 
Next,  time-scale  is  coarse  grained  under  an  assumption  that  environmental  excitations  decay 
over  times  which  are  not  resolved.  In  other  words,  the  reservoir  correlation  is  much  smaller 
than  the  time  scale  over  which  the  qubit  state  varies  appreciably.  This  allows  us  to  replace 
s  by  t  —  s  and  send  the  upper  limit  of  integration  to  oo  in  the  above  equation.  Such  an 
approximation  is  often  known  as  the  Markovian  approximation  and  it  leads  to  the  following 
master  equation: 

do  1  /*°°  r  r  - 

=  -jp]  Trfi  [t P}(t),  [Vi(t  -  s),  Pi(t)  <8>  WjR\  J  ds  .  (2.10) 

So  far  we  have  performed  the  Born-Markov  approximation.  However,  these  approxima¬ 
tions  do  not  guarantee  a  quantum  dynamical  semigroup  evolution.  A  further  approximation 
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involving  the  averaging  over  the  rapidly  oscillating  terms  in  the  master  equation  must  be 
performed.  Such  an  approximation  is  often  known  as  the  rotating  wave  or  secular  approximation. 
In  order  to  demonstrate  this,  we  first  decompose  the  interaction  Hamiltonian  into  eigenop- 
erators  of  the  system  Hamiltonian: 

Vj(t)  =  ^2  e~ZU}tAa{w)  <8>  Ba(t) ,  where 

a, CU 

Aa(uj)  =  ^2  n(e)HQfl(e/)  and 

e'—e=c j 

Ba(t)  =  eiflRt/n  Bae-^Rt/n  .  (2.11) 


Here  Aa,Ba  denote  operators  belonging  to  the  system  and  reservoir  in  the  Schrodinger 
picture.  n(e)  is  the  projection  operator  to  the  energy  eigenstate  of  the  system  with  energy 
e.  Here  we  assumed  discrete  energy  spectrum  of  the  system  Hamiltonian.  An  immediate 
consequence  of  the  above  definition  is  that  the  following  relation  is  satisfied:  [Hq,  Aa(u>)]  = 
—ujAa{pj),  that  is,  Aa(uj)  lowers  the  qubit  energy  by  u  while  A'( u)  raises  it  by  u.  Inserting 
this  form  of  interaction  Hamiltonian  into  the  Born-Markov  master  equation,  Eq.  [2.10],  we 
obtain: 


dpi 

dt 


K 2 


EE 

lj,cj'  a, (3 


Xa,p( u') 


Ap(uj)pi(t)  Al(u')  -  ii(w')  Ap(u>)  Pi(t ) 


+  h.  c , 

(2.12) 


where  Xa,p{A)  is  the  one-sided  Fourier  transform  of  the  reservoir  correlation  functions  given 
by: 


Xa,p(A)=  /  dselU}' s  {Bl(t)  Bp(t  -  s)) ,  where 


{Bl(t)  Bp(t  -  s))  =  TrR[Bl(t)  Bp{t  -  s)  Wj 


f] 


(2.13) 


If  W P  is  a  stationary  state  of  the  reservoir,  [Vi,Wp]  =  0,  then  the  reservoir  correlation 
functions  are  homogeneous  in  time 


(BUt)Bp(t-s))  =  (Bl(s)Bp(0)). 


(2.14) 
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In  the  above  equation,  \oj— u/|_1,  where  oj  /  a/,  defines  the  typical  time  scale  associated  with 
the  intrinsic  evolution  of  the  qubit  device.  If  the  characteristic  time  scale  is  short  compared 
to  the  relaxation  time  of  the  qubit  device,  the  non-secular  terms,  that  is,  those  for  which 
oj  ^  oj'  may  be  neglected  since  they  oscillate  very  rapidly  during  the  relaxation  time  when 
the  qubit  state  varies  appreciably.  This  is  called  the  rotating  wave  approximation.  With 
this,  the  evolution  equation  becomes: 

,/3  M  As  MA  (f )  AL  M  -  Ai  (u)  Ap  (u 

In  the  above  equation,  the  term  Xa,p  can  be  rearranged  as: 


)pi(t)  +h.c,  (2.15) 


With  these  rearrangements,  the  equation  of  motion  of  the  density  matrix  of  the  qubit  can 
be  written  as: 


Hls  ,  Pi(t)  +  V[pi(t)\ ,  where 
*ls  =  E£  At  (w)  Ap(oj) ,  (2.17) 

w  a, (3 

is  called  the  Lamb  shift  since  it  leads  to  the  Lamb  type  renormalization  of  the  unperturbed 
energy  levels  of  the  qubit  device  due  to  coupling  of  the  qubit  to  its  environment.  The  term 
T>[pf(t)\  is  called  the  dissipator  and  takes  the  following  form: 

v\pl(t)\  =  EE  TM  ^Ap(u)pi(t)  At(w)  -  ^|it(u;)A^(w),p/(t)|^  .  (2.18) 

Since  the  term  7ai(g  in  the  dissipator  is  the  Fourier  transform  of  the  the  homogeneous 
reservoir  correlations  functions,  it  is  always  positive  and  can  be  diagonalized.  With  that 
the  evolution  equations  take  the  form  of  the  standard  Lindblad  equation.  We  emphasize 
that  the  physical  assumptions  underlying  the  Lindblad  form  of  the  master  equation  are  the 


dpi  _  1 
dt  ih 
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Born  (weak  coupling),  Markov  (memoryless)  and  rotating  wave  approximation  (fast  system 
dynamics  compared  to  the  relaxation  time).  Lindblad  equation  is  suitable  for  problems 
where  the  above  approximations  are  valid.  In  the  next  section,  we  relax  the  Markovian 
assumption  and  derive  a  new  master  equation. 

Caldeira-Leggett  master  equation 

In  this  section,  we  consider  a  “microscopic”  model  of  the  environment,  which  is  modeled 
by  a  set  of  harmonic  oscillators.  Modeling  environment  by  a  bath  of  harmonic  excitations 
above  a  stable  ground  state  was  first  proposed  by  Caldeira  and  Leggett  [6T] .  We  derive  the 
equation  of  motion  of  the  reduced  density  matrix  of  the  qubit,  following  Ref.  (62J .  In  contrast 
to  the  Lindblad  equation,  here  we  do  not  make  any  Markovian  approximation.  We  assume 
that  system- reservoir  coupling  is  a  linear  function  of  the  system  and  bath  coordinates.  We 
also  assume  that  the  qubit  is  weakly  coupled  to  each  bath  mode.  The  weak  perturbation 
of  any  individual  bath  modes  does  not  necessarily  imply  that  the  dissipative  influence  of 
the  environment  on  the  system  is  weak  as  well.  This  is  because  the  coupling  of  bath  modes 
add  up  and  the  number  of  modes  can  be  very  large,  too.  Here  we  treat  bath  modes  as 
mechanical  oscillators.  For  electrical  environment,  the  complex  impedance  can  be  modeled 
by  simply  replacing  mechanical  oscillators  by  the  LC'-circuits  [63]. 

The  Hamiltonian  of  the  global  system  is  given  by 

H  =  Hq  +  Hr  +  V  ,  where 
Hq  =  tf/2 M  +  U{q) , 

N 

hr  =  y. 

a=l 
N 

V  =  la  q  Xg  •  (2.19) 

OL=l 

Here  Hq  and  Hr  are  the  qubit  system  and  reservoir  Hamiltonians,  respectively,  and  V  is 
the  interaction  Hamiltonian  between  the  qubit  and  its  environment.  It  is  also  necessary  to 
include  a  counter  term  in  the  form  Hc  =  q 2  'Ya/‘^maUJa  in  the  bare  potential  U(q)  to 


Pc 


2  mr 


■  L  2-2 
+  -  ma  uQ  xa 
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avoid  a  constant  shift  in  the  energy  due  to  coupling  of  a  qubit  device  to  the  bath  of  harmonic 
oscillators.  However,  this  term  by  itself  does  not  contribute  to  the  dynamics  of  the  system, 
hence  we  ignore  this  term  in  the  following  discussion.  Assuming  a  weak  coupling  between  the 
system  and  its  environment,  we  can  treat  the  system-reservoir  interaction  perturbatively. 


This  allows  us  to  solve  the  Von-Neumann  equation  2.6  for  the  total  density  matrix  of  the 
entire  system  iteratively  in  powers  of  the  coupling  parameter.  The  solution  can  be  written 
in  the  following  form: 


„  „  i  r*  r  ~ 

W!(t)  -  w>(o)  =  -  J  [VKO,  W,{ 0) 


dt'  +  — 


1 


t  rt' 


iKy 


o  j  o 


V/CO,  VKO-WH*" 


dt'  dt"  . 
(2.20) 


Finally,  expanding  it  to  all  orders  of  the  coupling  parameter,  we  obtain 


w^t)  -  Wi{ 0)  =  Erl  dt'. 


dtr‘ 


n>  1 


[ih)r 


w'), . ,  vdn.wd o) 


(2.21) 


where  T  is  a  time  ordering  operator.  Considering  only  first  two  terms  in  the  above  expan¬ 
sion,  we  have 


wT{t)  =  w>( 0)  +  ^  J*  [ rn 0,  Wi(0) 


t  ft 


'o  Jo 


vho,  W),wm 


dt '  dt" 


(2.22) 


In  order  to  obtain  the  equation  of  motion  of  the  reduced  density  matrix  of  the  qubit,  we 
trace  the  environment  out  followed  by  differentiation  with  respect  to  time.  This  leads  to 


dpi  _  1 
dt  ih  Tir 


VT(t),  W}(0) 


~kL  TrR 


Vi{t'),Wj{ 0) 


dt' 


(2.23) 


It  is  reasonable  to  approximate  that  the  initial  state  of  the  qubit  is  decoupled  from  that  of 
the  reservoir.  Therefore,  the  density  matrix  of  the  entire  system  at  t  =  0  can  be  factored 
into: 

W}(0)  =  W? (0)  ®  ,5/(0) .  (2.24) 

Tracing  the  reservoir  out  and  rearranging  terms  in  |2.22  one  can  immediately  obtain  an 


expression  for  p/(0).  Using  this  expression  for  pj( 0)  along  with  2.24  and  keeping  terms  only 
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up  to  second  order  in  the  coupling  parameter,  Eq.  2.23  reduces  to  the  following  form: 

pt 


dPl  _  1  rp, 

dt  ih TrR 

kl,TlR 


l 


-WoT- 


Vi(t),W?(  0)®p!{t) 

Vi(t ),Wf(0)®Trfl  VKO.P/W 


M*),  [Vf(O,Wf(0)®P/WjJ  ^'+ 

dt' .  (2.25) 


Eq.  |2.25  is  the  master  equation  written  in  the  interaction  picture.  We  can  transform  this 
back  into  the  Schrodinger  picture  by 


dp  1 
dt  ih 


Hq(t),p(t) 


(  -iHqt  \  dpi 

+expHH^exp 


iHqt 

h 


(2.26) 


Assuming  an  initial  state  of  the  reservoir  to  be  in  equilibrium,  that  is,  [WR(0),  Hb]  =  0 
and  (xa(0))  =  0,  which  implies  Tr/j  [V  W^(0)]  =  0,  the  master  equation  in  the  Schrodinger 
picture  can  be  simplified  to  the  following  form: 

~^  =  Hq(t),p(t)  +£[p(t)],  where 
£\P(t)\  =  ~\JQ  ~t)  p{t)dt'  +  ^ K(t' ,  t)  q  p(t)  q{t'  —  t)  dt' 


1 


1 


+  -  /  K(t,t')  q(t'  -  t)  p(t)  qdt'  -  -  /  K (t' ,  t)  p(t)  q(t'  -  t)  q  dt'  (2.27) 
h  Jo  h  Jq 


This  can  also  be  written  in  a  more  compact  form  as 


17  r  r 

=  dt' rii(t  -  t')\q , \q(t'-t),  p(t) 


+  li. 


dt' 


'm(t  - 0  - 1),  /5(t)| 


(2.28) 


Here  the  position  operator  q(t)  is  the  Heisenberg  operator.  In  Eq.  2.28  the  time-dependent 


kernel  K(t,  t')  that  appeared  in  2.27  has  been  split  into  two  parts  and  772 ,  which  are  often 
referred  to  as  quantum  noise  and  dissipation  kernel,  respectively  pJ5|.  These  quantities  can 
be  calculated  analytically  for  the  harmonic  heat  bath: 

K\t',t)  =  K(t,t')  =  K(t  —  t' )  =  771  (t  —  t')  —  irj2(t  —  t') 

N  ,  ,  ^  ,  /-00 

J(w)  [1  +  2 IV  (w)]  cos (ui{t  —  t'))  duj  , 
J[uj)sm.(uj{t  —  t'))  du  .  (2.29) 


Pl{t-t') 

2h  £  7“  ( 

a=l 

{  Xa(t) 

.  N 

m  {t-t') 

=  2h  S  ( 

a=l 

xa(t) 
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Here  the  position  operator  of  the  bath  mode  xa(t)  is  the  Heisenberg  operator.  The  Planck’s 
function  N(u>)  =  l/[exp (hu/3)  —  1]  defines  an  average  excitation  number  of  environment 
modes  with  frequency  oj  and  /3  =  .  where  kB  is  the  Boltzmann’s  constant  and  T  is  the 

temperature  of  the  heath  bath.  Here  J(u )  is  the  spectral  density  of  the  reservoir,  which  is 
given  by 


JM  =  2Z 


a=l 


^S(uj-uja) 
2  fflyi  OJ  /  , 


(2.30) 


It  may  also  be  useful  to  write  the  time  dependent  kernel  K(t,t')  in  the  following  form: 


K(t,  t’)  =  K(t  -t’)=  /  J{u)  +  {1  +  N(u)} 


du 


(2.31) 


Jo  L  J 

From  this  wee  see  that  the  kernel  K(t  —  t')  is  simply  a  Fourier  transform  of  the  spectral 

density  J(uj)  weighted  by  an  average  number  of  bath  excitations  for  emission  and  absorption. 
For  Ohmic  heat  bath,  the  spectral  density  is  given  by 


J(w)  =  2  M  7  uj  exp(— u/uis) , 


(2.32) 


where  7  is  frequency-independent  damping  rate  and  us  is  the  cut-off  frequency.  In  order  to 
obtain  a  memoryless  heat  bath,  that  is,  K(t  —  t')  — >  S(t.  —  t')  (approximately),  the  time-scale 
at  which  the  bath  changes  appreciably  |t  —  t'\  must  be  much  shorter  than  the  time-scale 
at  which  the  qubit  system  changes  appreciably.  In  other  words,  |t  —  t1  (7  <<  1.  Since  the 
maximum  frequency  contribution  of  the  Ohmic  bath  is  ujs ,  we  obtain  a  memoryless  bath  if 
ujs  »  7. 

We  also  note  that  for  the  Ohmic  heat  bath  and  in  the  limit  of  high  temperature  ksT  » 
hcoo,  Caldiera-Leggett  master  equation  for  the  quantum  Brownian  motion  (of  a  particle  in 
a  harmonic  potential)  reduces  to  that  of  Lindblad  form,  with  Lindblad  operator  given  by: 


,  y/AMkBT  .  ,  . 

A  =  - - - q  +  % 

n 


1 


4  MkBT 


P 


(2.33) 


2.5  Discussion 

In  this  chapter,  we  considered  different  sources  of  decoherence  that  affect  superconducting 
and  semiconducting  qubits.  Some  of  the  mechanisms  for  decoherence  in  superconducting 
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circuits  and  quantum  dots  are  external  in  nature,  while  the  others  are  purely  intrinsic  and 
specific  to  the  qubit  device  in  question.  We  also  discussed  theoretical  methods  used  to 
study  the  dynamics  of  an  open  quantum  system.  In  the  following  chapters,  we  use  methods 
discussed  in  this  chapter  to  study  the  dynamics  of  a  qubit  coupled  to  its  environment. 
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Chapter  3 

Effect  of  the  Ohmic  environment 
on  optimally  controlled  flux-biased 
phase  qubit 

3.1  Introduction 

Superconducting  circuits  containing  Josephson  junctions  are  promising  candidates  for  scal¬ 
able  quantum  information  processing  [6d] .  However,  small  separations  between  successive 
quantum  energy  states  in  these  circuits  [65,  66]  do  not  permit  selective  manipulation  of  a 
qubit  in  a  two  dimensional  subspace  and  results  in  a  dynamical  leakage  of  the  quantum 
state  to  a  broader  Hilbert  space  of  the  circuit.  To  reduce  this  leakage,  Motzoi  et  al.  m 
proposed  Derivative  Removal  by  Adiabatic  Gate  (DRAG)  method,  which  reduces  the  gate 
error  to  ICR5  for  an  experimentally  optimal  gate  time  of  6  ns.  This  error  is  well  below  the 
required  error  threshold  of  ICR3  for  fault  tolerant  quantum  computation  [68]. 

In  addition  to  the  dynamic  leakage,  any  realistic  model  of  a  qubit  must  also  address 
coupling  of  the  qubit  to  its  environment,  which  leads  to  further  destruction  of  qubit  states. 
Several  efforts  have  already  been  made  towards  the  study  of  accurate  control  of  a  qubit 
system  [mi  EDi  [nj.  However,  the  effect  of  environment  on  optimally  controlled  qubit  has 
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only  been  studied  in  a  phenomenological  model  EH,  which  leads  to  the  evolution  of  density 
matrix  of  the  qubit  in  the  Lindblad  form  m- 

In  this  chapter,  we  resort  to  a  microscopic  approach  to  the  modeling  of  the  environment. 
We  employ  the  Caldeira-Leggett  model  of  the  system-environment  coupling,  introduced  in 
chapter  [2j  to  describe  the  coupling  of  a  flux-biased  phase  qubit,  which  is  driven  by  the 
DRAG  pulses,  to  the  electromagnetic  environment.  Numerically  solving  equation  of  motion 
for  the  qubit  density  matrix,  we  study  the  dependence  of  the  gate  error  on  temperature, 
gate  time  and  environmental  coupling  strength.  Although  numerous  potential  sources  of 
decoherence  in  phase  qubits  have  been  identified  experimentally  mmm,  in  this  chapter, 
we  focus  on  decoherence  due  to  the  Gaussian  noise  from  the  surrounding  circuitry,  which  is 
introduced  to  the  qubit  system  within  the  Caldeira-Leggett  model.  The  study  of  the  effect 
of  ubiquitous  low  frequency  1  //  noise  on  the  gate  error  is  out  of  the  scope  of  this  chapter. 

We  specifically  study  the  role  of  dissipation  in  the  gate  error  during  the  NOT  gate  oper¬ 
ation.  We  find  that  for  phase  qubits  with  relaxation  time  T\  ~  700  ns  m,  two-quadrature 
DRAG  pulses  proposed  in  Ref.  m  result  in  the  gate  error  exceeding  7  x  10  3,  which  is  too 
high  for  fault  tolerant  quantum  computation.  We  then  address  the  limitation  posed  by  the 
environmental  coupling  on  two-quadrature  pulses.  Here  we  find  that  for  optimal  DRAG 
pulses  m,  the  coupling  to  environment  must  be  reduced  nearly  by  a  factor  of  6  to  suppress 
the  gate  error  below  the  required  threshold.  We  also  investigate  the  gate  error  for  simple 
pulse  shaping  where  the  pulse  amplitude  of  the  first  quadrature  vary  smoothly  according  to 
a  Gaussian-shaped  function  while  the  amplitude  of  the  second  quadrature  is  proportional 
to  the  derivative  of  the  first.  In  this  case,  however,  the  microwave  drive  frequency  is  held 
constant.  For  this  pulse  shape,  we  find  that  the  gate  error  reduces  to  10~3  for  a  gate  time 
~  7  ns  when  the  coupling  to  environment  is  reduced  by  an  approximate  factor  of  10  from 
the  coupling  in  currently  used  phase  qubits. 
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3.2  Model 


A  flux-biased  phase  qubit  consists  of  a  Josephson  junction  ( JJ)  embedded  in  a  superconduct¬ 
ing  loop.  The  schematic  of  a  flux-biased  phase  qubit  is  show  in  Fig.  |3.1[  Finite  resistance  of 


the  JJ  (see  Fig.  3.1  (b))  results  in  dissipation  processes  in  the  qubit  and  can  be  accounted 
for  by  the  Caldeira-Leggett  model,  discussed  in  the  previous  chapter.  The  full  Hamiltonian 
of  the  qubit  and  the  environment  is 


H  =  Hq  +  P(t)  +  Hr  +  V. 


(3.1) 


The  Hamiltonian  of  a  flux-biased  phase  qubit  Hq  is  written  in  terms  of  operators  Q  and  5, 
the  charge  and  phase  difference  of  the  JJ  respectively: 


Q 


1 


4>n 


2  C  2  L 


2t r 


/q^O 

2vr 


cos  5 , 


(3.2) 


where  L  ( C )  is  the  loop  inductance  (junction  capacitance),  is  the  external  magnetic 
flux  applied  to  the  phase  qubit,  Iq  is  the  critical  current  of  the  JJ,  and  <f>o  =  h/2e  is  the 
flux  quantum.  The  qubit  is  capacitively  coupled  to  a  microwave  current  source,  which  is 
used  to  induce  coherent  transitions  between  the  qubit  states.  This  coupling  introduces  a 
time-dependent  part  in  the  Hamiltonian: 


P(t) 


m 

2v r 


(3.3) 


Here  I(t)  =  Ix(t)  cosw^f  +  Iy(t)  sincu^f  is  the  microwave  current  with  frequency  Ud-  The 
coupling  between  the  qubit  system  and  the  reservoir  is  bilinear  in  the  JJ  phase  5  and 
oscillator  displacements  xa: 

V  =  \2'yaxaq,  q  =  5 - ^ -  ,  (3.4) 

where  parameters  7 a  determine  the  coupling  strength  between  the  qubit  and  reservoir  mode 
a.  Our  goal  is  to  describe  the  time  evolution  of  the  qubit  density  matrix  p{t).  The  qubit  is 
initially  prepared  in  a  pure  state,  corresponding  to  the  density  matrix  p(0).  Assuming  that 
the  environment  is  in  a  thermal  equilibrium  at  temperature  T,  the  master  equation  for  p(t) 
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Figure  3.1:  (a)  Schematic  of  a  flux-biased  phase  qubit,  and  (b)  Resistively  and  capacitively 
shunted  junction  model  (RCSJ) 


takes  the  following  form: 


dp{t) 


1  - 

ih  - 


dt  Hq(t),  p(t)  +£[p(t)] 


(3.5) 


where  the  Caldeira-Leggett  (CL)  dissipative  term  C[p(t)\  was  derived  in  the  previous  chap¬ 


ter.  Although  the  CL  master  equation  2.27,  is  local  in  time,  it  contains  time-dependent 


coefficients,  which  capture  memory  effects  of  the  noise  due  to  the  heat  bath.  Here  we 
consider  the  spectral  density  for  the  Ohmic  environment,  which  is  given  by 


C 

J(w)  =  £,^h2Lu0uj  exp(-u/uis)  , 


(3.6) 


where  £  is  a  dimensionless  coupling  parameter,  hoj o  is  the  energy  difference  between  the 
qubit  states  and  us  is  a  cut-off  frequency  that  exceeds  all  other  frequency  scales  of  the 
qubit  system. 

In  typical  experiments  with  phase  qubits  mm,  the  “potential”  part  of  Hq  ini 


3.2 


has 


one  deep  minimum  and  another  very  shallow  minimum  that  disappears  at  the  critical  flux 
<LC.  External  flux  $5  is  chosen  in  such  a  way  that  only  a  few  levels  are  localized  in  the 
shallow  well,  but  these  levels  are  still  separated  from  levels  localized  in  the  deep  well  by 
impenetrable  barrier  f  As  a  result,  we  truncate  the  qubit  Hamiltonian,  Eqs. 


3.2 


and 


3.3 


lrOie  levels  in  the  deep  well  can  also  be  accounted  in  the  present  model,  however,  our  numerical  results 
indicate  that  the  gate  error  does  not  change  significantly  if  those  levels  are  also  included  in  the  calculation 
for  chosen  values  of  parameters. 
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to  three  localized  levels  and  obtain  the  following  Hamiltonian: 

2 

Hq(t)  =  h  ^  [ajj-illj  +  aXjd-j-  +  a*Xjaj 


+  Hn 


(3.7) 


3= 1 

where  IF,  =  \j)(j\  is  the  projector  for  the  jth  level,  =  \j)(j  —  1 1  is  the  raising  operator, 
a  =  (Ix  —  Hy)eluldt/2  is  the  amplitude  of  microwave  drive,  A j  =  ^o(j\S\j  —  1) / 27 xh  is  the 
matrix  element  of  the  phase  operator,  ojj  =  (Ej+i  —  £o)/fr,  £j  is  an  energy  eigenvalue  of 
time-independent  Hamiltonian  Hq  and  Hnr  contains  non-resonant  terms.  In  this  three-level 
model,  the  lower  two  energy  levels  comprise  qubit  states  while  the  third  level  accounts  for 
a  leakage  level. 


3.3  Gate  error  and  DRAG  method 


In  order  to  quantify  the  error  during  gate  operation,  we  use  gate  fidelity  averaged  over  two 
initial  input  states  in  a  two  dimensional  Hilbert  space,  similar  to  one  defined  in  Ref.  m- 


F*  =  l 


Tr  U ideal  p\  )  u}deal  Pj  ( tg  ) 


(3.8) 


3= 1 

Here  Uideal  represents  an  ideal  evolution,  pj(t)  is  an  actual  density  matrix  of  the  qubit 
system  with  Pj(0)  =  Pj°\  and  pp  represents  two  initial  axial  states  in  a  Bloch  sphere.  The 
gate  error  E  is  defined  as  E  =  1  —  Fg. 

A  simple  approach  to  minimize  leakage  of  quantum  information  from  the  qubit  subspace 
is  to  use  a  single-quadrature  Gaussian  envelope  pulse  given  by 

4 it)  =  4(t)  =  ie-442)2/2-2—  B ,  Iy(t)  =  0  ,  (3.9) 


where  tg  is  a  gate  time  and  a  =  tg/2.  For  a  NOT  gate  operation,  which  we  choose  to  focus 
on  without  any  loss  of  generality,  constant  B  is  chosen  so  that  the  Gaussian  pulse  starts 
and  finishes  off  at  zero  and  A  is  defined  by 

In(t)  dt  =  7r.  (3.10) 

This  pulse  shape  results  in  a  large  gate  error,  due  to  population  of  the  third  level,  even  for 
reasonably  short  pulses. 
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The  DRAG  method  reduces  the  gate  error  to  order  of  10  5  for  a  gate  time  of  6  ns  EH  by 
using  two  quadratures  and  time-dependent  detuning  Ai(f)  =  ujo  —  cod  =  (A2  —  4)/2(f) / 4<5an, 
where  the  anharmonicity  parameter  5an  =  —  2wo,  and  A  measures  relative  strength 

of  0  — >  1  and  1  — >  2  transitions,  that  is,  A  =  A2/A1.  The  DRAG  protocol  uses  second 
quadrature  pulse  to  suppress  the  population  of  the  third  level  generated  by  single  quadrature 
drive.  This  is  evident  from  the  Hamiltonian  written  in  the  adiabatic  frame  (see  below).  We 
note  that  the  laboratory  frame  is  more  suitable  for  the  solution  of  the  reduced  density 
matrix  of  the  qubit  coupled  to  its  environment.  Therefore,  we  preserve  the  form  of  the 
quadrature  amplitudes  as  in  Ref.  m 


4 


(A2  -4 )J3 
8  S2an 


(3.11) 


and  obtain  the  following  equation  for  the  microwave  driving  frequency  for  the  Hamiltonian 


Eq.  3.7  in  the  laboratory  frame 


tujd(t)  +  ud{t)  =  u0  -  Ai (t)  ,  wd(0)  =  w0  . 


(3.12) 


Although  the  DRAG  correction  is  successful  in  reducing  the  gate  error  below  the  required 
threshold,  practical  implementation  may  not  be  feasible  due  to  stringent  requirement  to  vary 
microwave  frequency.  For  this  reason,  we  also  consider  two-quadrature  pulses  with  fixed 


driving  frequency  ujd  =  wq  2 .  We  transform  the  Hamiltonian  Eq.  3.7  to  a  frame  rotating 


with  frequency  ud  with  respect  to  the  laboratory  frame  and  obtain: 

2 

|  Ajiij  + 


HR  =  h 
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4(^)  ,  ~x  |  4(^)  \  -y 

~A3aj-hj  +  9 


(3.13) 


3= 1 


where  A2  =  5an  +  2  Ai,  and  for  ojd  =  wo,  the  detuning  Ai  =  0.  We  introduce  operators 

Oj,k  =  \k)(j\  +  and  &lk  =  -  Aj){kV 

To  analyze  the  dynamics  of  rotating  frame  Hamiltonian  HR,  it  is  convenient  to  perform 


an 


adiabatic  transformation  [6?]  D(t)  =  exp  — ilx(t )  <Jq  1  +  A cr\ /2 5t 


,  which  pre¬ 


serves  the  form  of  the  gate,  if  Ix(t)  starts  and  finishes  off  at  zero.  This  condition  is  satisfied 


2 We  also  made  constant  detuning  of  the  driving  frequency  from  u>o,  but  did  not  see  any  improvement 

compared  to  toa  =  u>o  case. 
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by  our  choice  of  Ix{t )  (see  Eq.  3.9).  The  parameter  a  appearing  in  I?  is  a  dimensionless 


scaling  parameter.  After  performing  the  transformation,  the  Hamiltonian,  to  first  order  in 
Ix / 5am  takes  the  following  form: 
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We  then  require  resonant  condition  for  the  microwave  n— pulse  in  the  qubit  subspace  and 
also  eliminate  the  imaginary  inertial  term  from  the  subspace,  that  is,  require  Hi  and  <7q  1 


terms  in  Eq.  3.14  to  vanish,  to  obtain: 


A2  _  / , \  t  /  \  t  / , \  —ctl-xit) 

a  =  ~r  >  Ix(t)  =  I-irif)  1  Iy(t)  =  7 


(3.15) 


where  In(t)  is  defined  by  Eqs.  3.9  The  contributions  to  the  gate  error  due  to  transitions 
to  the  third  level  come  from  the  second  and  third  lines  of  Eq. 


3.14 


except  for  Ht  term, 


which  is  not  directly  responsible  for  the  gate  error.  Using  the  above  expression  for  Iy(t) 
and  Eq.  3.10[  we  estimate  the  magnitude  of  these  terms  as 
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These  estimates  show  that  the  error  due  to  <7,  2  and  CTq  2  terms  are  comparable  and  results 
in  the  leading  contribution  to  the  gate  error.  In  the  case  of  a  time-dependent  detuning,  the 
choice  of  pulses  is  such  that  it  eliminates  the  error  associated  with  cr\  2  term,  and  rescaling 


of  the  pulse  intensity  (I3  term  in  Eq.  3.11)  removes  the  contribution  to  the  gate  error  due 
to  <7q2  term.  This  elimination  of  ~  l/(8ant2)  explains  high  effectiveness  of  variable  driving 
frequency  DRAG  pulses.  For  fixed  frequency  pulses,  the  pulse  rescaling  only  marginally 
reduces  the  gate  error. 
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Figure  3.2:  Gate  error  vs.  gate  time  in  log-normal  scale  with  (thick  lines)  and  without 
(thin  lines)  dissipation  for  a  single  quadrature  Gaussian  (a  =  0.5  tg)  pulse  (solid  black), 
the  Gaussian  ( a  =  0.5  tg)  pulse  with  first  order  DRAG  correction  and  dynamical  detuning 
(dashed-dot  blue),  and  the  Gaussian  (a  =  0.5  tg)  pulse  with  fixed  driving  frequency  uid  =  wo 
and  a  =  0.5  (dashed  red),  all  in  the  laboratory  frame.  For  the  dissipative  case,  temperature 
T  =  0.1  hu o,  the  cutoff  frequency  cos  =  IOcjq  and  the  coupling  parameter  £  =  2. 


3.4  Results 


Qubit  parameters  used  below  in  our  simulation  are  typical  of  phase  qubits:  C  =  1  pF, 
Iq  =  1.5  /i  A,  /3l  =  2ttI0L/^0  =  3.2  and  <!>{,  =  0.955  d>c  ,  where  is  a  critical  flux.  Numerical 
simulation  indicates  that  small  variations  of  qubit  parameters  do  not  incur  any  noticeable 
change  in  the  gate  error  as  long  as  there  are  at  least  three  energy  levels  in  the  shallow 
well  of  the  potential.  For  these  experimental  parameters,  we  numerically  solve  the  time- 


independent  Schrodinger’s  equation  with  the  Hamiltonian  given  by  Eq.  3.2  From  this 
simulation,  we  obtain  the  following  numerical  values  (rounded  up  to  two  decimal  places): 
ujq/2tt  =  6.2  GHz,  A  =  1.42  and  <5an/27r  =  —380  MHz. 


In  Fig.  3.2 ,  we  plot  the  gate  error  for  the  DRAG  pulses  with  and  without  time-dependent 
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detuning  for  an  ideal  phase  qubit  without  environment.  We  find  that  pulses  with  two 
quadratures  and  fixed  driving  frequency  (thin  dashed  red)  perform  much  better  than  single 
quadrature  Gaussian  pulses  (thin  solid  black) ,  but  are  not  as  effective  as  pulses  with  double 
quadratures  and  time-dependent  driving  frequency  (thin  dashed-dot  blue). 

We  verify  numerically  that  the  fixed  frequency  DRAG  pulses  give  the  minimal  gate  error 


for  the  choice  of  parameter  a  according  to  Eq.  3.15  As  shown  in  Fig.  |3.3[a),  minimum 
value  of  the  error  occurs  at  around  a  =  0.5  for  different  gate  times,  cf.  dashed  blue  curve 


for  tgU o  =  250  and  solid  black  curve  for  tguo  =  350.  This  result  is  consistent  with  Eq.  3.15 


since  for  the  phase  qubit  A  =  1.42,  which  implies  a  =  0.5.  For  transmon  qubits,  discussed 
in  Ref.  [T9] .  a  =  0.4  owing  to  different  value  of  A. 

In  order  to  study  the  effect  of  dissipation  on  the  DRAG  pulses,  we  integrate  the  master 


equation  4.4  numerically  using  the  fourth  and  fifth  order  Runge-Kutta  method.  First, 
we  consider  relaxation  of  the  qubit  from  the  first  excited  state  to  the  ground  state  in  the 
absence  of  a  microwave  drive,  shown  in  Fig.|3.3[b).  For  this  simulation,  we  choose  the  cut-off 
frequency  us  =  10 wo  (throughout  this  paper),  temperature  T  =  0.1/kJoj  and  the  coupling 
parameter  £  =  2  so  that  the  relaxation  time  corresponds  to  experimentally  observed  decay 
time  of  Ti  ~  700  ns  for  phase  qubits  m-  We  note  that  the  spontaneous  relaxation  rate  of 
the  first  excited  state  can  also  be  evaluated  from  the  master  equation  4.4  as 


1 


r  =  —  =  2irhujQ^\q01\2,  q0i  =  (0|g|l). 


T\ 


(3.16) 


For  the  above  choice  of  dimensionless  coupling  parameter  £,  temperature  T  and  the  cut-off 
frequency  uis,  we  study  the  effect  of  dissipation  on  two-quadrature  pulses. 


In  Fig.  3.2  we  observe  a  non-monotonic  behavior  of  the  gate  error  with  gate  time  for 
pulses  with  the  DRAG  corrections.  We  find  that  for  shorter  gate  times,  two-quadrature 
pulses  with  time-dependent  driving  frequency  are  weakly  affected  by  dissipation  (thick 
dashed-dot  blue).  However,  for  longer  gate  times,  dissipation  has  a  substantial  effect  on 
two-quadrature  pulses.  For  instance,  for  a  gate  time  of  u>o tg  =  250  (tg  ~  6  ns),  the  gate 
error  increases  from  10~5  to  higher  order  of  10~ 3  for  two-quadrature  pulses  with  dynamical 
detuning,  when  dissipation  is  taken  into  account.  This  increase  in  the  gate  error  is  due  to 
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Figure  3.3:  (a)  Gate  error  vs.  alpha  for  gate  times  uJotg  =  250  (dashed  blue)  and  cootg  =  350 
(solid  black).  The  phase  qubit  is  driven  by  two-quadrature  pulses  with  driving  frequency 
uid  =  ojq  hi  the  absence  of  environment,  (b)  Probability  vs.  time  for  temperature  T  = 
0.1  hwo,  coupling  parameter  £  =  2  and  cut-off  frequency  us  =  10  uiq.  The  microwave  pulse 
is  turned  off  and  the  qubit  is  initially  prepared  in  |1)  state  (solid  black),  which  relaxes  to 
|0)  state  (dashed  blue)  due  to  dissipation. 


relaxation  of  the  qubit  from  the  excited  state  to  the  ground  state,  which  becomes  prominent 
for  longer  gate  times.  For  comparison,  we  plot  the  gate  error  for  three  different  types  of 
pulses:  single-quadrature  Gaussian  pulse  (thick  solid  black),  the  Gaussian  pulse  with  first 
order  DRAG  correction  and  time-dependent  driving  frequency  (thick  dashed-dot  blue)  and 
the  Gaussian  pulse  with  two  quadratures  and  fixed  driving  frequency  (thick  dashed  red). 
One  can  conclude  from  these  plots  that  the  performance  of  two-quadrature  pulses  without 
detuning  is  comparable  to  the  DRAG  pulses  with  dynamical  detuning  when  dissipation  is 
included. 

Next,  we  study  the  effect  of  environmental  coupling  strength  on  the  gate  error.  In 


Fig.  3.4,  we  plot  the  gate  error  for  different  coupling  parameters  £  for  the  phase  qubit 
driven  by  two-quadrature  pulses  with  dynamical  detuning.  In  this  simulation,  we  consider 
the  temperature  T  =  O.lfojjQ,  and  coupling  parameters:  £  =  0  (thin  dashed-dot  blue), 
£  =  0.1  (thick  solid  green  with  circles),  £  =  0.3  (thick  solid  red)  and  £  =  2  (thick  dashed- 
dot  blue).  At  £  =  0  the  gate  error  originates  entirely  due  to  microwave-induced  leakage 
of  the  qubit  state  from  the  lowest  two  level  subspace.  The  environment-induced  transition 


rates  increase  with  increase  in  the  environmental  coupling  strength  £  (see  Eq.  3.16).  As 
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Figure  3.4:  Gate  error  vs.  gate  time  in  log-normal  scale  for  the  Gaussian  (a  =  0.5  tg) 
pulse  with  variable  frequency  DRAG  correction  for  temperature  T  =  0.1  huo  and  cut-off 
frequency  us  =  10  uq.  The  environmental  coupling  parameters:  £  =  0  (thin  dashed-dot 
blue),  £  =  0.1  (thick  solid  green  with  circles),  £  =  0.3  (thick  solid  red)  and  £  =  2  (thick 
dashed-dot  blue). 


a  result,  the  gate  error  also  increases,  which  is  corroborated  by  Fig.  3.4  One  can  infer 
from  the  plot  that  two-quadrature  pulses  with  time-dependent  driving  frequency  suppress 
the  gate  error  to  10~3  for  £  =  0.3  and  gate  time  iv^tg  ~  200  (tg  ~  5.5  ns).  This  indicates 
that  an  increase  in  the  relaxation  time  nearly  by  a  factor  of  6  from  the  currently  observed 
value  is  necessary  to  suppress  the  error  below  the  threshold  for  fault  tolerant  computation. 
We  further  analyze  the  effect  of  environmental  coupling  on  fixed  frequency  two-quadrature 
pulses  for  a  range  of  gate  times.  For  this  case,  gate  errors  for  different  values  of  £  are  plotted 
in  Fig.  |3.5|,  where  the  temperature  is  the  same  as  above  and  coupling  parameters  are:  £  =  0 
(thin  dashed  red),  £  =  0.1  (thick  solid  green  with  circles),  £  =  0.2  (thick  solid  blue),  £  =  0.5 
(thick  solid  black  with  triangles)  and  £  =  2  (thick  dashed  red).  These  plots  indicate  that 
the  DRAG  pulses  with  fixed  driving  frequency  can  effectively  suppress  the  gate  error  if  the 
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Figure  3.5:  Gate  error  vs.  gate  time  in  log-normal  scale  for  the  Gaussian  ( a  =  0.5  tg)  pulse 
with  fixed  frequency  DRAG  correction  for  temperature  T  =  0.1  huo  and  cut-off  frequency 
ojs  =  lOwo-  The  environmental  coupling  parameters:  £  =  0  (thin  dashed  red),  £  =  0.1 
(thick  solid  green  with  circles),  £  =  0.2  (thick  solid  blue),  £  =  0.5  (thick  solid  black  with 
triangles)  and  £  =  2  (thick  dashed  red). 


environmental  coupling  strength  is  weakened  and  gate  times  are  slightly  longer  than  6  ns. 
More  specifically,  for  £  =  0.2  and  a  gate  time  of  uotg  ~  300  (tg  ~  7  ns),  the  gate  error  is 
close  to  10~3.  Therefore,  we  conclude  that  the  relaxation  time  must  be  nearly  a  factor  of  10 
longer  than  the  currently  observed  value  to  attain  the  threshold  of  the  gate  error  for  fixed 
frequency  DRAG  pulses.  This  is  a  much  better  improvement  compared  to  single-quadrature 
pulses  for  which  the  gate  error  never  reduces  to  the  threshold  for  a  reasonable  choice  of  gate 


times  even  in  an  ideal  case,  that  is,  £  =  0,  as  shown  in  Fig.  3.2  Finally,  we  investigate  the 


effect  of  the  temperature  on  the  gate  error.  In  Fig.  3.6,  we  plot  the  gate  error  normalized 


to  the  error  at  zero  temperature  for  two  different  gate  times:  uq tg  =  150  (dashed-dot  blue) 
and  uotg  =  350  (dashed  black).  The  plot  shows  a  monotonic  growth  of  the  gate  error  as 
the  temperature  increases  due  to  enhancement  in  the  relaxation  rate.  We  compare  results 
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Temperature  (T)  [hcoQ/27r] 

Figure  3.6:  Normalized  gate  error  vs.  temperature  for  the  Gaussian  (a  =  0.5  tg)  pulse 
with  varying  frequency  DRAG  correction.  In  numerical  simulation,  environmental  coupling 
parameter  £  =  2  and  cut-off  frequency  cos  =  10  wo  are  used  for  gate  times  uiotg  =  150 
(dashed-dot  blue)  and  uiQtg  =  350  (dashed  black).  Analytical  rate  equation  estimation  of 
the  normalized  gate  error  (solid  red). 


of  the  numerical  solution  of  the  master  equation  4.4,  and  a  simple  picture  of  the  error 
due  to  coupling  to  the  environment  in  terms  of  the  “Fermi-Golden  rule”  transition  rates. 
Considering  the  environment  at  zero  temperature  and  assuming  that  the  contribution  to  the 
error  E  from  the  environment  is  small,  E  <C  1,  we  can  evaluate  the  error  as  the  probability 
of  an  excitation  of  a  reservoir  mode  during  the  qubit  operation,  which  happens  with  rate  T : 
E(T  =  0)  =  Ttgpu{t),  where  pu(t)  =  /0<9  pn(t)dt/tg  is  the  time-average  of  the  probability 
of  the  qubit  being  in  the  first  excited  state.  At  finite  temperature,  processes  with  excitation 
of  environment  happen  with  rate  T(T)  =  T[1  +  Al(cu)].  In  addition,  the  qubit  can  absorb  an 
excitation  from  the  environment  with  rate  TN(u).  We  combine  the  qubit  excitations  from 
the  ground  to  first  excited  state  and  the  first  to  second  excited  state  with  the  relaxation 
from  the  first  to  ground  state  and  obtain  the  following  estimate  for  the  gate  error  due  to 
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coupling  to  the  environment: 


~  [{1  +  ^(wo)}  +  A2iV(u;i  -  wo)]  pn(i)  +  N(uo)poo{t) .  (3.17) 

Assuming  the  average  occupation  of  the  ground  and  first  excited  states  to  be  ~  1/2,  and 
for  a  weak  anharmonicity  of  the  qubit  system  |<5an|  <C  wq,  the  gate  error  reduces  to 


E(T) 

m 


1  +  41V ( wo)  ■ 


(3.18) 


The  estimated  normalized  gate  error  (solid  red)  is  plotted  in  Fig.  3.6  together  with  the 
gate  error  obtained  from  the  numerical  simulation.  The  rate  equation  estimate  of  the  error 
is  fairly  close  to  the  error  obtained  from  direct  numerical  simulation  for  a  longer  gate  time 
(dashed  black).  However,  for  a  shorter  gate  time  (dashed-dot  blue),  the  estimated  error 
deviates  from  the  exact  numerical  simulation  considerably  suggesting  that  the  rate  equation 
description  may  not  be  valid  for  shorter  gate  times,  as  well  as  for  higher  temperatures. 


3.5  Discussion 

In  this  chapter  we  compared  possible  choices  of  microwave  pulses  for  the  NOT  gate  oper¬ 
ation  in  fluxed-biased  qubits.  Particularly,  we  considered  three  options:  single-quadrature 
pulses  and  two-quadrature  microwave  (control)  pulses  with  both  variable  and  fixed  fre¬ 
quencies.  Two-quadrature  pulses  led  to  significant  suppression  of  the  gate  error  compared 
to  single-quadrature  pulses.  However,  the  presence  of  dissipative  environment  increased 
the  gate  error  even  for  two-quadrature  pulses  significantly  above  the  required  threshold  for 
fault  tolerant  quantum  computation  in  currently  available  phase  qubits.  We  further  inves¬ 
tigated  how  the  environmental  coupling  strength  affects  the  gate  error  and  found  that  an 
improvement  of  the  qubit  relaxation  time  is  crucial  for  effectiveness  of  the  DRAG  pulses. 
We  determined  that  two-quadrature  pulses  with  fixed  driving  frequency  suppressed  the  gate 
error  below  the  required  threshold  for  a  reasonable  gate  time  of  7  ns,  but  for  qubits  with 
the  relaxation  time  ten  times  longer  than  the  currently  observed  relaxation  times.  Simi¬ 
larly,  our  analysis  indicated  that  two-quadrature  pulses  with  dynamical  detuning  can  also 
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effectively  reduce  the  gate  error  below  the  required  threshold  if  the  relaxation  time  is  longer 
by  at  least  a  factor  of  6.  In  comparison  to  fixed  frequency  DRAG  pulses,  this  is  a  moderate 
improvement  over  the  longer  relaxation  time  requirement,  yet  not  impressive  enough  to 
overcome  the  difficulty  associated  with  implementing  control  pulses  with  variable  driving 
frequency.  We  expect  that  in  a  trade-off  between  complicated  driving  frequency  and  longer 
relaxation  times,  the  DRAG  pulses  with  fixed  frequency  are  viable  alternatives  for  reducing 
the  gate  error.  We  emphasize  that  for  single-quadrature  pulses,  reduction  of  the  gate  error 
below  the  error  threshold  of  10~3  is  not  possible  for  reasonable  gate  times,  even  in  an  ideal 
case  without  any  dissipation. 

In  addition,  we  observed  a  monotonic  increase  in  the  gate  error  with  the  temperature, 
which  is  due  to  increase  in  the  relaxation  rate  with  the  temperature.  We  found  that  the 
temperature  dependence  of  the  gate  error  for  longer  pulses  can  be  captured  by  a  simple 
error  estimation  based  on  the  rate  equations.  Nonetheless,  these  estimates  of  the  error  for 
shorter  pulses  differs  from  the  gate  error  obtained  from  direct  numerical  solution  of  the 
reduced  density  matrix.  Therefore,  full  density  matrix  solution  is  necessary  to  calculate  the 
error  for  shorter  gate  times. 
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Chapter  4 

Microwave  photon  detector  based 
on  a  current-biased  Josephson 
junction 

4.1  Introduction 

Quantum  optical  photodetectors  have  contributed  significantly  to  the  development  of  quan¬ 
tum  optics  and  atomic  physics  [80]  and  now  play  an  essential  role  in  optical  quantum  in¬ 
formation  applications  such  as  quantum  computing  and  quantum  key  distribution  [Ell  182] . 
Recently,  circuit  quantum  electrodynamics  (cQED)  has  emerged  as  a  novel  paradigm  for  the 
study  of  radiation-matter  interaction  in  mesoscopic  systems  [|83,  84,  E5|.  Moreover,  cQED 
is  an  attractive  candidate  for  scalable  quantum  computing  and  transmission  of  quantum  in¬ 
formation  [S3  EH-  Following  the  original  proposal,  a  variety  of  cQED  architectures  demon¬ 
strating  strong  coupling  between  single  photons  and  superconducting  integrated  circuits 
have  been  realized  experimentally  [881 189j.  This  work  has  paved  the  way  for  the  develop¬ 
ment  of  a  superconducting  microwave  photon  detector  [90,  .91]  with  possible  applications 
to  quantum  information  processing  and  communication  [92j. 

The  microwave  photon  detector  is  based  on  the  current-biased  Josephson  junction  (JJ): 
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Figure  4.1:  (a)  Schematic  circuit  diagram  of  JJ-based  microwave  photon  detector  coupled 
to  a  resonator,  (b)  Potential  energy  landscape  of  the  detector  when  bias  current  is  close  to 
the  critical  current  of  the  JJ  .  The  junction  is  initialized  in  the  | g)  state  and  upon  absorbing 
an  incident  photon  transitions  to  the  |e)  state,  which  rapidly  tunnels  to  the  continuum. 

the  JJ  is  biased  so  that  absorption  of  a  single  microwave  photon  induces  a  transition  to  the 
voltage  state,  producing  a  large  and  easily  detected  classical  signal.  While  these  detectors 
are  straightforward  to  operate  and  show  potential  for  scalability,  performance  is  degraded 
by  spurious  dark  counts  due  to  quantum  tunneling  events  in  the  absence  of  an  absorbed 
photon;  moreover,  energy  relaxation  within  the  detector  results  in  photon  loss  and  leads  to 
a  reduction  in  measurement  efficiency. 

In  this  work,  we  theoretically  determine  the  quantum  efficiency  of  a  microwave  photon 
detector  based  on  a  current-biased  JJ.  Previous  analysis  [93;  ;9lj  of  this  system  was  focused 
on  a  wave-packet  formulation  of  the  photon  field  in  a  transmission  line  coupled  to  the 
detector.  Here  we  study  the  probability  of  photon  detection  by  a  JJ  coupled  to  a  microwave 
cavity  mode  that  is  loaded  with  a  fixed  number  of  photons  [95].  We  solve  the  equation  of 
motion  for  the  density  matrix  of  the  cavity- JJ  system  to  obtain  detector  efficiency  for 
different  values  of  operation  time,  current  bias,  and  relaxation  time  of  the  junction.  Our 
results  indicate  that  a  JJ  with  decay  time  around  10  ns  can  detect  a  single  microwave 
photon  in  the  cavity  with  an  efficiency  greater  than  80%,  for  readily  achievable  circuit 
parameters.  We  also  find  that  the  detector  efficiency  increases  significantly  with  increasing 
energy  relaxation  time  T\  of  the  junction,  suggesting  that  a  highly  efficient  single  microwave 
photon  detector  is  attainable  for  moderate  improvements  in  junction  quality. 
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4.2  Josephson-j unction  based  photon  detector 


The  circuit  diagram  of  the  JJ-cavity  system  is  shown  in  Fig.  |4.1[a).  The  JJ  is  biased  with  a 
current  I  close  to  the  critical  current  I q.  The  junction  Hamiltonian  can  be  written  in  terms 
of  the  charge  operator  Q  and  the  operator  6  of  the  phase  difference  across  the  JJ  |96|: 

Hjj  =  ^  +  U(5) ,  U(6)  =  (cos  5-  j-j')  .  (4.1) 

Here  C  is  the  junction  capacitance  and  4>o  =  h/2e  is  the  magnetic  flux  quantum.  For 
I  <  70,  the  potential  energy  landscape  U(5)  takes  on  a  “tilted  washboard”  shape,  with  a 
few  discrete  energy  levels  in  shallow  minima  separated  from  the  continuum  by  a  barrier, 
see  Fig.  |4.1|(b).  We  truncate  the  junction  Hamiltonian  to  the  ground  | g)  and  first  excited 
states  |e)  and  obtain  the  following  Hamiltonian  for  the  JJ: 


Hjj 


(4.2) 


where  ne  =  |e)(e|  is  the  projection  operator  to  the  excited  state  and  ojeg  =  (ee  —  eg)/h. 

The  coupling  of  the  cavity  with  the  JJ  is  modeled  by  the  Jaynes- Cummings  (JC)  Hamil¬ 
tonian  m- 


Hjc  = 


1\  -  Kl 

nu>r  I  a  'a  +  -  I  +  niuegIle  +  — (a1  +  aa+) 


(4.3) 


where  u:r  is  the  cavity  resonance  frequency,  P  is  the  vacuum  Rabi  frequency,  and  d^(a), 
<r+((T_)  are  the  creation  (annihilation)  operators  of  the  cavity  and  the  junction,  respectively. 

The  time  evolution  of  the  density  matrix  p(t)  of  the  cavity- JJ  system  coupled  to  its 
environment  is  governed  by  the  following  equation: 


dp(t ) 
dt 


ih  . 


Hjc,  pit)  +  ^j[p(t)\  +  ^n[p{t)]  +  C-T[p{t)] , 


(4.4) 


where  /l7[p(t)]  and  CK{p{t)\  are  superoperators  that  capture  damping  in  the  JJ  and  the 
cavity  at  low  temperatures  T  <C  hojeg,  hur  [98] : 


CK[p{t)\  =  k  (dpd)  -  ]^{a]a,p}  \  , 


(4.5a) 


£y\P(t)]  =  7  (  d-P&+  ~ 


(4.5b) 
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To  account  for  switching  of  the  ground  and  the  first  excited  states  of  the  JJ  to  the 
voltage  state,  we  introduce  the  tunneling  superoperator  T/j/il  /)  !JT?).  1.00. 11011 1 1 02] : 

— 0C  +  PST) 

(aC  +  PT) 

(4.6) 

where  Tej9  are  the  tunneling  rates  from  the  ground  (\g))  and  first  excited  (|e))  states  of  the 
junction.  If  we  approximate  the  potential  in  Eq.|4.1|by  a  cubic  potential,  then  the  tunneling 
rate  of  the  ground  and  the  first  excited  states  of  the  cubic  potential  can  be  computed  by 
WKB  approximation: 

Tj  =  up/2Tr['i32AU/hu;p\j+1/2 /nj/2  exp[-36AU/5frwp\  , 

where  Tj=o  =  Ts  and  Tj=i  =  Te  represent  the  tunneling  rates  from  the  |  g)  and  |e)  states  of 
the  JJ,  respectively.  The  ratio  T e/T9  ~  250A U/hwp.  Here  A U  =  4/o4>o/3\/27r(l  —  I/Iq )3/2 
is  the  barrier  height  and  ojp  =  21/4  ^2irIo/C<&o(l  —  J/Io)1^4  is  the  plasma  frequency  of 
the  cubic  potential.  The  junction  frequency  uieg  is  related  to  the  plasma  frequency  by 
uieg  ~  ujp{  1  —  5Hlup/36AU).  The  tunneling  rate  of  the  first  excited  state  of  the  junction  is 
then  given  by  Te  ~  500  T9  =  7.3  x  10' s^1  for  AU/frwp  ~  2. 

4.3  Quantum  efficiency 

The  system  is  originally  prepared  in  a  pure  state  pn( 0)  =  \n,g)(n,g\  with  n  photons  in  the 
cavity  and  the  junction  in  the  ground  state  |  <7) .  We  assume  n  photons  are  loaded  into  the 
cavity  in  a  manner  similar  to  that  described  by  Hofheinz  et  al.  [25],  with  loading  rate  faster 
than  the  interaction  rate  between  the  JJ  and  the  cavity.  With  this  assumption,  the  detector 
efficiency  is  unaffected  by  how  photons  are  loaded  into  the  cavity.  We  numerically  solve 
the  above  equation  for  the  time  evolution  of  the  density  matrix  to  compute  the  occupation 
probabilities  of  the  cavity  and  junction  states.  The  probability  that  the  JJ  has  switched  to 
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Figure  4.2:  Switching  probability  P\  (Pq)  vs.  photon  detection  time  for  initial  states  |l,g) 
( 1 0,  <7) ) .  Parameters  used  in  this  plot  are:  junction  T\  =  10ns,  barrier  height  AU/hup  = 
2,  vacuum  Rabi  frequency  Q/2ir  =  200  MHz,  and  cavity  decay  time  l/n=  1  ps.  Here  the 
detuning  between  the  cavity  and  the  junction  A  =  0.  The  solid  blue  curve  is  the  quantum 
efficiency  r/i  =  P\  —  Po  (see  text).  The  maximum  quantum  efficiency  of  the  detector  is  50% 
for  an  optimal  detection  time  of  45  ns  (the  optimal  point  is  marked  by  the  filled  blue  circle). 


the  voltage  state  at  time  t  is  given  by 


Pn(t)  =  1  -  Tr  [pn(t)]. 


(4.7) 


4.4  Result 

We  consider  the  following  set  of  parameters  for  the  JJ-cavity  system,  typical  of  those  realized 
in  experiments  m-  JJ  frequency  ujeg/2ir  =  4.8  GHz,  junction  decay  rate  7  =  108s  x, 
cavity  decay  rate  k  =  106s-1,  and  vacuum  Rabi  frequency  Q/2n  =  200  MHz.  In  Fig. 
we  plot  switching  probabilities  P\(t),  Po(t)  of  the  junction  for  initial  states  |l,g)  (solid 
red)  and  1 0,  <7)  (solid  black),  respectively.  In  this  simulation,  the  parameters  we  consider 
are  AU/fiWp  =  2  and  we  set  the  detuning  between  the  cavity  and  the  JJ  to  zero,  i.e. , 


4.2 
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A  =  uir  —  ueg  =  0.  The  switching  probability  P\  in  Fig.  4.2  features  steps  whose  periodic 
occurrence  is  a  manifestation  of  Rabi  oscillations  of  the  JJ  with  period  2n/Cl  =  5  ns.  This 
result  is  consistent  with  the  picture  that  switching  of  the  JJ  halts  momentarily  when  the 
junction  returns  to  the  ground  state  in  course  of  the  Rabi  oscillations. 


Next,  we  discuss  the  presence  of  the  wide  plateau  of  Pi  in  Fig.  4.2  The  occurrence 
of  this  plateau  can  be  understood  from  the  fact  that  switching  of  the  junction  is  briefly 
frozen  when  the  junction  relaxes  to  the  ground  state  due  to  dissipation.  The  JJ  then 
switches  to  the  voltage  state  after  time  ~  1/T g,  the  characteristic  time  scale  for  switching 
of  the  junction  in  the  case  of  zero  photons.  The  height  of  this  plateau  can  be  estimated  as 
re/  (re  +  7)  ~  0.5,  which  agrees  with  the  numerical  result  in  Fig. 
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In  order  to  determine  the  quantum  efficiency  of  the  detector,  we  must  properly  treat  dark 
counts  due  to  quantum  tunneling  from  the  |0,g)  state  in  the  absence  of  photon  absorption. 
The  quantum  efficiency  rji  of  the  detector  is  defined  as  the  difference  between  the  switching 
probabilities  for  an  initial  state  with  one  photon  Pi(t),  and  for  an  initial  state  with  no 


photons  Po(t):  i]  1  =  Pi(t)  —  Po(t).  The  quantum  efficiency  is  shown  in  Fig.  4.2  by  the 
solid  blue  curve.  For  our  choice  of  parameters  Te  ~  7,  the  detector  has  maximum  efficiency 
of  about  50%  for  the  optimal  detection  time  td  around  45  ns.  Next,  we  demonstrate  that 
the  bandwidth  of  the  Josephson  microwave  photon  detector  is  broadened  due  to  the  finite 
lifetime  of  the  junction  excited  state.  Here,  we  vary  the  frequency  ur  of  the  cavity  and 
compute  the  quantum  efficiency  of  the  detector  for  the  optimal  detection  time  td  obtained 
at  zero  detuning  A  =  0.  The  detector  bandwidth  is  then  given  by  the  detuning  at  which 
the  quantum  efficiency  of  the  detector  is  reduced  to  half  the  efficiency  obtained  at  zero 
detuning.  For  a  dissipation-free  junction,  the  detector  bandwidth  is  approximately  given 
by  the  vacuum  Rabi  frequency.  However,  in  the  presence  of  dissipation  and  tunneling, 
the  first  excited  state  of  the  junction  is  broadened  by  ~  7  +  re.  This  broadening  of  the 
energy  level  roughly  accounts  for  the  increased  bandwidth  of  the  detector.  We  find  that 
bandwidths  are  factors  of  1.6,  2  and  2.3  larger  than  the  vacuum  Rabi  frequency  for  bias 


points  A U/hup  =  2,1.9  and  1.8,  respectively,  as  shown  in  Fig.  4.3  As  we  lower  the 
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A/a 


Figure  4.3:  Quantum  efficiency  vs.  detuning  A/12  for  optimal  detection  time  obtained 
at  A  =  0  for  various  bias  points:  AU/hujp  =  2  (solid  blue),  1.9  (dashed-dot  black),  and  1.8 
(dashed  red).  The  bandwidths  of  the  detector  are  1.611  (solid  blue),  211  (solid  black)  and 
2. 311  (solid  red),  respectively.  The  remaining  parameters  are  as  in  Fig.  4.2 


ratio  AU/hujp ,  the  tunneling  rate  Te  of  the  first  excited  state  of  the  junction  increases. 
This  in  turn  causes  further  broadening  of  the  junction  excited  state  thereby  increasing  the 
bandwidth  of  the  detector.  Next,  we  analyze  the  effect  of  dissipation  and  bias  point  on  the 
efficiency  of  the  detector.  In  Fig.  |4.4[a),  we  plot  the  quantum  efficiency  of  the  detector  for 
different  values  of  the  junction  relaxation  time  T\  from  10  ns  to  500  ns,  keeping  all  other 


parameters  the  same  as  in  Fig.  4.2  We  find  that  the  quantum  efficiency  increases  with 
increasing  junction  relaxation  time  T\  and  reaches  94%  for  T\  =  500  ns  and  for  a  detection 
time  around  95  ns. 

The  change  in  bias  current  I  of  the  junction  modifies  the  ratio  of  barrier  height  A U 
to  the  junction  plasma  frequency  c op.  Taking  different  values  of  this  ratio,  we  compute 
the  efficiency  of  the  detector  at  fixed  relaxation  time  Tj  =  10  ns;  the  results  are  shown 


in  Fig.  4.4  (b).  Upon  decreasing  the  ratio  AU/huip,  the  potential  well  becomes  shallower, 
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Figure  4.4:  (a)  Quantum  efficiency  t] \  vs.  photon  detection  time  for  AU/hup  =  2  and  for 
various  decay  times  of  the  junction:  T\  =  10ns  (solid  blue),  T\  =  20ns  (dashed  green),  Tj 
=  50ns  (dash-dotted  red),  and  T\  =  500ns  (dotted  black).  For  T\  =  500  ns,  the  maximum 
quantum  efficiency  is  94%.  (b)  Quantum  efficiency  rji  vs.  photon  detection  time  for  junction 
Tj  =  10ns  for  various  bias  points  of  the  junction:  AU/hu>p  =  2  (dotted  blue),  1.9  (dash- 
dotted  green),  1.8  (dashed  red)  and  1.7  (solid  black).  The  maximum  quantum  efficiency  is 
84%  for  AU/hiVp  =  1.7.  The  remaining  parameters  are  as  in  Fig. 


4.2 


leading  to  enhanced  tunneling  out  of  the  first  excited  state  and  increased  efficiency  of 
the  detector.  Our  simulation  results  indicate  that  a  significant  improvement  in  detector 
efficiency  is  achieved  when  the  tunneling  rate  exceeds  the  dissipation  rate  of  the  junction. 
We  find  that  for  the  bias  point  AU/hojp  =  1.7,  the  efficiency  of  the  detector  is  about  84% 
for  a  detection  time  around  9  ns.  Finally,  we  analyze  the  efficiency  of  the  detector  to  detect 
single  photons  when  the  cavity  is  loaded  with  n  >  1  photons.  Generalizing  the  previous 
case  of  a  single  photon  in  the  cavity,  the  efficiency  to  detect  a  single  microwave  photon 


in  a  cavity  loaded  with  n  photons  is  given  by  rjn  =  Pn[t)  —  Po(t).  In  Fig.  4.5  we  plot 


the  efficiency  at  fixed  bias  point  A  U/hujp  =  2  and  T\  =  10  ns  for  different  numbers  of 
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Figure  4.5:  Efficiency  r]n  to  detect  a  photon  vs.  detection  time  for  n  =  1  (solid  blue),  2 
(dashed  blue),  and  3  (dash-dotted  blue)  photons  in  the  cavity.  The  rest  of  the  parameters 
are  as  in  Fig.  |4.2|  For  three  photons  in  the  cavity,  the  efficiency  to  detect  a  photon  is  85%. 


photons  in  the  cavity:  n  =  1  (solid  blue),  2  (dashed  blue),  and  3  (dash-dotted  blue).  We 
find  that  detection  efficiency  increases  with  the  increasing  number  of  photons  in  the  cavity 
and  reaches  85%  for  three  photons  in  the  resonator.  This  result  is  consistent  with  previous 
studies  [93l  l9Tj  that  reported  a  higher  detection  efficiency,  for  the  same  parameters  as  above, 
when  a  continuous  flux  of  photons  was  incident  on  the  detector.  For  the  case  of  a  single 
photon  in  the  cavity,  the  detector  returns  to  the  ground  state  after  the  photon  is  absorbed 
by  the  environment,  and  no  further  excitation  of  the  junction  is  possible.  However,  when 
multiple  photons  are  present  in  the  cavity,  other  photons  are  available  to  induce  excitation 
if  the  junction  relaxes,  thereby  increasing  the  probability  of  photon  detection.  We  note  that 
for  multiple  photons  in  the  cavity,  the  measured  efficiency  r]  can  also  be  used  to  estimate 
the  average  number  of  photons  in  the  cavity. 

Recently,  Peropadre  et  al.  [94]  proposed  a  phenomenological  model  for  the  JJ-based 
microwave  photon  detector  that  fails  to  address  specific  microscopic  details  of  experimen¬ 
tally  realized  detectors  m-  Specifically:  (1)  Peropadre  et  al.  treat  tunneling  from  the 
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excited  state  of  the  junction  by  a  non-hermitian  term  in  the  junction  Hamiltonian;  this  is 
not  consistent  with  the  standard  form  of  quantum  tunneling.  (2)  These  authors  do  not 
consider  tunneling  from  the  low-energy  state  of  the  junction,  which  is  responsible  for  dark 
counts  of  the  detector.  (3)  Finally,  their  model  does  not  take  into  account  the  relatively 
strong  relaxation  from  the  excited  to  the  ground  state  of  the  junction.  This  relaxation  cor¬ 
responds  to  a  Xj  time  of  order  of  a  few  nanoseconds  in  present  devices,  which  are  strongly 
coupled  to  a  50  P  readout  line,  and  is  responsible  for  a  significant  suppression  of  detector 
efficiency  m-  If  the  relaxation  time  were  above  500  ns,  the  efficiency  would  reach  nearly 
100%,  see  Fig.  4(a). 

4.5  Input-output  theory  and  reflection  coefficient 

For  practical  realization  of  a  photon  detection  scheme  discussed  in  the  previous  section,  it  is 
important  to  address  several  issues  that  one  may  encounter  in  a  real  experiment.  One  of  such 
concerns  that  directly  impacts  the  detection  efficiency  is  the  power  reflected  by  the  detector. 
Since  our  goal  is  to  maximize  the  detection  efficiency,  it  is  crucial  to  minimize  the  power 
reflected  by  the  detector.  In  this  section,  we  investigate  the  condition  for  minimizing  the 
power  reflected  by  the  photon  detector.  We  apply  input-output  theory  m  to  calculate  the 
reflection  coefficient  of  the  detector  connected  to  a  microwave  source  through  a  transmission 
line,  taking  into  account  switching  behavior  of  and/or  dissipative  effect  in  the  junction. 

Consider  a  Josephson  junction  (JJ)  biased  with  a  current  I  close  to  the  critical  current 
Jo-  The  Lagrangian  of  the  junction  is  given  by: 

£[<!>,  <F]  =  where  t/(4>)  = -Xjcos(^)  -  «F/  and  <F0  =  V2e- 

(4.8) 

If  we  expand  the  potential  function  around  the  local  minimum  <3?mjn  =  ^  sin_1(I//o),  we 
obtain 

U($)  =  - ~  ^min)2  +  0($3)- 


(4.9) 
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Here  Cj  and  Lj  are  the  junction’s  capacitance  and  inductance,  respectively.  We  will  assume 
that  the  JJ  is  positioned  at  x  =  0  and  is  connected  to  a  transmission  line  (left)  whose 
characteristic  impedance  Zi  =  where  Li  and  Ci  are  the  inductance  and  capacitance 

per  unit  length  of  the  transmission  line.  Similarly,  we  model  the  Ohmic  loss  in  the  junction 
(tunneling  in  JJ  can  also  be  modeled  this  way)  by  connecting  the  junction  in  parallel  to  a 
semi-infinite  transmission  line  (right)  whose  characteristic  impedance  Zr  =  where  Lr 

and  Cr  are  the  inductance  and  capacitance  per  unit  length  of  the  right  transmission  line. 

Far  away  from  the  junction  (i.  e.,  x  >  0  or  x  <  0),  the  flux  field  <h^r(x,t)  in  the  left  or 
right  transmission  line  satisfies  the  following  wave  equation: 


d  2®i!r(x,t) 


2  d2®i,r(x,t) 

Vl,7 


dt 2  l’r  dx2 

The  solution  to  the  wave  equation  is  given  by  10  lj 
/  frZir  du 


=  0  ,  where  the  phase  velocity  vi)T  =  Ll^Cl  '■  (4.10) 


$2,r(M)  = 


2  ‘htyfu  Y 


ai 


+  a^(w)ei(=Ffcj,r*-«t)  +  /j.c.  (4.11) 


where  \k^r\  =  The  continuity  of  the  flux  field  at  x  =  0  requires 


$2(0-,*)  =  $r(0+,t). 


(4.12) 


The  discontinuity  in  the  derivative  (w.r.t.  x)  of  the  flux  field  at  x  =  0  leads  to  the  Kirchhoff’s 
current  law,  which  is  equivalent  to  the  Euler-Lagrange  equation  of  motion  derived  from  the 
total  Lagrangian  density  of  the  system  (or  Hamilton  equation  of  motion  in  Hamiltonian 
formulation) : 


i  a$i(o -,*)  i  3<Mo+,t)  ,  ^  52^(o ,t)  ,  $,(o,*)  n 

+  W - 7^ - 1 - - -  =  0  . 


Li  dx 


dx 


dt2 


Lj 


(4.13) 


Plugging  in  the  solution  of  the  wave  equation  into  the  boundary  conditions,  Eqs.  4.12  and 


4.13,  we  obtain: 


^Zia,r\u)  -  y/Z^a^iu)  =  -^fz^Yu)  +  V^r4"C w) 

Xa\u\ u)  +  Ya°ru\u)  =  -X*ajn(uj)  +  Ya™(uj) ,  where 


X  = 


and  Y  =  —i 


Zr 


(4.14) 
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If  we  ignore  the  right  input  amplitude  for  a  moment,  we  get  a°ut(uj )  =  R( u>)  ajn(w),  where 
the  reflection  coefficient  R(oj)  is  given  by: 


upZ\Zr  —  io2CjZjZiZr  +  iio(ZjZr  —  ZjZi) 
u>pZiZr  —  ^CjZjZiZr  —  iu(ZjZr  +  ZjZi) 


Here  up  =  \J\/LjCj  and  Z j  =  \J~LJJC~j  are  the  plasma  frequency  and  the  impedance  of 
the  junction,  respectively.  From  the  above  expression  for  reflection  coefficient,  we  see  that 
for  Z[  =  Zr  and  oj  =  cop,  we  have  R  =  0. 

This  calculation  indicates  that  the  power  reflected  by  the  detector  is  minimal  when 
the  impedance  of  the  transmission  line  is  equal  to  the  impedance  due  to  tunneling  or  any 
internal  loss  in  the  detector.  It  also  requires  a  resonance  condition,  namely,  u  =  ojp. 


4.6  Power  absorbed  by  the  detector 

So  far  we  only  considered  a  detection  model  where  the  resonator  was  loaded  with  a  fixed 
number  of  photons,  which  were  then  detected  by  a  detector  that  was  capacitively  coupled  to 
the  resonator.  Alternatively,  in  a  more  realistic  setting,  one  can  have  a  classical  microwave 
source  directly  connected  to  the  detector,  or  a  microwave  driven  resonator  connected  to  the 
detector.  In  this  section,  we  present  a  rigorous  method  to  compute  the  power  absorbed 
by  the  detector  in  these  alternative  set-ups.  This  is  an  important  calculation  since  one 
can  show  that  the  power  absorbed  by  the  detector  is  directly  related  to  the  efficiency  of 
the  photon  detector.  The  schematics  of  circuit  diagrams  of  a  microwave  driven  resonator 
connected  to  a  JJ  and  a  microwave  driven  JJ  are  shown  in  Fig.  |4.6[ 

Formalism 

Here  we  present  a  formalism  to  compute  the  power  absorbed  by  the  detector  in  the  presence 
of  decoherence.  Since  the  photon  detector  we  consider  here  switches  from  a  supercurrent 
state  to  the  voltage  state  upon  absorbing  photons,  we  would  like  to  characterize  the  wait 
time  between  the  reset  of  the  junction  and  switching  to  the  voltage  state.  The  wait  time  is 
often  specified  by  a  probability  distribution  which  is  known  as  the  waiting  time  distribution. 
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Figure  4.6:  Schematic  circuit  diagram  of  JJ-based  microwave  photon  detector  (a)  driven  by 
a  classical  microwave  source,  (b)  coupled  to  a  microwave  driven  resonator. 

This  distribution  allows  us  to  compute  the  average  switching  rate  of  the  detector,  which 
is  proportional  to  the  power  absorbed  by  the  detector.  In  the  following,  we  provide  a 
derivation  for  the  waiting  time  distribution  [lOo] . 

For  the  sake  of  generality,  we  will  not  specify  the  Hamiltonian  explicitly  and  simply 
denote  it  by  Hsys,  which  can  represent  either  of  the  detection  schemes  shown  in  the  circuit 
diagrams  above.  The  time  evolution  of  the  density  matrix  of  this  system  is  characterized 
by  a  Liouvillian  superoperator  C-totai > 

=  £total[p(t)\  =  A)[/0(O]  +Je[p(t)\  +  Jg[p(t)\ ,  where  (4.16) 

£o[p{t)]  =  -r^[Hsys,  p{t)}  +  £dissip[p{t)\  ■  (4.17) 

Here  Cdissip  is  the  Lindblad  superoperator  that  characterizes  the  dissipation  in  system  in 
consideration  and  Ji  is  the  tunneling  jump  operator  of  the  JJ  from  the  initial  metastable 
state  |i  =  g,e)  (or  supercurrent  state)  to  the  final  voltage  state  |r)  (or  resistive  state), 
defined  as 

Ji[p{t)\  =TiSip(t)  s\  ,  where  s*  =  )T)Jn,r)(n,i|.  (4-18) 

r y  e  are  tunneling  rates  from  the  ground  and  exited  states,  respectively  and  n  is  the  number 
of  photons  in  the  resonator,  if  the  resonator  is  present  in  the  model.  Since  the  tunneling 
rate  from  the  ground  state,  which  contributes  to  the  dark  count  of  the  detector,  is  orders 
of  magnitude  smaller  compared  to  the  tunneling  rate  from  the  excited  state,  we  neglect  the 


50 


tunneling  rate  from  the  ground  state  in  the  following  discussion.  With  this  simplification, 
Co  can  be  re-defined  in  terms  of  effective  non-Hermitian  Hamiltonian  by  absorbing  the 
tunneling  jump  operator  Je  into  it  in  the  following  way: 

A) [p(t)\  =  ( Heffp(t )  -  p(t)H\f^  +  Cdissip[p(t)\ ,  where  (4.19) 

^  ^  i  hr 

Heff  =  Hsys-  2%se.  (4.20) 

For  time-independent  Liouvillian,  the  formal  solution  of  the  density  matrix  equation  can 
be  written  as 

00  ft  ft.n  ft2 

m  =  eCo+Je  Pin  =  V  /  dtn  /  dtn- 1 .  /  dti  Je  eCodn-tn-i)  Je... 

n= 0  J°  J° 

. Jee^2^  Jeelotl  Pin.  (4.21) 

One  can  define  a  conditioned  density  matrix  (un- normalized)  as 

pc({h  ■  ■  ■  tn})  =  e^d-tn)  Je  eCo{tn-tn-i)  je  ■ ■  ■  Je  e4(t2~tl}  Je  etotl  pin  .  (4.22) 


The  normalized  conditioned-density  matrix  is  then  given  by 

7  e^°r  hfOj 

Pnorm  =  6  - i  where  p(0)  =  \n,g)(n,g\.  (4.23) 

Here  W(t )  is  the  normalizing  factor  and  is  also  known  as  the  waiting  time  distribution, 
which  is  given  by 


W{t)  =  re  Tr[se  q(t)  sj]  ,  and  ^  =  £0[£(*)]  with  p(°)  =  P(0)-  (4.24) 

With  this  definition  of  waiting  time  distribution,  we  are  now  ready  to  compute  the  average 
switching  rate  of  the  JJ,  which  is  given  by 

1  W(r)dr .  (4.25) 

T 

The  average  power  absorbed  by  the  detector  can  now  be  defined  in  terms  of  the  switching 
rate  in  the  following  way: 


Pahs  h  Cg 


1 


eg\ 


(4.26) 


Here  uieg  is  the  frequency  of  the  photon  absorbed  by  the  detector,  which  is  equal  to  the 
characteristic  frequency  of  the  JJ. 
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Figure  4.7:  Waiting  time  distribution  of  a  microwave  driven  resonator  for  (n)  =  5,  12/27T  = 
50  MHz,  7  =  0. 


Waiting  time  distribution 

Below  we  compute  and  briefly  discuss  the  waiting  time  distribution  for  two  different  detec¬ 


tion  schemes  shown  in  the  circuit  diagrams  above.  The  system  Hamiltonian  for  Fig.  4.6 
(a)  is  given  by 


H-sys  —  hAesIIe  T  -—((7^  + 


(4.27) 


where  Aeg  =  \uid  —  uieg\  is  the  detuning  between  the  drive  frequency  ujd  and  the  JJ’s  char¬ 
acteristic  frequency  ueg.  Ue  is  the  projection  operator  for  the  excited  state  of  the  junction, 
ne  =  |e)(e|,  and  /  is  the  drive  amplitude.  We  choose  to  write  /  =  \J (n)Q,  anticipating  its 
connection  with  an  alternative  detection  scheme  based  on  a  microwave  driven  resonator,  in 
which  case  (n)  is  the  average  number  of  photons  in  the  resonator  and  12  is  vacuum  Rabi 
frequency.  However,  for  a  JJ  driven  by  a  classical  microwave  source  with  drive  amplitude 
/,  the  quantity  /  also  characterizes  the  Rabi  frequency. 


In  Fig.  4.7 ,  we  plot  the  waiting  time  distribution  of  a  microwave  driven  JJ.  Here  we 
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Figure  4.8:  Waiting  time  distribution  of  a  microwave  driven  resonator  for  (n)  =  5,  £l/2n  = 
200  MHz,  7  =  100  MHz. 


take  (n)  =  5,  H/27T  =  50  MHz  and  relaxation  rate  in  the  junction  7  =  0  and  Te  =  7.3  x  10' 


-,-1 


The  first  peak  on  the  waiting  time  distribution  always  occurs  at  half  the  Rabi  period, 


which  is  given  by  Tr  =  27 r//  and  the  subsequent  peaks  are  separated  by  a  full  Rabi  period. 
This  is  due  to  the  fact  that  the  junction  is  most  likely  to  switch  when  it  reaches  the  excited 
state  and  it  only  takes  half  the  Rabi  period  to  reach  the  excited  state  for  the  first  time, 
but  a  full  Rabi  period  in  all  future  times.  The  number  of  peaks  that  appear  on  the  waiting 
time  distribution  are  characteristics  of  the  Rabi  frequency.  In  addition,  for  non-zero  7, 


these  peaks  are  not  well  resolved,  as  shown  in  Fig.  4.8  In  this  figure,  7  =  100  MHz  and 
H/2vr  =  200  MHz. 

Next,  we  consider  the  waiting  time  distribution  for  a  microwave  driven  resonator  con¬ 
nected  to  a  JJ.  Here  the  system  Hamiltonian  is  given  by 

Hsys  =  -^-(a  +  a))  +  hAr(a)  a  +  — )  +  hAegtle  H — —  (a^  +  o  cr+) ,  (4.28) 

where  the  detuning  Ar  =  \ojd  —  wr|,  ojr  is  the  resonance  frequency  of  the  resonator  and  H  is 
vacuum  Rabi  frequency.  The  drive  amplitude  /  is  given  by  /  =  y/ (n)(4Aj!  +  k2),  where  k 
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Figure  4.9:  Waiting  time  distribution  of  a  microwave  driven  resonator  connected  to  a  JJ 
for  two  different  input  powers:  (n)  =  5  (solid  red  line)  and  (n)  =  11  (dashed  black  line). 
Other  parameters  are  VL/2tt  =  50  MHz,  7  =  0  and  k  =  1  MHz. 


is  the  inverse  lifetime  of  photons  in  the  resonator. 


In  Fig.  4.9,  we  plot  the  waiting  time  distribution  of  a  microwave  driven  resonator  con¬ 
nected  to  a  JJ  for  two  different  input  powers  of  the  drive.  We  find  that  the  waiting  time 
distribution  for  this  detection  scheme  is  very  different  from  that  of  a  microwave  driven  JJ. 
For  a  coherent  state  in  the  resonator,  we  notice  a  far  fewer  oscillatory  peaks  on  W(r). 
This  is  due  to  the  beating  between  the  Rabi  oscillations  with  different  frequencies,  which  is 
caused  due  to  different  number  of  photons  in  the  resonator. 

Next,  we  compute  the  average  switching  rate  of  the  detector.  As  shown  in  the  previous 
section,  the  average  switching  rate  is  given  by 

rP 


dr  W (r)  = 


dr —  Tr 

T 


sP  e 


Cqt 


m 


(4.29) 


We  numerically  diagonalize  Cq  and  obtain  the  eigenvalues  and  eigenvectors,  which  can  be 
written  compactly  as 


D  Cq  D-1  =  A 


(4.30) 
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Figure  4.10:  Average  switching  rates  for  two  different  detection  schemes:  a  microwave 
driven  resonator  connected  to  a  JJ  (red  solid  line)  and  a  microwave  driven  JJ  (solid  black 
line).  Here  n/27r  =  50  MHz,  7  =  0  and  k  =  1  MHz. 


where  D  is  the  matrix  of  eigenvectors  and  A  is  the  diagonal  matrix  of  eigenvalues  of  £0  • 
If  we  plug  this  into  the  expression  for  switching  rate,  we  obtain 


-  )  =  Te  Tr 


/  r°°  pAr  . 

D  (  y  —dTj'D-1Q(0)s 


(4.31) 


Since  DD  1  =  I,  Tr  [se  D  D-1  £(0)  si]  =  Tr  [se  £(0)  s£] ,  which  is  precisely  0.  We  can 
therefore  introduce  an  exponent  of  an  arbitrary  complex  number  A  with  a  negative  real 
part  in  the  time  integrand,  which  allows  us  to  perform  the  time  integration  analytically: 
r(eAT  —  eAr) 


f 


dr  =  —  log[ — A]  +  log[ — A] ,  where  Re  [A,  A]  <  0  . 


In  Fig.  4.10,  we  plot  the  average  switching  rates  for  two  different  detection  schemes. 
We  find  that  average  switching  rates  for  these  detection  schemes  coincide  as  we  increase 
the  number  of  photons  in  the  resonator,  as  expected.  The  deviation  in  the  switching  rate 
from  the  expected  value  for  a  large  photon  numbers  is  simply  a  numerical  artifact  due  to 
truncation  of  the  full  Hilbert  space. 
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4.7  Discussion 

In  conclusion,  we  have  presented  a  model  to  determine  the  quantum  efficiency  of  a  microwave 
photon  detector  based  on  a  current-biased  JJ.  We  have  demonstrated  that  the  efficiency 
to  detect  a  single  photon  loaded  in  a  cavity  has  maximal  value  re/(Te  +  7).  We  have  also 
determined  that  the  bandwidth  of  the  detector  is  characterized  by  the  sum  of  the  vacuum 
Rabi  frequency  and  the  broadening  of  the  first  excited  state  of  the  JJ  due  to  tunneling  and 
relaxation  processes.  Our  simulations  indicate  that  for  currently  used  JJ  photon  detectors, 
the  quantum  efficiency  is  about  50%  for  the  bias  point  A U/hup  =  2  and  about  85%  for 
A U/hujp  =  1.7.  We  have  finally  investigated  the  probability  to  detect  a  photon  in  the 
case  of  a  multiphoton  initial  resonator  state  and  have  found  that  the  detection  efficiency 
quickly  approaches  100%  as  the  initial  number  of  photons  increases,  consistent  with  previous 
studies  [Ml  EEJ  of  a  continuous  flux  of  photons  incident  on  the  detector. 

We  also  considered  alternative  photon  detection  schemes  and  studied  several  issues  that 
are  of  importance  for  analyzing  the  efficiency  of  the  detector  based  on  such  schemes.  First, 
we  derived  a  condition  to  minimize  the  power  reflected  by  the  photon  detector.  We  then 
presented  a  rigorous  method  to  accurately  compute  the  power  absorbed  by  the  detector. 
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Chapter  5 

Effect  of  decoherence  on  the  Berry 
curvature 

5.1  Introduction 

Quantum  computation  based  on  geometric  phases  is  another  emerging  paradigm  that  has 
become  increasingly  popular  due  to  its  resilience  against  errors.  Geometric  phases  emerge 
naturally  in  quantum  mechanics  during  the  adiabatic  evolution  of  the  quantum  mechanical 
state  for  instance,  the  Berry  phase  [106].  Since  the  Berry  phase  depends  on  the  evolution 
path  and  not  on  the  detailed  dynamics,  which  suggests  an  inherent  fault  tolerance,  it  can 
be  a  useful  resource  for  quantum  computation  m-  The  past  measurement  technique 
of  the  Berry  phase  in  superconducting  circuits  was  based  on  an  interference  experiment 
wherein  the  drive  parameter  was  changed  adiabatically  [29].  A  recent  approach  towards 
the  detection  of  the  Berry  curvature  (and  hence  the  Berry  phase)  relies  on  the  dynamical 
response  of  the  generalized  force  to  quench  velocity  of  the  generalized  displacement  m- 
This  approach  is  more  feasible  than  previous  approaches  since  it  does  not  require  a  strict 
adiabatic  condition,  which  is  hard  to  achieve  in  generic  experimental  settings.  Although 
quantum  computation  based  on  geometric  phases  is  robust  against  error,  the  entanglement 
of  the  qubit  system  with  its  environment  makes  it  difficult  to  separate  the  geometric  and 
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dynamical  contributions  to  quantum  gates.  Thus,  it  is  necessary  to  analyze  the  behavior 
of  the  Berry  phase  and  curvature  in  the  presence  of  decoherence.  Here  we  present  results 
for  the  Berry  curvature  measurement  in  the  presence  of  decoherence  using  density  matrix 
formalism  within  linear  response  theory.  We  also  study  the  effect  of  the  third  level  on  the 
measurement  of  the  Berry  curvature  in  phase  or  transmon  qubits. 

5.2  Berry  phase  and  curvature 

The  Berry  phase  is  defined  as  the  phase  accumulated  by  the  wave  function  during  the  adi¬ 
abatic  evolution  around  a  closed  path  in  a  parameter  space  defined  by  I (t) .  This  geometric 
phase  7 b  can  also  be  obtained  by  integrating  the  Berry  connection  A,  which  is  analogous 
to  a  vector  potential  in  the  parameter  space,  along  a  closed  path  V: 


l 


7 b  =  f  A  ■  dl ,  where 


(5.1) 


A  =  i(T(/)|V|T(/))  and  V  =  (dIa,dIy,i 9/J  . 


(5.2) 


The  curl  of  the  Berry  connection  then  gives  the  Berry  curvature,  which  can  be  thought  of 
as  a  magnetic  field  in  the  parameter  space: 


Ffj,v  =  d/jAu  -  dyA^  ,  where  d M  =  d ixt Iyjz. 


(5.3) 


From  the  Stoke’s  theorem,  it  then  follows  that  the  Berry  phase  is  given  by  integrating  the 
Berry  curvature  over  an  area  enclosed  by  the  closed  path,  namely 


(5.4) 


One  can  also  show  that  the  Berry  curvature  is  directly  related  to  the  geometric  tensor  x,  also 
known  as  Fubini’s  metric  on  quantum  rays,  which  is  often  used  to  describe  the  manifolds 
of  adiabatically  connected  wave  functions  T  (I)  |108| : 


(5.5) 
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It  is  straightforward  to  show  that  F^v  =  —2\m.Xtiv  In  this  chapter  we  will  only  focus  on 
the  Berry  phase  or  curvature  of  the  ground  state  d'g  of  the  system  Hamiltonian  H(I),  for 
which  the  geometric  tensor  can  be  written  as 


x%  = 


*g(I)\d^H\Vn(I))^n(I)\duH\yg(I) 

(Eg  -  En)2 


(5.6) 


Here  En  is  the  energy  eigenvalue  of  the  nth  eigenstate  yVn  of  the  Hamiltonian  H(I).  We 
use  this  formula  to  compute  the  analytical  Berry  curvature  for  two-level  and  three-level 
systems.  We  note  that  the  drive  field  I(t)  corresponds  to  the  magnetic  field  in  case  of  a 
spin  qubit  and  to  the  microwave  drive  in  other  solid-state  charge  qubits. 


5.3  Berry  curvature  in  a  two-level  system 

The  Hamiltonian  of  a  spin-1/2  particle  in  the  magnetic  held  /  is  given  by  : 

=  (5.7) 

where  I(t)  =  (Ix(t),  Iy(t),  Iz(t))  is  a  dimensionless  held  amplitude  with  Iz(t)  =  1  and  Larmor 
frequency  A.  For  a  non-zero  magnetic  held  along  the  x-direction  and  a  zero  held  along  the 
y-direction,  we  can  write  the  Hamiltonian  as 

H(t)  =  ~^(az  +  axIx(t))  .  (5.8) 

The  x-component  of  the  magnetic  held  drive  is  set  to  0  at  t  =  0,  which  is  then  ramped 
quadratically  with  time  while  the  y-component  of  the  held  is  set  to  zero  at  all  time,  i.e. , 
Ix(t)  =  1/2 v^t2  with  IX(T)  =  0.5  or  any  other  value,  and  Iy(t)  =  0.  Here  vx  is  the  ramp 
velocity  of  the  x-component  of  the  held  at  the  end  of  the  drive  or  sweep  time  T.  This 
quadratic  protocol  is  slightly  different  from  the  linear  quenching  of  the  x-component  of  the 
drive  proposed  in  Ref.  where  the  quench  velocity  is  held  at  a  constant  value.  However, 
for  the  quadratic  protocol,  the  ramp  velocity  need  not  be  held  constant,  a  convenience  that 
makes  this  protocol  more  attractive  experimentally. 
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The  response  of  the  system  to  changing  magnetic  field  drive  in  x-direction  is  obtained 
by  calculating  the  y-component  of  the  magnetization  (my  ={ay}/2)  for  different  ramp  ve¬ 
locities  vx.  Within  the  linear  response  regime,  the  slope  of  the  curve  my(vx )  determines  the 
Berry  curvature  of  the  system  m-  The  qubit  is  first  prepared  in  the  ground  state  of  the 
Hamiltonian  H(t)  at  t  =  0.  Subsequently,  the  x-component  of  the  drive  is  ramped  quadrat- 
ically  with  time.  After  time  t  =  T,  the  y-component  of  the  magnetization  is  computed  from 
the  density  matrix  as  my(T)  =  Tr[cr  p(T)]/2  =  Imp(T). 


Results 


In  order  to  calculate  the  magnetization,  we  first  need  to  compute  the  time  evolution  of 
density  matrix  of  the  qubit.  First,  we  consider  the  evolution  of  the  system  in  the  absence 
of  environment  and  study  the  effect  of  decoherence  later.  The  density  matrix  equation  in 
this  case  is  simply: 


dp(t) 

dt 


1 

i  h 


(5.9) 


We  perform  a  transformation  (rotation)  of  the  two  level  system  that  follows  the  eigenstate 
basis,  which  we  refer  to  as  “adiabatic  frame” .  When  the  field  is  in  x  —  z  plane,  we  can  write 
the  unitary  transformation  as 


Hu (t)  =  U(t )  H(t)W  (t)  —  ih  U{t)  i  where  U(t )  =  exp[i  6(t)  ay/ 2]  ,  (5.10) 

and  8(t)  =  arctan [Ix{t)/Iz)  is  the  angle  formed  by  the  the  x-component  of  the  drive  with 
respect  to  the  z-direction.  In  this  basis,  the  Hamiltonian  Hu  is  given  by 


Hu(t)  =  -—I(t)az+N(t), 

■r, ,  hA  •. ,  „  hA  ( 

Af{t)  =  — —  0{t)  ay  =  — —  ( 


where 

Iz  dlx  (t)  \ 

(pz+m)  ^  ay)  ■ 


(5.11) 


The  form  of  the  equation  of  motion  of  the  density  matrix  in  the  adiabatic  frame  remains 
unchanged  and  is  given  by 


dpu{t) 


1 

ih 


HU(t),pu(t) 


dt 


(5.12) 


60 


Due  to  non-linear  time-dependent  nature  of  the  rotation  operator  U ,  the  Hamiltonian  picks 
up  an  additional  term  N(t)  in  the  adiabatic  frame.  This  term  is  negligibly  small  only  in 
the  adiabatic  limit,  that  is,  when  6  <C  /(t).  From  this  new  expression  for  the  Hamiltonian, 
it  is  easy  to  see  that  the  y-component  of  the  magnetization  is  given  by 


(ay{t))  =  TV 


h  dlx(t ) 
2  P(t)  dt 


(5.13) 


and  the  Berry  curvature  Fxy  is  given  by  the  quantity  inside  the  bracket  evaluated  at  the 
end  of  the  drive  time  T,  i.e.,  Fxy  =  J2/2/3(T). 

We  now  proceed  to  calculate  the  Berry  phase  by  integrating  the  Berry  curvature  evalu¬ 
ated  at  each  point  on  a  circular  disk  of  radius  Ip(t)  located  Iz  above  the  x  —  y  plane  in  the 
parameter  space  defined  by  I(t): 

r  r 2?r  rFd)  (  j  \ 

1B  =  jFxyixiy  =  ^  (5.14) 


where  Ip(t)  is  the  radial  component  of  the  field,  that  is,  Ip(t)  =  y  J|(f)  +  Iy(t).  This  expres¬ 
sion  for  the  Berry  phase  is  identical  to  a  well  known  expression  obtained  from  the  adiabatic 
measurement  protocol  [23  (1091  illOj.  that  is,  7 b  =  vr(l  —  cos(a)),  where  a  is  the  angle 
between  I(t)  and  the  z-axis  in  the  parameter  space  and  is  given  by  a  =  tan _1(/p(t)/I2). 
Hence,  7 b  can  be  viewed  as  the  solid  angle  subtended  by  the  circular  disk  of  radius  Ip  at 
the  origin  in  the  parameter  space. 

The  Berry  curvature  measurement  discussed  above  requires  measurement  of  the  magne¬ 
tization  in  the  y-direction.  However,  experimentally,  it  is  much  more  convenient  to  measure 
the  expectation  value  of  (<r2).  Therefore,  we  also  consider  the  Berry  curvature  measure¬ 
ment  by  first  measuring  the  (<t2).  In  order  to  achieve  this,  we  let  the  qubit  or  spin  to  freely 
precess  for  time  r  in  the  presence  of  a  constant  magnetic  field  by  turning  off  the  velocity. 
We  then  measure  the  difference  between  the  maximum  and  minimum  values  of  (<t2(t))  for 
each  ramp  velocity  and  extract  the  Berry  curvature.  We  find  that  the  following  relation 
holds  true: 


&{oz) 


max((<T2(r )))  -  min((d2(r))) 


sin 0 Iz  dlx[t) 
2/3(t)  dt 


(5.15) 
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Figure  5.1:  The  left  panel  is  a  free  precession  of  (az)  for  two  different  ramp  velocities: 
vx  =  0.00314  ns-1  (top)  and  vx  =  0.0974  ns-1  (bottom).  The  right  panel  is  the  difference 
between  the  maximum  and  minimum  of  ( az(r ))  vs.  the  ramp  velocity  vx  of  a  quadratic 
drive  with  detuning  A/27 r  =  100  MHz:  analytical  result  without  decoherence  (red  solid 
line),  and  a  full  numerical  solution  (blue  solid  line)  of  the  density  matrix  equation  without 
decoherence. 


On  the  left  panel  of  Fig.  5.1,  we  plot  the  time  evolution  of  (&z(t))  for  various  ramp  veloc¬ 
ities.  On  the  right  panel,  we  plot  the  difference  between  the  maximum  and  minimum  of 
those  oscillations  and  find  excellent  agreement  between  the  numerical  and  analytical  results. 


Using  Eq.  5.15,  the  Berry  curvature  can  then  be  extracted  from  the  slope  of  this  line. 


Decoherence 

Next  we  consider  the  effect  of  decoherence  on  the  Berry  curvature.  We  can  naturally 
introduce  pure  dephasing  and  relaxation  mechanisms  with  the  Lindblad  superoperator  in 
the  eigenbasis  frame  where  the  system  Hamiltonian  is  diagonal.  The  adiabatic  frame  is 
not  suitable  for  this  purpose  since  the  Hamiltonian  Hu  in  this  frame  is  not  truly  diagonal. 
We  note  that  the  adiabatic  frame  would  be  a  true  eigenbasis  frame  if  9  <C  /(i).  However, 
this  is  not  applicable  in  our  case  where  0/I(t)  <  1.  Therefore,  we  perform  an  additional 
transformation  such  that  the  final  Hamiltonian  is  truly  diagonal  or  approximately  diagonal 
with  error  roughly  ^  or  which  we  assume  to  be  negligibly  small.  The 
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unitary  transformation  is  given  by: 

HD(t )  =  T>(t)  Hu {t)T>\t)  —  ih  V(t )  '  wliere  =  exp[— i/3(f)  dx/2] ,  (5.16) 

and  /3(f)  =  arctan(#(f)//(f))  is  the  angle  of  the  x-component  of  the  drive  with  respect  to  the 
z-direction  in  the  adiabatic  frame.  After  the  transformation,  the  Hamiltonian  HD  becomes 


HD(t)  =  ~hY  \Jp{t)  +  Q2(t)  (7 z  +  Q(t) , 


(5.17) 


where  the  “drag  term”  Q(t)  is  given  by 


Q(t)  =  ~y  and  /3(f) 


Iz  I(t)  -  J, 


dt 2 


dlx(t)\2 
dt  ) 


2i 


(5.18) 


Due  to  non-linear  time-dependent  nature  of  the  rotation  operator  T>(t ),  the  Hamiltonian 
picks  up  an  additional  term  Q(t),  which  is  negligibly  small.  However,  for  the  sake  of  clarity, 
we  choose  the  following  sinusoidal  drive 


Ix(t )  =  I  sin(#(f))  and  Iz(t)  =  I  cos (0(f)) 


for  which  the  rotation  operator  T>  is  truly  time  independent: 


V  =  exp 


arctan(fl/J)  ax/2 


(5.19) 


(5.20) 


Here  9  =  Qt  and  I  is  the  amplitude  of  the  sinusoidal  drive,  which,  unlike  in  the  case 
of  a  quadratic  drive,  is  time-independent.  Consequently,  /3(f)  =  0  and  the  Hamiltonian 
becomes  truly  diagonal.  In  the  eigenbasis  frame,  the  evolution  of  the  density  matrix  of  the 
qubit  coupled  to  its  environment  is  governed  by  the  following  master  equation,  which  is 
characterized  by  a  phenomenological  decay  rate  7  ( T\  =  1/7)  and  a  pure  dephasing  rate  7^ 
(Tfi  =  I/7 0)  for  the  qubit: 


dpDit) 

dt 


1 

ih 


HD(t ) , pu{t)  +  £7[/3w(t)]  +  C^[pn{t)}  , 


Dt 


t.D( 


where 


A [pD{t)]  =  7  {^-9°°+  ~  \{v+°-,pD}\  and 

KlpDm  =  io>  ^n0 pd nj  -  ^{nSn0,pD})  ,with  n0  =  |o)(o| . 


(5.21) 
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Figure  5.2:  Magnetization  ( dy{t ))  in  the  eigenbasis  frame  vs.  drive  time  (t)  for  a  sinusoidal 
drive  with  detuning  A/2ir  =  100  MHz  and  amplitude  I  =  \/2:  analytical  magnetization  at 
the  end  of  the  drive  time  T  in  the  absence  of  decoherence  (red  solid  line)  and  a  full  numerical 
solution  (blue  solid  line)  of  the  density  matrix  equation  in  the  presence  of  decoherence  with 
rates  7  =  1  X  107s-1  and  7^  =  1  x  108  s-1.  Here  three  panels  correspond  to  three  different 
drive  velocities  H  =  0.00314  ns-1  (top),  0.053407 ns-1  (middle),  and  0.8513 ns-1  (bottom). 


We  fist  investigate  the  dynamics  of  the  magnetization  for  several  decoherence  parameters. 
The  y-component  of  the  magnetization  in  the  eigenbasis  frame  is  given  by 


{&y  (0)  =  Tr  pD{t)  V  fry  V 


(5.22) 


In  Fig.  5.2,  we  plot  the  y-component  of  the  magnetization  in  the  eigenbasis  frame  as  a 


function  of  time  for  three  different  velocities  of  the  sinusoidal  drive.  Here  we  choose  the 
amplitude  of  the  drive  to  be  I  =  y/2  and  the  angle  between  the  x-component  of  the  drive 
at  the  end  of  the  drive  and  2-axis  to  be  6(T)  =  7t/4.  For  our  choice  of  pure  dephasing 
and  relaxation  times  of  10  ns  and  100  ns,  respectively,  we  observe  rapid  oscillations  in 
magnetization  at  timescales  shorter  than  the  pure  dephasing  time.  Dephasing  essentially 
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Figure  5.3:  Magnetization  {&y  (T))  in  the  eigenbasis  frame  vs.  the  angular  velocity  17  of  a 
sinusoidal  drive  with  detuning  A/2n  =  100  MHz  and  amplitude  I  =  y/2:  analytical  result  at 
the  end  of  the  drive  time  T  in  the  absence  of  decoherence  (red  solid  line)  and  a  full  numerical 
solution  (blue  solid  line)  of  the  density  matrix  equation  in  the  presence  of  decoherence  with 
rates  7  =  1  X  10' s^1  and  7^  =  1  x  108  s_1.  The  plot  in  the  inset  is  of  the  y-component  of 
the  magnetization  in  the  absence  of  decoherence. 


suppresses  these  oscillations  as  time  progresses  while  the  relaxation  forces  the  spin  to  return 
to  the  ground  state. 

Next,  we  plot  the  y-component  of  the  magnetization  vs.  the  drive  velocity  of  a  sinu¬ 


soidal  drive  in  the  eigenbasis  frame,  as  shown  in  Fig.  5.3  Our  numerical  result  indicates 
that  decoherence  has  negligible  effect  on  the  magnetization  and  hence  the  Berry  curvature. 


The  inset  in  Fig.  5.3  displays  the  y-component  of  the  magnetization  in  the  absence  of  de¬ 
coherence.  Here  we  find  that  the  response  is  not  linear  with  respect  to  the  drive  velocity. 
Thus,  we  conclude  that  decoherence  is  essential  for  the  linear  response  of  the  magnetization 
to  the  drive  velocity,  at  least  for  a  sinusoidal  drive. 


The  effect  of  decoherence  on  the  adiabatic  response  of  a  driven  open  quantum  system  was 
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also  previously  studied  by  Avron  et.  al.  m •  They  considered  the  effect  of  pure  dephasing 
on  the  Berry  curvature  in  the  adiabatic  limit  assuming  that  instantaneous  stationary  states 
move  continuously  with  the  control  parameter.  With  this  assumption,  they  introduced  the 
Lindblad  superoperator  to  take  into  account  pure  dephasing  of  the  quantum  system  in  the 
adiabatic  frame.  However,  it  is  important  to  note  that  the  system  Hamiltonian  is  not  truly 
diagonal  in  this  frame.  Nevertheless,  it  is  possible  to  derive  their  result  if  we  look  for  a 
steady  state  solution  of  the  density  matrix  in  the  adiabatic  frame,  in  the  following  form: 

p(t )  =  po  +  5p(t )  ,  where 
p0  =  ^(l-az)  and 

5p(t)  =  A(t)  gz  +  B(t)  ay  +  C(t)  gz  .  (5.23) 

Here  A,  B  and  C  are  unknown  coefficients  whose  first  time  derivatives  are  zero.  In  order 
to  simplify  notations,  here  we  dropped  the  superscript  U  and  chose  to  use  simply  p  to 
represent  the  density  matrix  in  the  adiabatic  frame.  We  then  plug  the  above  ansatz  into 
the  master  equation  and  compare  coefficients  to  obtain  the  following  set  of  equations: 

\  ^  +  hAI(t)  B  —  C  8  =  0 

B'y 

-HAI (t)  A  -  B^ - —  =  0 

AB-'iiC- 1)  =  0.  (5.24) 


We  solve  these  equations  for  B  and  obtain  the  y-component  of  the  magnetization: 


{ay)  =  Tr [/3(f)  Gy]  = 


A2  I2(t) 


[A2  I2(t )  +  (7/2 +  70)2 


2/3(t)j 


1  dlx(t) 


dt 


(5.25) 


From  this,  we  can  easily  read  off  the  Berry  curvature,  which  is  given  by  the  quantity  inside 
the  curly  bracket  in  the  above  expression.  We  find  that  the  correction  due  to  the  decoherence 
on  the  Berry  curvature  of  an  ideal  system  is  given  by  A2  T2{t)/\A2  I2{t)  +  (7/2  +  7(/>)2] , 
which  is  consistent  with  the  result  derived  in  Ref.  uni  that  used  a  different  approach  and 
only  considered  the  effect  of  pure  dephasing. 

In  Fig.  |5.4[  we  plot  the  y-component  of  the  magnetization  for  different  ramp  velocities 
of  a  quadratic  drive  Ix(t )  =  v2  t2/ 2.  We  find  excellent  agreement  between  numerical  and 
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Figure  5.4:  Magnetization  ( cry{T ))  vs.  the  ramp  velocity  vx  of  a  quadratic  drive  with 
detuning  A/27T  =  100  MHz:  analytical  result  in  the  absence  of  decoherence  (red  solid  line), 
analytical  result  with  the  relaxation  rate  7  =  1  GHz  and  the  pure  dephasing  rate  7^  =  1 
GHz  (black  dashed  line) ,  and  a  full  numerical  solution  (blue  solid  line)  of  the  density  matrix 
equation  in  the  presence  of  decoherence. 


analytical  results  in  the  presence  of  decoherence.  Nonetheless,  we  believe  that  the  effect  of 
decoherence  can  be  accurately  described  in  the  eigenbasis  frame  |lT2j.  where  we  discovered 
that  decoherence  does  not  hinder  the  curvature  measurement,  but  instead  assists  in  the 
linear  response  of  the  magnetization,  in  particular  for  the  sinusoidal  drive. 


5.4  Berry  curvature  in  superconducting  qubits 


In  this  section,  we  consider  the  Berry  curvature  measurement  in  superconducting  qubits.  In 
real  qubits  like  phase  or  transmon,  the  presence  of  the  third  level  cannot  be  ignored  and  its 
effect  on  the  measurement  of  the  Berry  curvature  should  be  carefully  studied.  Nonetheless, 
the  results  presented  in  the  previous  section  is  also  applicable  to  a  three-level  system  with 
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anharmonicity  much  larger  than  detuning.  For  definiteness,  here  we  only  consider  a  phase 
qubit.  The  Hamiltonian  of  a  flux-biased  phase  qubit  Hq  is  written  in  terms  of  operators  Q 
and  5,  the  charge  and  phase  difference  of  the  JJ,  respectively: 


Q 


1 


Tn 


2  C  2  L 


2vr 


/0T0 
2t r 


cos  <5 , 


(5.26) 


where  L  (C)  is  loop  inductance  (junction  capacitance),  T;,  is  the  external  magnetic  flux 
applied  to  the  phase  qubit,  Iq  is  the  critical  current  of  the  JJ,  and  To  =  h/2e  is  the  flux 
quantum.  The  qubit  is  capacitively  coupled  to  a  microwave  current  source,  which  is  used  to 
drive  the  qubit.  This  coupling  introduces  a  time-dependent  part  in  the  qubit  Hamiltonian: 


(5.27) 


Here  I(t)  =  Ix(t)  cosw^t  +  Iy(t)  suiuyt  is  the  microwave  current  with  frequency  aid-  Total 
Hamiltonian  of  the  driven  qubit  is  then  given  by  H(t)  =  Hq  +  P(t).  We  truncate  the  qubit 
Hamiltonian  Eq.  5.26,  along  with  the  time-dependent  term  Eq.  |5.27[  to  three  localized  levels 
and  obtain  the  following  Hamiltonian: 


H(t)  =  t  +  a(t)*Xj&j 

j= 1 


+  Hn 


(5.28) 


where  H,  =  \j){j\  is  the  projector  for  the  jt  l  level,  =  \j)(j  —  lj  is  the  raising  operator, 
a(t')  =  (Ix(t)  —  iIy(t))elUdt /2  is  the  amplitude  of  microwave  drive,  A j  =  To (j\S\j  —  l)/27r h 
is  the  matrix  element  of  the  phase  operator,  Uj  =  (e^+i  —  £o)/fr,  £j  is  an  energy  eigenvalue 
of  time-independent  Hamiltonian  and  Hnr  contains  non-resonant  terms,  which  we  neglect. 


We  transform  the  Hamiltonian  Eq.  5.28  to  a  frame  rotating  with  frequency  Ud  with 


respect  to  the  laboratory  frame  and  obtain: 

HR(t)  =  hYj[ 

3= 1 


A  tt  ,  Iy(t)  y 

j — ^~A3aj- lj  — 2rA3aj-i,3 


(5.29) 


where  detuning  Ai  =  A  =  ujq  —  Ud,  A2  =  5an  +  2  A,  and  the  anharmonicity  parameter 
San  =0J\-  2ujo-  We  define  operators  =  \k)(j\  +  \j)(k\  and  ay-k  =  i\k)(j\  -  i\j){k\. 

We  consider  a  setup  where  the  drive  current  Ix(t)  increases  quadratically  with  time. 
Here  Ix(t)  =  v^t2/ 2  and  we  require  that  IX(T)  =  0.5.  The  y-quadrature  drive  Iy(t )  is  set 


to  zero  at  all  times. 
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Figure  5.5:  The  difference  between  the  maximum  and  minimum  of  (<x~(t))  vs.  the  ramp  ve¬ 
locity  vx  of  a  quadratic  drive  with  detuning  A/27 r  =  100  MHz  and  anharmonicity  5an/ = 
200  MHz  (left  panel)  and  1000  MHz  (right  panel).  The  solid  blue  line  is  the  result  obtained 
for  a  three-level  system  from  a  full  numerical  solution  of  the  density  matrix  equation,  the 
solid  red  line  is  the  analytical  result  for  a  strictly  three-level  system,  and  the  black  dashed 
line  on  the  left  panel  and  black  circles  on  the  right  panel  correspond  to  the  analytical  result 
for  a  strictly  two- level  system. 


Results 


For  a  three-level  system,  it  is  convenient  experimentally  to  perform  a  direct  measurement  of 
(az)  instead  of  the  y-component  of  the  magnetization,  in  part  due  to  difficulty  in  performing 
7r/2-rotation  in  the  presence  of  the  third  level.  For  the  sake  of  simplicity,  we  do  not  consider 
decoherence  and  focus  on  the  effect  of  the  third  level  on  the  Berry  curvature  measurement. 


In  Fig.  5.5,  we  plot  the  result  obtained  from  {az}  measurement  as  a  function  of  ramp  ve¬ 
locities  for  two  different  anharmonicities.  For  comparison,  we  also  plot  the  Berry  curvature 
obtained  analytically  for  strictly  two-  and  three-level  systems.  We  find  that  as  we  increase 
the  anharmonicity  5an,  all  three  plots  coincide,  as  expected.  However,  for  anharmonicity 
close  to  that  of  a  transmon  qubit,  <5an/27r  =  200  MHz,  we  find  a  few  percentage  difference 
between  the  result  obtained  from  a  full  numerical  simulation  that  takes  the  third  level  into 
account  and  that  of  analytical  results  for  a  strictly  two-  or  three-level  system. 
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5.5  Discussion 

We  have  presented  both  theoretical  and  numerical  analysis  of  the  effect  of  decoherence  on  the 
measurement  of  the  Berry  curvature  in  superconducting  qubits.  We  find  that  decoherence 
forces  magnetization  to  respond  linearly  to  the  drive  velocity  in  case  of  linear  drives,  thereby 
assisting  in  the  Berry  curvature  measurement.  In  addition,  we  presented  a  formula  for  the 
correction  due  to  dephasing  and  relaxation  to  the  Berry  curvature  in  a  two-level  system 
in  the  adiabatic  frame.  This  result  is  consistent  with  previous  studies  of  the  effect  of  pure 
dephasing  on  the  curvature.  However,  we  clarified  that  the  adiabatic  frame  is  not  well 
suited  for  studying  the  effect  of  decoherence  on  the  curvature.  Furthermore,  we  studied  the 
effect  of  the  third  level  on  the  Berry  curvature  measurement  in  phase  or  transmon  qubits, 
which  are  attractive  candidates  for  experimental  measurement  of  the  Berry  curvature. 
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Chapter  6 

Qubit  relaxation  from 
evanescent- wave  Johnson  noise: 
Dipole  approximation 


6.1  Introduction 


So  far  in  this  thesis,  we  considered  a  phenomenological  model  of  the  noises  and  studied 
their  effect  on  fidelity  of  control  and  measurement  in  superconducting  qubit  devices.  In 
this  chapter,  we  turn  our  attention  to  microscopic  description  of  one  of  such  noises  that  can 
be  a  cause  of  concern  in  semiconducting  qubit  devices. 

Semiconducting  quantum  dots  are  promising  candidates  for  scalable  quantum  informa¬ 
tion  processing.  Quantum  dots  are  realized  in  a  variety  of  experimental  setups,  including 
a  Si  and  GaAs  two-dimensional  electron  gas  rainsim  semiconductor  nanowires  m 
and  carbon  nanotubes  uni-  Several  experiments  performed  on  laterally  coupled  double 
quantum  dots  (DQDs)  have  demonstrated  precise  and  rapid  control  of  the  coupling  be¬ 
tween  electronic  charge  states  and  coherent  manipulation  of  trapped  electrons,  leading  to 
realization  of  a  DQD  as  a  qubit. 


A  schematic  of  the  cross-section  of  such  an  architecture  is  shown  in  Fig.  6.1  These  de 
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vices  are  fabricated  from  heterostructures  of  GaAs  and  AlGaAs  grown  by  molecular  beam 
epitaxy.  Free  electrons  are  introduced  by  doping  the  AlGaAs  layer  with  Si,  which  accumu¬ 
late  at  the  GaAs/ AlGaAs  interface,  deep  down  in  the  minimum  of  the  vertical  potential, 
that  provides  strong  confinement  of  the  electrons  along  the  growth  direction.  However, 
these  electrons  are  free  to  move  along  the  interface,  where  they  form  a  two  dimensional 
electron  gas  (2DEG)  and  can  have  a  high  mobility  due  to  a  relatively  low  electron  density. 
The  quantum  dots  are  then  defined  by  etching  through  the  2DEG  or  by  depositing  metallic 
gates  on  top  of  the  heterostructure.  Applying  a  negative  voltage  on  these  gates  will  deplete 
the  2DEG  below  them,  leading  to  electron  confinement  in  a  desired  location. 

In  almost  all  of  these  implementations,  confinement  and  manipulation  of  an  electron  in 
a  quantum  dot  is  achieved  by  applying  an  electrostatic  potential  through  metallic  gates. 
While  the  metallic  gates  are  crucial  for  qubit  control,  they  can  also  act  as  a  source  of 
decoherence  during  qubit  operations.  This  creates  a  control  -  isolation  dilemma:  connec¬ 
tions  from  the  outside  world  are  what  make  the  devices  useful,  but  they  are  also  sources  of 
decoherence. 

Several  other  decoherence  mechanisms,  such  as  hyperfine  coupling  of  the  trapped  elec¬ 
tron  spin  to  host  lattice  nuclear  spins  in  spin-based  qubits  m  and  electron  coupling  to 
phonon  modes  [119,  1201150]  in  charge-based  qubits,  have  been  previously  studied  in  an  effort 
to  identify  the  major  source  of  decoherence  in  semiconductor  qubits.  A  more  recent  study 
investigated  decoherence  due  to  voltage  fluctuations  in  the  metallic  gates  using  the  lumped 
circuit  model  of  a  DQD  charge  qubit  m-  In  almost  all  of  these  studies,  dm  nu  m  eh 
the  estimated  energy  relaxation  rate  is  at  least  an  order  of  magnitude  smaller  than  the  rate 
observed  experimentally  1  11.  113] .  suggesting  that  a  different  decoherence  mechanism  is 
dominant  in  current  experimental  setups  for  charge  qubits. 

The  relaxation  of  a  charge  or  spin  qubit  can  also  be  induced  by  the  thermal  and  quantum 
fluctuations  of  electromagnetic  fields.  The  fluctuations  of  the  electromagnetic  fields  are 
greatly  enhanced  in  the  vicinity  of  conductors  because  of  the  evanescent  waves  [121,  11221 


[T2Hl[TMIT251IT2fi].  An  illustration  of  the  concept  is  shown  in  Fig.|6.2|  This  evanescent- 


wave 
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Figure  6.1:  Schematic  of  the  cross-section  of  a  quantum  dot  architecture  with  surface 
metallic  gates  (green)  and  back  gates  (blue).  Figure  courtesy  of  Prof.  Mark  Eriksson’s 
group. 


Johnson  noise  (EWJN)  has  been  shown,  both  theoretically  (127)  and  experimentally  [128], 
to  be  an  important  source  of  decoherence  for  atomic  qubits  near  the  metallic  walls  of  a 
trap.  In  this  work  we  investigate  the  effects  of  metallic  device  elements  in  solid-state  qubit 
architectures.  Similar  investigations  have  been  carried  out  previously  using  lumped  -circuit 
calculations  of  Johnson  noise  |59l  129[  51  j .  Here  we  do  the  noise  calculations  taking  into 
account  the  detailed  spatial  dependence  of  the  fields  and  the  important  effects  of  non-local 
corrections  to  the  electromagnetic  response  functions  (1301  131  j. We  focus  particularly  on 
situations  likely  to  be  of  interest  to  experimenters,  such  as  an  electron  in  a  quantum  dot 
close  to  a  gate  electrode. 

We  pause  here  to  mention  that  the  EWJN  originates  from  properties  of  the  metal  near 
its  surface.  The  expressions  for  the  electromagnetic  field  fluctuations  presented  in  this  paper 
are  for  a  conducting  half  space.  However,  our  treatment  can  be  extended  to  calculate  the 
strength  of  EWJN  in  the  vicinity  of  a  conducting  slab  of  finite  thickness.  In  this  case,  the 
electromagnetic  fluctuations  are  independent  of  the  thickness  of  the  slab  as  long  as  this 
thickness  is  significantly  larger  than  the  skin  depth  of  the  metal.  As  such,  we  anticipate 
that  EWJN  may  be  alternatively  interpreted  as  arising  from  overdamped  surface  plasmon 
excitations  which  exist  within  a  skin  depth  of  the  surface  of  the  metal  |126j . 
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Figure  6.2:  An  illustration  of  evanescent-waves  leaking  outside  of  a  conducting  surface. 
This  figure  is  adapted  from  Ref.  m- 

6.2  Formalism 

The  interaction  Hamiltonian  between  the  fluctuating  electric  and  magnetic  fields  and  charge 
or  spin  qubits  within  the  dipole  approximation,  which  we  will  assume  throughout  this 
chapter,  is  given  by  : 

Hm  =  ~d  ■  E(r,  t ) ,  (charge  qubit)  (6.1) 

Hint  =  -A*  •  B(r,  t ) ,  (spin  qubit) ,  (6.2) 

where  d  and  ft  are  electric  and  magnetic  dipole  moments,  respectively.  To  describe  the 
decoherence  of  qubits  resulting  from  the  evanescent  electromagnetic  fields  surrounding  a 
conducting  gate,  it  is  necessary  to  compute  spectral  densities  StJ  of  the  electric  and  magnetic 


74 


field  fluctuations,  which  are  obtained  from  the  cross-correlation  tensor  for  the  fields: 


1  poo 

SFj{r,r',uj)  =  ^  JQ  e^T  ({ (^r)  >  Ei  (^’0)})  dr  (6-3) 

i  r°° 

Sfj  (r,  r ',  u)  =  -  J  elUT  ({ Bt  (r,  r)  ,  Bj  (r  0)})  dr  (6.4) 

Fermi’s  golden  rule  imply  that  the  relaxation  rate  1  /T\  c  s  of  a  (charge  or  spin,  respectively) 
qubit  transition  of  a  particular  frequency  will  be  proportional  to  the  spectral  densities 
at  that  frequency,  evaluated  at  coincident  spatial  arguments  r  =  r'  (within  the  dipole 
approximation) : 


1  d2  qE(-  -  \ 

—  =  wSu(r,r,u) 

1  qB( ->  ->  \ 

—  =  ^(r,r)W) 


(6.5) 

(6.6) 


These  spectral  densities  are  connected  to  the  electric  and  magnetic  susceptibilities  or  re¬ 
sponse  functions 


•  ,oo 

Xij(r,r',u )  =  -J  9(T)elUT  ([Ei(r,T)  ,  Ej  (r',0)])  dr  (6.7) 

•  pOO 

X?j(r,r',uj)  =  -J  6(t)Eut  {[Bi{r,r)  ,  Bj  (r',0)])  dr  (6.8) 

by  fluctuation  dissipation  theorem  m- 

=  h  coth  •  (6-9) 

As  a  result,  the  relaxation  time  TijC  of  a  charge  qubit  with  dipole  moment  d  pointing  in  the 
zth  direction  at  position  r  and  level  separation  u  will  be  given  by 

ijr-  =  j~lmxFi(r,r,uj)  (1  +  N(u,T))  ,  (6.10) 

and  TijS  of  a  spin  qubit  with  magnetic  dipole  moment  fl  in  the  ith  direction  at  position  r 
and  level  separation  lo  will  be  given  by 

J-  =  ^rI™xFi(.r,r,uj)  (i  +  N(uj,T))  ,  (6.11) 

J-  l.s  h 


where  N(lo,T )  =  l/[exp (hui/ksT)  —  1]  is  the  Planck’s  function.  We  deal  first  with  fields 
that  have  been  averaged  over  distances  of  the  order  of  aatomic  and  lmean,  where  aatomic 
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is  an  interatomic  distance  and  lmean  is  the  mean  free  path  in  the  metal.  As  a  result  the 
dielectric  function  e(r,  w)  may  be  treated  as  a  local  function  of  space.  This  approximation 
breaks  down  in  the  near  vicinity  of  the  conducting  surface,  and  the  influence  of  a  nonlocal 
dielectric  response  on  the  electric  and  magnetic  field  fluctuations  is  addressed  below.  The 
conducting  electrodes  in  semiconductor  qubit  architectures  are  commonly  constructed  of 
copper.  Therefore,  all  numerical  results  presented  here  use  the  electronic  properties  of 
copper  near  OK.  In  this  case,  a  local  dielectric  function  may  be  used  for  distances  much 
larger  than  the  Fermi  wavelength  \p  =  yj 4n2h2 /mEp  ~  0.4  nm.  We  work  in  a  gauge 
where  the  scalar  potential  (f>  =  0,  so  that  for  harmonic  fields  we  have  Eu  =  iujAu. 

The  electric  and  magnetic  susceptibilities  can  be  shown  [1331  126]  to  be  directly  related 
to  the  equilibrium  retarded  photon  Green’s  function  by  the  relations 

,2 


or 


Xij{r,r  ,u)  = 
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e0c 
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,2  Dij 


XijiT 5  ^  i  k-0  —  £  (p-  ^ikm^jnpdkdnDmp  (j* >  ^ 
where  i,j  are  Cartesian  indices  that  run  over  x,y,z,  and  Dij  satisfies 


(6.12a) 

(6.12b) 


-Sa  Vz  + 


2  u2e(r,oj) 


+  didj 


Dik  (r,r')  =<53  (■ r-r’)6 _ 


'jk 


(6.13) 


The  geometry  of  a  particular  problem  is  expressed  through  the  function  e  (r,  w).  Eqs.  6.10 


and  |6.11|  assume  the  charge  or  spin  qubit  can  be  adequately  approximated  as  a  point  dipole. 
The  effect  of  a  qubit  with  an  extended  spatial  distribution  will  be  considered  in  future  work. 
The  task  of  computing  in  a  particular  geometry  is,  in  general,  a  complicated  problem 
in  electrodynamics.  In  this  chapter,  we  shall  limit  ourselves  to  the  situation  where  the 
separation  of  the  qubit  from  the  metal  surface  is  much  less  than  any  radius  of  curvature  of 
the  surface  so  that  the  surface  can  be  thought  of  as  flat. 

Let  z  be  the  distance  from  the  surface.  The  result  for  the  spectral  density  of  the  electric 
field  for  local  electrodynamics  has  been  obtained  by  Henkel  et.  al  m- 
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(6.14a) 
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1  f00  n3 

xfz(z,z,u)  =  - - Re/  —dpe2iqzrp(p),  (6.14b) 

47reo  Jo  Q 

where  q  =  u2/c 2  —  p2  for  p2  <  oj2 /c2  and  q  =  i-\/p2  —  u2 /c2  for  p2  >  oj2/c2  is  the  z- 
component  of  the  wavevector,  and  p  is  the  transverse  component.  Our  notation  follows 
that  of  Ford  and  Weber  530]. 


and 


rs  (p)  = 


rP  (p)  = 


q  —  ^/< jj2e/c 2  —  p2 
q  +  ^oj2e/c2  —  p2 

eq  —  \J oj2e/c 2  —  p2 


(6.15a) 


(6.15b) 


eq  +  yju2e/c2  —  p2 

are  the  Fresnel  reflection  coefficients.  The  corresponding  expressions  for  the  spectral  den¬ 
sities  of  the  magnetic  field  are  identical  to  Eqs.  6.14  if  we  multiply  by  1/c2  and  make  the 
replacement  rs  ■H-  rp. 

We  are  interested  in  separations  that  are  sufficiently  small  so  that  retardation,  and  hence, 
radiation  of  the  electromagnetic  field  may  be  neglected.  This  is  known  as  the  quasistatic 
approximation,  and  it  is  formally  employed  by  taking  the  limit  c  — oo.  This  results  in  the 
greatly  simplified  expressions 


Xzziz-.Z^)  =  2XxX(z,z,uj)  =  167r.1eoZ3Im7qr \  >  (6.16a) 

UJ2 

Xz. z(z,z,u)  =  2 Xxx{z,z,u)  =  — - Im(e  —  1) .  (6.16b) 

lb7reom2: 

Employing  this  approximation  has  eliminated  the  functional  difference  between  the 
transverse  and  longitudinal  components  of  the  field  fluctuations.  The  quasistatic  expres¬ 
sions  differ  from  the  exact  values  by  less  than  1%  when  z  <  <5/10,  where  <5  is  the  skin  depth 
of  the  metal.  For  copper  near  absolute  zero,  6  ~  3  ^m.  Eqs.  |6.16|  diverge  as  z  — >  0,  but 
this  divergence  is  not  physical:  it  is  an  artifact  of  treating  the  dielectric  function  as  local 
at  distances  comparable  to  the  interatomic  spacing.  In  any  differential  equation  satisfied 
by  the  spatial  Fourier  components  of  Dij,  a  local  dielectric  function  will  be  independent 
of  the  wavevector  while  a  nonlocal  one  will  have  a  nontrivial  wavevector  dependence.  It  is 
conventional  to  represent  the  spatial  Fourier  components  of  the  nonlocal  dielectric  function 
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as  a  tensor  quantity  in  the  form 

eij(k,uj)  =  ei(k,uj)^-  +  et(k,u;)^6ij-^-S\  ,  (6.17) 

where  we  have  separated  the  function  into  its  longitudinal  ei  and  transverse  et  components. 


Nonlocality  in  the  dielectric  function  changes  the  reflection  coefficients.  In  the  qua¬ 


sistatic  approximation  rp  is  [130] 


i  2p  f°°  j  1 

7T  Jo  k^{k,u) 

rp  , ,  2p  r,  1 

vr  J0  K  k2ei(k,  uj) 

and  rs  becomes 

r  =  w/o°°  dit/{et(k,u)  -  c2  k2/u2}  -  u/q 
1&  fo°°  dK/{et(k,uj)  -  c2  k2/u2}  -  u>/q 

Here  k2  =  p2  +  n2, 

3uj2  {u +  iv)fl({u +  iv)/kvF) 
1  ’  k2Vp  uj  +  +  iv)/kvp) 

u2 

et{k,u)  =  1 - - — -—ft((uj  +  iu)/kvF) 

uj(uj  +  IV) 


(6.18a) 


(6.18b) 


(6.19a) 

(6.19b) 


fi(x)  =  1  -  |  ln(x  +  l)/(x  -  1)  (6.20a) 

ft(x)  =  ^x2  -  ^ x(x 2  -  1)  In  (a;  +  1)/ (x  -  1),  (6.20b) 


v  is  the  electron  collision  frequency,  ujp  =  (ne2/me o)1/2  is  the  plasma  frequency,  and  vf  is 
the  Fermi  velocity. 


Although  Eqs.  6.14  are  derived  assuming  locality,  it  is  a  convenient  fact  |130|  that  these 
equations  are  also  valid  in  the  nonlocal  regime,  as  long  as  the  nonlocal  form  of  the  Fresnel 


coefficients  Eqs.  6.18  are  used. 


6.3  Results 

In  Fig.  |6.3[  we  present  the  zero-temperature  results  for  the  relaxation  time  T±  from  EWJN 
for  a  qubit  with  an  electric  dipole  moment  of  magnitude  |e|aB,  where  clb  is  the  Bohr  radius. 
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Figure  6.3:  Plot  of  the  T\  time  of  a  charge  qubit  computed  from  Eq.  (1)  using  the  local 
approximation  (dashed  blue)  and  the  full  nonlocal  theory  (solid  red).  We  used  the  values 
Ep  =  7  eV,  oj  =  6i r  x  108  s_1,  v  =  n  x  1013  s~3,  and  ujp  =  1.6  x  1016  s-1,  appropriate  for 
a  copper  surface  and  a  device  operating  in  the  GHz  range.  The  dipole  moment  is  taken  as 
d  =  |e|oB,  where  |e|  is  minus  the  charge  on  the  electron  and  ap  is  the  Bohr  radius.  These 
results  are  for  zero  temperature.  The  dashed  horizontal  line  in  the  left  figure  represents  the 
strength  of  the  electric  field  fluctuations  inside  the  bulk  of  a  uniform  metal. 


Both  the  local  and  nonlocal  results  are  shown.  It  is  seen  that  the  correct  nonlocal  dielectric 
function  eliminates  the  unphysical  divergence  of  1/Ti  at  z  =  0.  For  separations  z  ~  A p,  the 
differences  are  very  significant,  while  for  z  >  lmean  =  vf/v  ~  200Air,  the  local  and  nonlocal 
results  nearly  coincide.  It  is  interesting  to  note  that  for  the  electric  field  fluctuations  there 
is  a  crossover  region  where  the  nonlocal  result  becomes  larger  than  the  local  result  in  the 
range  A f  <  z  <  lmean  (see  Fig.  6.3),  in  alignment  with  the  results  of  Volokitin  et.al  |125| 
who  showed  an  enhancement  of  the  nonlocal  result  above  the  local  result.  We  see  that  at 
T  =  0  and  GHz  operations,  T\  from  spontaneous  emission  is  of  the  order  of  seconds  at 
separations  z  ~  30 A p-  These  results  are  directly  applicable  to  atomic  qubits,  but  the  rate 
1/Ti  is  proportional  to  the  square  of  the  dipole  moment,  so  rates  for  other  charge  qubits 
are  easily  deduced. 
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Figure  6.4:  Plot  of  the  T\  time  vs.  frequency  w  of  a  charge  qubit  computed  from 
at  different  temperatures  (solid  red  at  0  K  and  dashed  red  at  25  rnK)  at  fixed  2 
Ep  =  7  eV,  v  =  7r  x  1013  s-1,  up  =  1.6  x  1016  s_1,  and  d  =  \e\ap,  as  in  Fig.  6.3 


Eq.  (1) 
=  30Af. 


Fig.  6.4  shows  that  T\  falls  off  slowly  at  higher  frequencies.  The  divergence  of  T\  as 
oj  — y  0  is  removed  by  including  a  small  finite  temperature.  Fig.  |6.4|  is  plotted  using  the 
non-local  expression  for  the  z-component  of  the  field  fluctuations. 


Fig.  6.5  gives  the  analogous  results  for  magnetic  EWJN  on  a  spin  qubit  with  a  mag¬ 
netic  dipole  moment  of  1  Bohr  magneton.  Nonlocal  corrections  are  somewhat  stronger  for 
this  case,  and  persist  to  larger  distances.  Interestingly,  the  falloff  with  distance  of  T\  is 
slower  for  magnetic  EWJN  than  for  electric  EWJN.  However,  magnetic  relaxation  times 
are  typically  somewhat  larger  than  electric  relaxation  times.  The  crossover  of  the  local 


and  nonlocal  results  is  not  present  in  the  magnetic  case.  Fig.  6.6  shows  that  the  frequency 
and  temperature  dependence  of  magnetic  EWJN  is  similar  to  the  electric  EWJN  shown  in 


Fig.  6.4  Brief  mention  should  be  made  of  the  dip  that  is  observed  as  oj  — >  0  in  Fig.  6.6 


The  reflection  coefficients  rs  contribute  to  the  magnetic  field  fluctuations.  This  contrasts 
with  the  electric  case,  where  only  rp  contributes  to  leading  order  inw/c.  The  dip  is  a  result 
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Figure  6.5:  Plot  of  T\  time  for  a  spin  qubit  at  zero  temperature  computed  from  Eq.  (2)  in 
the  local  approximation  (dashed  blue)  and  the  full  nonlocal  theory  (solid  red).  Ep  =  7  eV, 
oj  =  6n  x  108  s_1,  v  =  n  x  1013  s  ,  and  ojp  =  1.6  x  1016  s-1,  appropriate  for  a  copper  surface 
and  a  device  operating  in  the  GHz  range.  We  have  taken  g  =  gfiR.  appropriate  for  a  single 
electron,  where  the  gyrometric  factor  of  an  electron  g  ~  2.  The  rate  1/Ti  is  proportional  to 
the  square  of  g,  so  rates  for  other  local  magnetic  qubits  can  be  easily  deduced.  The  dashed 
horizontal  line  in  the  left  figure  represents  the  strength  of  the  magnetic  field  fluctuations 
inside  the  bulk  of  a  uniform  metal. 


of  the  contributions  of  rs  to  the  field  fluctuations.  The  rs  term  in  xfi  is  linear  in  w  as  w  ->  0 
and  negative,  while  the  rp  term  is  cubic  in  oj  as  oj  — >  0  and  positive.  We  should  also  mention 
that  for  extremely  small  oj  the  factor  of  N(oj,T)  cancels  out  the  linear  oj  dependence  and 


the  dip  flattens  out  as  oj  — >  0  (not  observable  in  the  resolution  of  Fig.  6.6). 

The  limit  of  the  nonlocal  quasistatic  field  fluctuations  as  z  — >  0  should  be  of  the  same 
order  of  magnitude  as  the  value  of  these  fluctuations  inside  the  metal.  To  check  this,  we 
calculate  the  electromagnetic  Green’s  function  inside  the  bulk  of  a  uniform  metal  using  a 
nonlocal  dielectric  function.  The  result  is 

iv'l  Ibj  IX/'IIXjj 


Dij(k,oj )  = 


1 


c xj2et/c2  —  k2 


Sij  ~  ^P-  +  P(et  ~  et) 
oj2€i  k~ei 


(6.21) 
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Figure  6.6:  Plot  of  the  T\  time  vs.  frequency  of  spin  qubit  computed  at  different  tem¬ 
peratures  (solid  red  at  0  K  and  dashed  red  at  25  mK)  at  fixed  z  =  30A f-  Ep  =  7  eV, 
v  =  7r  x  1013  s-1,  and  ujp  =  1.6  x  1016  s-1,  appropriate  for  a  copper  surface,  /i  =  g/ap  as  in 


Fig.  6.5 


Dij  [f  —  r',u)  = 


1 


(2vr) 


d3 


(6.22) 


Numerical  evaluation  of  Eq.  6.22  when  r  =  r'  gives  (imaginary  parts)  Dxx  =  Dzz  ~ 
3  x  1019  m  .  An  evaluation  of  the  Green’s  function  outside  the  metal  in  the  nonlocal 
quasistatic  regime  for  z  — >  0  gives  (imaginary  part)  Dxx  ~  1.25  x  1019  m”1  and  Dzz  ~ 
2.47  x  1019  m_1,  slightly  less  than  in  the  bulk  of  a  uniform  metal,  as  expected. 

We  see  that  there  are  three  relevant  distance  regimes.  For  z  <  Xf,  a  quasistatic  approx¬ 


imation  to  Henkel’s  results  6.14  using  the  nonlocal  expressions  for  rs  and  rp  will  accurately 
describe  the  field  fluctuations.  For  intermediate  distances  A p  <  z  <  lmean,  there  is  a  sig¬ 
nificant  enhancement  in  the  electric  field  fluctuations  of  the  nonlocal  expression  compared 
to  the  local  expression.  For  distances  lmean  <  z  <  10  5,  local  forms  Eq.  |6.14|  will  accurately 


describe  the  held  fluctuations,  while  for  z  »  5 ,  retardation  effects  must  also  be  taken  into 


account. 
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Comparison  of  our  results  with  experimental  measurements  of  Xj  times  for  spin  relax¬ 
ation  in  single-electron  quantum  dots  supports  the  notion  that  relaxation  from  EWJN  may 
sometimes  constitute  the  dominant  relaxation  mechanism  in  present  semiconductor  qubit 
architectures.  In  particular,  we  reference  the  low  B  saturation  of  T\  for  a  single  electron  in 
Fig.  |6.6p)  of  [134|  by  Xiao  et.al.  In  their  measurements,  T\  saturates  to  ~  40  ms.  Plug¬ 
ging  their  values  of  z  =  50  nm  and  hu  =  0.4  meV  into  Eq.  6.16b  gives  a  value  of  Xj  ~  15 
ms.  That  our  expression  gives  a  shorter  relaxation  time  than  what  was  measured  can  be 
understood  because  we  assume  the  gates  constitute  a  conducting  half-space,  rather  than 
the  more  sparse  gate  geometry  used  in  [134| .  It  is  also  instructive  to  compare  our  results 
with  the  measurements  of  Elzerman  et  al.  in  m-  They  measured  a  T\  of  0.55  ms  at  a 
magnetic  held  of  10  T.  Using  their  values  of  z  =  90  nm  and  hu  =  0.2  meV,  our  results 
predict  T\  ~  1  s.  EWJN  is  clearly  not  the  dominant  relaxation  mechanism  in  this  situation. 
Through  a  perusal  of  the  literature,  we  have  found  generally  that  EWJN  is  insufficient  to 
describe  the  spin  relaxation  rate  in  experiments  with  a  high  external  magnetic  held. 


6.4  Discussion 

The  density  of  photon  states  in  a  metal  is  very  high  owing  to  the  large  polarizability. 
For  blackbody  radiation,  this  high  density  of  states  does  not  matter,  since  total  internal 
rehection  reduces  the  outgoing  radiation  hux  to  its  universal  Stefan-Boltzmann  value.  In 
contrast,  the  evanescent  waves  are  strongly  enhanced  and  the  resultant  electromagnetic 
noise  just  outside  the  surface  can  be  intense.  This  is  a  concern  for  quantum  devices  operating 
close  to  metallic  objects.  This  paper  has  concentrated  on  the  frequency,  temperature,  and 
distance  dependence  of  the  noise,  and  on  the  effects  of  assuming  a  local  dielectric  function. 
We  conclude  that  the  effect  is  signihcant  for  charge  qubits  with  large  dipole  moments  such 
as  double  quantum  dots.  The  EWJN  relaxation  may  be  the  limiting  decoherence  effect  in 
designs  which  involve  close  proximity  to  bulk  metals.  For  magnetic  qubits  the  effects  are 
smaller,  but  as  can  be  seen  by  comparison  to  m,  they  can  still  be  significant.  We  found 
that  nonlocal  effects  are  very  important  at  short  distances  -  indeed,  local  calculations  can 
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produce  spurious  divergences.  At  distances  large  compared  to  the  Fermi  wavelength  of  the 
metal,  local  approximations  work  well. 

We  have  not  considered  extended  qubits  for  which  the  off-diagonal  function  Da(f,r') 
at  r  /  f'  is  required.  The  extended  qubit  geometries  will  be  considered  in  detail  in  the 
next  chapter. 
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Chapter  7 

Qubit  relaxation  due  to 
evanescent-wave  Johnson  noise: 
Beyond  dipole  approximation 

7.1  Introduction 

In  this  chapter,  we  continue  our  discussion  about  the  decoherence  in  a  quantum  dot  due 
to  electromagnetic  field  fluctuations  near  the  metallic  gates.  In  Fig.  \7.1\  we  show  the 
schematic  cross-section  of  the  type  of  device  we  consider  in  this  chapter  [136] .  The  picture 
of  the  actual  device  is  shown  on  the  right  panel.  The  previous  chapter  [6j  as  well  as  other 
theoretical  estimates  [127]  of  the  effect  of  Johnson  noise  in  atomic  and  quantum  dot  based 
qubits,  use  the  dipole  approximation,  which  is  a  valid  approximation  if  the  distance  from 
the  metallic  gate  to  the  qubit  is  much  larger  than  the  size  of  the  qubit.  However,  it  may  be 
necessary  to  go  beyond  the  dipole  approximation  in  the  case  of  EWJN  in  a  quantum  dot. 

Here  we  present  our  study  of  the  energy  relaxation  of  a  single  electron  charge  qubit 
in  a  DQD  system  and  a  single  electron  spin  qubit  in  a  single  quantum  dot.  We  assume 
that  the  primary  source  of  field  fluctuations  are  the  metallic  top  gates  of  the  quantum  dot 
architecture.  Back  gates  are  typically  a  distance  on  the  order  of  a  micron  from  the  qubits, 
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Figure  7.1:  (a)  Schematic  cross-section  of  the  device,  (b)  Scanning  electron  micrograph  of 
the  actual  device.  Pictures  reprinted  with  permission  from  Ref.  [136] .  Copyright  2011,  AIP 
Publishing  LLC. 


which  is  too  far  to  experience  significant  EWJN  enhancement.  We  consider  the  detailed 
spatial  variation  of  the  electromagnetic  field  fluctuations  and  present  results  beyond  the 
dipole  approximation  which  take  into  account  the  finite  size  of  the  quantum  dot.  We  show 
that  this  extension  of  the  dipole  approximation  removes  the  unphysical  divergence  in  the 
relaxation  rate  at  the  metallic  surface. 

This  chapter  is  organized  as  follows:  In  Section  |7.2|  we  present  our  formalism  for  cal¬ 
culating  the  relaxation  rate  of  a  charge  qubit  beyond  dipole  eapproximation.  Results  are 


presented  for  a  DQD  geometry.  Section  7.3  presents  the  formalism  and  results  for  the  re¬ 


laxation  rate  of  a  spin  qubit  in  a  single  QD.  In  Section  7.4  we  derive  an  enhancement  of 
the  noise  spectrum  that  results  as  the  thickness  of  the  metallic  gate  is  decreased.  Finally, 
Section  [7.51  summarizes  our  results.  Our  results  indicate  that  EWJN  is  the  dominant  cause 
of  energy  relaxation  in  some  spin  qubit  experiments,  particularly  those  performed  in  a  small 
external  magnetic  field,  and  is  comparable  in  effect  with  previously  studied  noise  sources 
in  charge  qubits. 
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7.2  Charge  qubit 

We  consider  a  charge  qubit  realized  in  a  gated  lateral  DQD  in  an  AlGaAs/GaAs  het¬ 
erostructure  where  electron  confinement  along  the  2  direction  is  much  smaller  than  in  the  x 
or  y  directions,  so  that  we  can  safely  decouple  the  dynamics  along  x  and  y  directions  from 
the  z  direction  1 .  The  total  Hamiltonian  of  the  charge  qubit  and  its  interaction  with  the 
electromagnetic  environment  is  given  by 


H  —  Hq  +  Hint 


(7.1) 


where  Hq  is  the  Hamiltonian  of  the  charge  qubit  in  a  DQD,  which  we  model  in  the  basis  of 
the  localized  charge  states  {|L),  | R)}  as  Hq  =  e/2(\L){L\  —  |i?)(i?|)  +  A/2(|L)(i?|  +  \R)(L\). 
e  is  the  bias  energy  between  the  two  dots,  and  A  is  the  tunneling  amplitude.  In  the  energy 
eigenbasis  this  Hamiltonian  reduces  to 

Hq  =  (7.2) 


where  az  is  the  Pauli  matrix,  and  Hu  =  \J e1  +  A2.  For  all  our  calculations  except  those 


in  Fig.  7.3,  we  will  set  e  =  0.  The  interaction  Hamiltonian  Hjnt  may  be  expressed  in  this 
same  basis  as 


Hint  —  /  dv 


ax  Mr(r)  +  az  M^r)  ■  A(f,t) 


(7.3) 


where  A(r,  t )  is  the  vector  potential  of  the  fluctuating  field.  Mr  and  are  associated  with 
energy  relaxation  and  pure  dephasing  in  the  charge  qubit,  respectively  and  are  defined  as 


Mr(r)  = 


2m  - 


ipX (r)  pip+(r)  —  ?/>*  (r)  p ip-  (r) 


ieh 
4  m 


%p*+  (r)  ip+  {f)  V  —  ip*_(r)  V’- (r)  V 


Here  m  is  the  effective  mass  and  p  is  the  momentum  operator  of  the  trapped  electron. 
Because  we  are  operating  within  the  weak  field  limit,  the  term  proportional  to  A2  has  been 
dropped  from  the  interaction  Hamiltonian.  We  choose  the  gauge  where  the  scalar  potential 

1In  the  ID  case  of  a  DQD  realized  in  carbon  nanotube  Ena,  or  semiconducting  nanowire  we 

assume  the  electron  dynamics  along  the  x  direction  are  decoupled  from  that  of  the  y  and  z  directions. 
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cj)  =  0  so  that  E  =  —dtA.  The  expression  for  Hmf  derives  from  an  interaction  in  terms  of 
operator  quantities  of  the  form 

Hmt  =  -  ^  ( A(r,  t)-p  +  p-  A(r,  t)j  .  (7.4) 

This  symmetrized  version  of  the  vector  potential  is  not  strictly  necessary  in  our  case  since  we 
are  only  interested  in  the  near- field  contribution  to  the  Greens  function  for  which  V  •  A  =  0, 
but  we  included  it  to  keep  our  results  more  generally  applicable.  The  relaxation  rate 
IT,c  =  1  /T\,c  can  be  calculated  using  the  following  expression,  which  follows  directly  from 
the  Fermi  golden  rule: 

J  d3ff  J  d3r'  M*l{r)M3r{r')Sij{r,r'  ,uj)  ,  (7.5) 

v 

where  S{j  is  the  spectral  density  of  the  vector  potential: 

1  f°°  ■ 

Sij(f,f',u )  =  -  eluJT  ({Ai(f,T)  ,  Aj  (r',0)})  dr  (7.6) 

The  spectral  density  is  related  to  the  response  function  or  susceptibility  of  the  vector  po¬ 
tential, 

■  poo 

Xij(f,f',u)  =  -  J  0(r)eluJT  ([Ai(r,T)  ,  Aj  (f',0)])  dr  ,  (7.7) 

via  fluctuation-dissipation  theorem  [|132j : 

Sij(r,r*,u)  =  h  coth  (tJ^)  ,u) ,  (7.8) 

where  Xij(r,rf ,lj)  is  the  imaginary  part  of  the  susceptibility  and  is  given  by 

Xij(r,f,u} )  =  Im  Xij{r,7*,u)  =  —  J  eluJT  ([Ai  (r,  r)  ,  Aj  (f1, 0)])  dr.  (7.9) 

Combining  the  Fermi’s  golden  rule  with  fluctuation-dissipation  theorem,  the  zero  temper¬ 
ature  relaxation  rate  can  be  written  as: 

Ti,c  =  t  y~]  J  d3f  J  d3f '  M™(r)Mi(r')Imxij(r,r',u) , 


(7.10) 
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At  finite  temperature,  the  emission  (transition  from  excited  to  ground  state)  and  absorption 
(transition  from  ground  to  excited  state)  rates  are  given  by 


r?iC  =  (1  +  N(ui,T))TitC  , 

r“c  =  n(u,  r)rljC  (7.11) 


The  Planck  function  N(u,T )  =  l/[exp(hu/kBT)  —  1]  gives  the  average  occupation  number 
of  environment  modes  with  frequency  uj  at  temperature  T.  The  response  function  or  sus¬ 
ceptibility  of  the  vector  potential  \ij  is  related  to  the  retarded  photon  Green’s  function  D^ 

by  [Ml 

Xij  (' r,f',u} )  =  Dij  ( r,r',u )  , 


where  i,j  are  Cartesian  indices  that  run  over  x,y,z  and  the  square  brackets  denote  the 
commutator.  D%j  is  obtained  by  solving 


Sij  [  V"  + 


2  u2e  (r,uj) 


+  didj 


Dik  (' r,r',u )  =  53  (f-  r')  Sjk. 


(7.12) 


Here  the  relative  permittivity  e(r,ui)  characterizes  the  geometry  of  a  particular  problem. 
In  this  section,  we  shall  limit  ourselves  to  the  case  where  the  metallic  top  gate  of  the  lateral 
DQD  is  approximated  by  the  half-space,  z  <  0.  Then  we  can  derive  an  analytical  expression 
for  Di:j  [Ml  \m 


Dij(r ,f' ,u)  =“j~2  J  e*k'*ll.Djj(fc,  z,z',uj) 


dk . 


Dxx(k,z,z',u )  =— e 


J^Jq(z+z') 


q2c2 

rs(k ,  ui)  cos2  9 - irrp{k, w)  sbr2  ^ 

UJZ 


where  rs  and  rp  are  Fresnel’s  reflection  coefficients  given  by 


(7.13) 

(7.14) 


rp(k,u) 


eg  -  q i 
eg  +  gi 


rs(k,u) 


g  ~  gi 
g  +  gi  ’ 


Here  k  =  (kx,  ky ),  l||  =  rj|  —  r  m  =  (x  —  x',y  —  y'),  6  is  the  angle  between  k  and  the  x-axis,  and 
q  =  yju i2 1  ft  —  k 2  and  q\  =  \J  ew2/c2  —  k 2  are  the  z-components  of  the  photon  wavevector 
in  the  vacuum  and  the  metal,  respectively.  All  other  components  of  Dij  can  be  derived  from 
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Dxx .  [133]  In  this  work  we  consider  e  ~  ia/uue o,  where  a  =  6  x  10'  S/m  is  the  conductivity 
of  the  copper  gate. 

We  pause  briefly  to  mention  that  typical  models  of  the  interaction  of  a  DQD  with  the 
electromagnetic  held  use  the  dipole  interaction  Hamiltonian 


Hint  =  ~E{r)  ■  d, 


(7.15) 


which  will  result  in  a  relaxation  rate  of 

IYC  = 


d2u 


(7.16) 


87 vhz3a 

in  the  quasistatic  approximation,  where  d  is  the  dipole  moment  of  the  qubit  and  E(r)  is 
the  strength  of  the  fluctuating  electric  held  evaluated  at  the  location  of  the  qubit.  Here 
we  also  assume  that  the  dipole  moment  points  along  the  z-direction.  This  expression 
approximates  that  the  electric  held  is  uniform  over  the  spatial  extent  of  the  qubit,  which 
is  equivalent  to  treating  the  qubit  as  a  point  dipole.  Here  we  also  .  As  such,  the  qubit  is 
able  to  couple  to  arbitrarily  small  wavelengths  of  the  electromagnetic  spectrum,  and  the 
relaxation  rate  is  seen  to  diverge  at  shorter  distances  as  ~  1  /z3  if  the  conductor  is  modeled 
with  a  local  dielectric  function  [1271  [33] .  Using  the  complete  electromagnetic  interaction 
Hamiltonian  |7.4|  accounts  for  fluctuations  of  the  held  over  the  spatial  extent  of  the  qubit. 
If  the  wavelength  of  a  particular  Fourier  component  of  the  held  fluctuations  is  smaller  than 
the  length  of  the  qubit  in  that  direction,  its  influence  on  the  electron  will  average  out  and 
it  will  not  contribute  to  qubit  relaxation.  The  exact  and  dipole  approximation  forms  of 


the  interaction  Hamiltonian,  Eqs.  7.4  and  7.15,  converge  when  the  distance  from  the  gate 
becomes  larger  than  the  spatial  extent  of  the  qubit. 

We  present  calculations  of  the  relaxation  time  for  charge  qubits  that  highlight  the  dif¬ 
ferences  between  these  two  forms  of  the  interaction.  First  we  consider  DQDs  in  a  one 
dimensional  nanowire,  which  are  realizable  in  semiconducting  nanowires  or  carbon 

nanotubes  m-  In  such  a  geometry,  the  wave  functions  of  trapped  electrons  in  quantum 
dots  have  appreciable  spatial  extent  in  only  one  direction.  We  model  the  confining  poten¬ 
tial  of  the  DQD  as  a  symmetric  double  square  well  potential  and  compute  the  lowest  two 
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Figure  7.2:  Energy  relaxation  time  T\  vs.  the  distance  from  the  metallic  gate  to  the  DQDs 
z  for  dot  geometry  of  d  =  30  nm  and  l  =  30  nm  (dashed  and  dash-dotted  blue  lines)  and  60 
nm  (solid  and  dashed  black  lines)  at  0  K,  uj/2it  =  1  GHz,  and  e  =  0.  Solid  and  dash-dotted 
lines  are  T\  times  for  the  exact  form  of  the  interaction  Hamiltonian,  whereas  dashed  lines 
are  for  the  dipole  form  of  the  interaction.  The  inset  in  the  figure  displays  confining  potential 
of  a  typical  DQD  in  a  one  dimensional  nanowire  and  the  corresponding  symmetric  ground 
state  and  antisymmetric  first  excited  state. 


eigenenergies  and  wavefunctions.  We  then  compute  the  relaxation  rate  between  these  two 
lowest  states  which  are  separated  by  a  fixed  transition  frequency  ijo/2tt  =  1  GHz.  A  plot  of 


the  wave  functions  and  the  shape  of  the  potential  is  shown  in  the  inset  of  Fig.  7.2 


We  plot  the  energy  relaxation  time  T\  vs.  the  distance  z  from  the  metallic  gate  to  the 


DQD  in  Fig.  7.2  In  this  plot,  we  choose  the  size  of  the  dot  in  the  x-direction  d  =  30  nm  and 


half  the  separation  between  the  dots  l  =  30  nm  (dash-dotted  and  dashed  blue  lines)  and 
60  nm  (solid  and  dashed  black  lines).  The  curves  that  are  shown  in  solid  and  dash-dotted 
lines  are  relaxation  times  for  the  exact  form  of  the  interaction  Hamiltonian,  whereas  those 
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Figure  7.3:  Ratio  of  the  T\  time  for  finite  bias  e  to  the  T\  time  at  zero  bias  vs.  s/A  for 
three  values  of  distances  2:  from  the  metallic  gate:  1  nm  (solid  blue  line),  10  nrn  (dashed 
red  line)  and  50  nm  (dash-dotted  black  line),  for  a  dot  geometry  of  d  =  30  nm  and  l  =  30 
nm  at  0  K  temperature. 


shown  in  dashed  lines  are  obtained  using  the  dipole  form  of  the  interaction.  The  curves 
show  significant  deviation  of  the  exact  relaxation  rate  from  the  dipole  relaxation  rate  at 
shorter  distances  and  convergence  of  the  two  results  at  longer  distances. 


In  Fig.  7.3,  we  present  the  ratio  of  T\  for  a  charge  DQD  qubit  at  bias  e  to  the  T\ 
obtained  at  £  =  0  vs.  the  ratio  s/A.  An  increase  in  bias  increases  the  level  splitting  and 
decreases  the  dipole  moment  of  the  DQD.  Since  the  relaxation  time  T\  ~  l/uid !2,  where 
d  =  21  sin  (arctan  (A/e))  is  the  dipole  moment  of  the  quantum  dot,  T\  increases  for  larger 


bias. 

Next,  we  present  results  from  the  relaxation  rate  calculation  for  a  DQD  in  a  two- 
dimensional  quantum  well.  In  this  treatment  we  label  the  z-axis  as  the  vertical  confinement 
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Figure  7.4:  Energy  relaxation  time  T\  vs.  the  distance  from  the  metallic  gate  to  the  DQDs 
z  for  dot  geometry  of  d  =  /  =  30  nm  and  l  =  30  nm  (dashed  and  dash-dotted  blue  lines) 
and  60  nm  (solid  and  dashed  black  lines)  at  0  K,  cj/27t  =  1  GHz,  and  e  =  0.  Solid  and 
dash-dotted  lines  are  T\  times  for  the  exact  form  of  the  interaction  Hamiltonian,  whereas 
dashed  lines  are  for  the  dipole  form  of  the  interaction. 


direction  and  do  not  consider  excitations  along  the  z-direction.  We  model  the  confining 
potential  by  a  symmetric  double  rectangular  well  in  2D  and  numerically  compute  the  lowest 
two  eigenenergies  and  wavefunctions.  We  then  compute  the  electron  relaxation  rate  between 
these  two  lowest  states.  The  results  are  qualitatively  similar  to  the  one-dimensional  case  and 


are  shown  in  Fig.  7.4,  where  we  plot  the  energy  relaxation  time  T\  vs.  the  distance  2  from 
the  metallic  gate  to  the  DQD.  In  this  plot,  we  choose  the  size  of  the  dot  in  the  x-direction 
d  =  30  nm,  the  size  in  the  y-direction  /  =  30  nm  and  half  the  separation  between  the  dots  l 
to  be  30  nm  (dashed  and  dash-dotted  blue  lines)  and  60  nm  (solid  and  dashed  black  lines). 
We  find  that  for  l  =  30  nm  and  z  =  90  nm,  T\  is  4  ys  while  for  l  =  60  nm,  T\  is  roughly 
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1.6  fj,s.  These  relaxation  times  are  somewhat  longer  than  the  experimentally  reported  value 
of  Ti  =  20  ns  in  DQD-based  charge  qubits.  |1141 1137]  We  note  that  the  relaxation  rate  for  a 
two-dimensional  DQD  is  shorter  than  for  a  one-dimensional  DQD  of  comparable  geometry 
by  about  a  factor  of  5.  A  two-dimensional  DQD  is  able  to  couple  to  obliquely  oriented 
wavevectors  in  addition  to  those  which  point  in  the  direction  of  separation  between  the 
dots,  and  this  can  be  reasonably  expected  to  enhance  relaxation  by  a  geometric  factor  of 
order  unity. 


7.3  Spin  qubit 


We  now  focus  on  the  calculation  of  the  relaxation  rate  for  a  single  electron  in  a  spin  qubit 
realized  in  a  single  quantum  dot.  Here  the  system  Hs  and  the  interaction  Hint  Hamiltonians 
are  given  by 


Hs  = -gnBa  ■  Bo/2  (7.17) 

Hint  =  ~9Hb  o  ■  B(f,  t)  (7.18) 


where  a  is  the  vector  of  Pauli  matrices,  g  is  the  gyrometric  factor  of  the  trapped  electron  in 
a  quantum  dot,  gs  is  the  Bohr  magneton,  Bo  is  the  externally  applied  magnetic  field  and 
B(r,t )  is  the  fluctuating  EWJN  field.  The  rate  of  spin  flip  from  excited  |  |)  to  ground  |  4) 
at  T  =  0  K  can  be  obtained  from  the  Fermi  golden  rule 

Ti,s  =  -jr  J  d3r  j  d3f'Mr>s(r)Mr}S(f,)eijkeij'k'lmXkk'{^^,^)njnj'  ■  (7-19) 


where  repeated  indices  are  summed  over,  and  nj  are  the  components  of  a  unit  vector  h  in 
the  direction  of  Bq.  The  effect  of  finite  temperature  on  the  transition  rates  is  the  same  as 


for  charge  qubits,  as  shown  in  Eq.  7.11  The  magnetic  response  function  or  susceptibility 


xfj  and  the  matrix  element  Mr_s(r)  are 


B  /  ->  ->/  \ 

Xij{r,r  ,w)  = 


£ikm.£jnp 

e0c2 


dkdnDmp(r,  r  ,  uj) 


Mr,s(f)  =  g^B |^o(70|2  • 


(7.20) 

(7.21) 
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Figure  7.5:  Energy  relaxation  time  T\  vs.  the  distance  from  the  metallic  gate  to  quantum 
dot  z  for  spin  qubit  (in  Si)  single  dot  geometry  of  d  =  60  nm  at  0  K  temperature  and 
gfiBBo/2Trh  =  50  GHz.  The  solid  line  represent  T\  time  for  the  exact  form  of  the  interaction 
Hamiltonian,  whereas  the  dashed  line  is  for  the  dipole  form  of  the  interaction. 


Here  the  spin  qubit  frequency  u  =  ggsBo/H  and  ipo (r)  is  the  spatial  part  of  the  the 
ground  state  wave  function  of  the  spin  qubit.  Eq.  |7. 19| is  a  generalization  beyond  the  dipole 
approximation  of  the  simpler  expression  [33] : 


ri-  = 


9  VBaB 0 

167r/i2eQC4z’ 


(7.22) 


which  has  been  obtained  by  using  the  quasistatic  limit  for  the  Green’s  function  7.20,  and 


assuming  it  is  constant  over  the  spatial  extent  of  the  qubit.  Eq.  7.22  also  assumes  the 
external  magnetic  field  Ho  points  in  the  ^-direction. 

A  plot  of  energy  relaxation  time  T\  vs.  the  distance  from  the  metallic  gate  z  for  a  spin 


qubit  is  displayed  in  Fig.  7.5  Here  we  consider  a  single  quantum  dot  of  diameter  d  =  60 
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nm  and  approximate  the  ground  state  spatial  wave  function  of  the  spin  qubit  by  the  ground 
state  wave  function  of  a  harmonic  potential.  We  assume  the  Zeeman  splitting  between  spin 
states  is  50  GHz,  typical  of  experiments  in  spin  qubits.  [135]  The  solid  line  is  the  T\  time 
obtained  using  the  non-local  magnetic  spectral  density  while  the  dotted  line  is  obtained  for 
a  local  spectral  density,  which  diverges  as  ~  \/z  as  one  approaches  the  metallic  gate.  The 
reason  for  saturation  of  the  T\  time  at  smaller  distances  is  similar  to  the  case  for  charge 
qubits.  There  is  a  slight  distinction  in  that  the  spin  case  involves  a  spatially  extended  dipole 
interaction,  as  opposed  to  the  charge  case  which  involves  qenuine  quadrupole  and  higher 
multipole  contributions.  This  distinction  is  largely  technical,  however,  and  a  saturation  of 
T\  as  2  — )•  0  is  observed  in  both  cases. 

We  find  that  the  T\  time  for  a  spin  qubit  in  a  GaAs  quantum  dot  with  an  external 
magnetic  field  of  10  T  and  z  =  90  nm  is  roughly  1  s  which  is  larger  than  the  experimentally 
reported  value  (135]  of  0.55  ms,  and  generally  EWJN  does  not  seem  to  be  the  dominant 
source  of  decoherence  for  semiconductor  devices  in  large  magnetic  fields.  GaAs  has  a  strong 
spin-orbit  interaction  (SOI),  which  mixes  the  Zeeman-split  spin  states  with  orbitally  excited 
states.  Spin  relaxation  can  then  occur  via  coupling  of  the  qubit  to  piezoelectric  phonon 
noise  in  the  2DEG  layer.  The  relaxation  rate  from  this  mechanism  scales  as  B5  and  is  the 
dominant  pathway  for  spin  relaxation  at  large  external  magnetic  fields  B  >  IT  [1381 1139] . 
Additionally,  Marquardt  and  Abalmassov  calculate  relaxation  of  spin  qubits  from 

electric  EWJN  via  the  SOI.  Again,  mixing  of  the  charge  and  spin  states  via  the  SOI  allows 
spin  relaxation  to  be  induced  from  electric  field  fluctuations.  They  estimate  the  power 
spectrum  of  the  Johnson  noise  using  a  lumped  circuit  model  and  found  a  B3  dependence 
of  the  relaxation  rate.  Our  treatment  involves  a  direct  coupling  of  the  fluctuating  magnetic 
field  from  the  top  gates  with  the  spin  states,  and  our  rate  scales  linearly  with  the  magnetic 
field.  We  therefore  expect  our  relaxation  pathway  to  dominate  at  low  magnetic  fields,  and 
indeed  while  we  predict  a  much  slower  relaxation  rate  than  measured  by  Arnasha  et.  al  [138] 
for  B  ~  7T,  at  .B  =  1  T  our  results  predict  T\  ~  5  s  which  is  comparable  to  their  measured 
value  of  T\  =  1  s.  Additionally,  in  Si  quantum  dots  with  a  2  T  external  magnetic  field  and 
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z  =  50  nm,  we  predict  a  T\  time  of  15  ms  which  is  smaller  than  the  experimentally  reported 
value  of  40  ms  HM|.  However,  it  must  be  kept  in  mind  that  we  have  so  far  considered 
the  simpler  top  gate  geometry  of  a  conducting  half-space  rather  than  the  thin  layer  of 
finger  gates  used  in  these  experiments.  In  the  next  section  we  address  modifications  to  our 
calculations  that  we  expect  from  more  realistic  gate  geometry. 


7.4  Thin  metallic  gates 


A  conducting  half-space  is  an  analytically  convenient  gate  geometry,  but  a  poor  approxi¬ 
mation  to  the  thin  top  gates  commonly  used  in  semiconductor  devices.  In  this  section  we 
present  an  exact  treatment  of  the  behavior  of  EWJN  in  the  vicinity  of  a  metallic  film  of 
finite  thickness.  Changing  the  half-space  to  a  thin  film  affects  EWJN  by  modifying  the 
reflection  coefficients  rs  and  rp.  The  power  spectrum  of  the  resultant  EWJN  is  obtained  by 
substituting  these  modified  reflection  coefficients  into  the  photon  Green’s  function  7.14|  and 
the  relaxation  time  of,  e.g.  a  charge  qubit,  is  obtained  by  plugging  Eq.  |7.14  into  Eqs.  7J5 
and  17.131  The  modified  reflection  coefficients  for  a  him  of  thickness  a  take  the  form 


rp(k,u,a)  = 


rs(k,  uj,  a)  = 


e2q2  -  q\ 


q2  +  e2q2  +  2iqq\ecot(qia) 
Q2  ~  Qi 


(7.23) 

(7.24) 


q2  +  q\  +  2,iqq\  cot (<71  a) 

They  differ  significantly  from  the  half-space  result  only  when  the  thickness  a  is  of  the  order 
or  smaller  than  the  skin  depth  <5,  and  they  reduce  to  the  half-space  result  for  a><5.  A 


derivation  of  Eqs.  7.23  and  7.24  is  given  in  the  Appendix.  Eqs.  7.23  and  7.24  are  exact,  but 
for  a  good  conductor  they  can  be  cast  into  a  simpler  approximate  form 


rp(k,  u>,  a)  «  (  1  +  ^  cot  (qia) 


rs(k,u,a)  «  -  (  1  - 


2c2q\k 


eu^ 


cot(<7ia) 


-1 


(7.25) 

(7.26) 


These  expressions  have  been  obtained  by  expanding  Eqs.  7.23  and  7.24  for  large  imag¬ 


inary  e  and  then  taking  the  quasistatic  approximation  q  — - >  ik.  The  first  approximation 
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Figure  7.6:  Ratio  of  energy  relaxation  time  T\  from  conducting  film  to  T\  time  from  half¬ 
space  vs.  thickness  of  the  film  a  for  a  DQD  charge  qubit  in  one  dimension  with  dot  geometry 
d  =  30  nm  and  l  =  60  nm  at  0  K  temperature.  We  take  the  exact  form  of  the  interaction 
Hamiltonian.  The  distance  z  from  the  film  or  half-space  is  chosen  as  follows:  z  =  10  nm 
(black  dash-dotted  line),  z  =  50  nm  (blue  dashed  line)  and  z  =  150  nm  (solid  red  line). 


Other  parameters  are  the  same  as  in  Fig.  7.2 


is  extremely  accurate  for  copper  near  zero  temperature  and  the  second  is  accurate  for  all 
distances  z  such  that  EWJN  is  appreciably  enhanced  above  blackbody  radiation  [33] .  The 


remarkable  feature  of  Eqs.  7.25  and  7.26  is  that  they  show  the  strength  of  the  fluctuating 
fields  outside  the  film  are  actually  amplified  relative  to  the  half-space  result.  This  can  be 
understood  by  analogy  to  the  behavior  of  a  particle  trapped  in  a  finite  one-dimensional 
potential  well.  For  a  given  width  of  the  well,  the  wavefunction  will  have  an  exponentially 
decaying  tail  in  the  forbidden  region.  As  the  confinement  is  increased,  the  particle  will 
be  squeezed  and  its  wavefunction  will  leak  farther  into  the  forbidden  region.  It  will  be 
interesting  to  see  if  this  enhancement  is  observable  in  the  Casimir  attraction  between  2  thin 
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conducting  plates. 

Using  the  modified  expression  for  the  reflection  coefficients,  we  compute  the  T\  time  of 


a  DQD  charge  qubit  in  one  dimension  due  to  the  metallic  film.  In  Fig.  7.6,  we  plot  the  ratio 
of  the  T\  time  obtained  for  the  film  to  the  time  computed  for  the  metallic  half-space  as  a 
function  of  the  film  thickness.  We  take  the  exact  form  of  the  interaction  Hamiltonian  for 
a  variety  of  distances  from  the  gate.  We  find  that  for  distance  z  >  a,  the  relaxation  time 
due  to  the  film  can  be  reduced  by  over  an  order  of  magnitude  relative  to  the  half-space.  It 
converges  to  the  half-space  result  as  z  becomes  smaller  than  the  thickness  of  the  film. 

Common  semiconductor  qubit  architectures  employ  thin  finger-shaped  top  gates  which 
are  more  sparse  than  the  films  considered  here.  An  exact  treatment  of  EWJN  from  a 
detailed  finger  gate  geometry  would  be  prohibitively  difficult,  but  we  expect  to  a  reasonable 
approximation  that  EWJN  from  finger  gates  will  be  reduced  by  a  factor  of  the  fraction  of 
the  top  gate  layer  that  is  not  composed  of  metal.  Our  results  should  then  overestimate 
the  relaxation  rate  by  a  geometric  factor.  We  note  however  that  newer  accumulation¬ 
mode  architectures  employ  a  second  top  gate  above  the  confinement  top  gates  (136} .  These 
accumulation  gates  are  solid  sheets  and  are  typically  around  100  nm  from  the  qubit,  so  our 
treatment  should  accurately  describe  their  contribution  to  relaxation. 


7.5  Discussion 

In  conclusion,  we  have  presented  a  detailed  study  of  the  effect  of  evanescent-wave  Johnson 
noise  on  energy  relaxation  of  quantum  dots  beyond  the  dipole  approximation.  We  have 
noted  that  previous  studies  of  charge  and  spin  qubits  which  use  the  dipole  approximation 
allow  contribution  from  infinitely  large  components  of  the  photon  wavevector  leading  to 
overestimation  and  divergence  of  the  energy  relaxation  rate  as  z  — >  0.  We  have  demon¬ 
strated  that  it  is  possible  to  remedy  this  spurious  divergence  by  taking  into  account  the 
finite  size  of  the  quantum  dot.  While  a  non-local  permittivity  of  the  surface  metal  will 
remove  the  divergence  in  the  field  fluctuations  at  the  boundary,  we  have  shown  that  the 
finite  size  of  the  dot  provides  an  alternative  normalization  mechanism  by  enforcing  a  finite 
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cutoff  in  the  magnitude  of  the  contributing  wavevector.  In  addition,  we  have  derived  a  novel 
enhancement  of  the  EWJN  field  fluctuations  that  occurs  outside  a  metallic  him,  relative  to 
the  held  outside  a  metallic  half-space. 

This  chapter  has  focused  exclusively  on  relaxation,  though  we  expect  dephasing  times 
from  EWJN  to  be  of  comparable  magnitude.  The  power  spectrum  of  EWJN  is  linear  in 
oj,  and  this  will  suppress  contribution  from  the  small  frequency  part  of  the  electromagnetic 
spectrum,  which  typically  enhances  dephasing  rates.  While  the  temperature  dependence  of 
the  relaxation  rate  is  simply  given  by  the  Planck  function,  we  do  expect  a  more  non-trivial 
temperature  dependence  of  the  dephasing  rate. 

Of  particular  interest  are  experimental  signatures  of  EWJN-induced  relaxation.  No¬ 
tably,  at  zero  temperature  the  charge  relaxation  rate  scales  linearly  with  the  qubit  transi¬ 
tion  frequency  and  as  the  inverse  cubic  power  of  the  distance  between  the  qubit  and  the 
metallic  top  gates.  The  zero  temperature  spin  relaxation  rate  scales  linearly  with  the  ex¬ 
ternal  magnetic  field  and  inversely  with  the  distance  to  the  gates.  Our  results  indicate 
that  EWJN  from  the  metallic  top  gate  is  not  a  dominant  source  of  relaxation  in  charge 
qubits,  but  can  be  the  dominant  noise  source  for  energy  relaxation  in  spin  qubits  held  at 
low  external  magnetic  field. 
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Chapter  8 

Conclusions 


In  this  thesis,  we  have  considered  several  aspects  of  superconducting  and  semiconducting 
qubits,  exploring  mainly  the  effect  of  decoherence  on  these  devices.  In  the  second  chap¬ 
ter  [2j  we  discuss  different  noise  sources  that  plague  semiconducting  and  superconducting 
qubits  and  present  different  methods  we  can  use  to  study  the  dynamics  of  open  quantum 
systems.  In  the  third  chapter  [3j  we  investigate  the  effect  of  the  Ohmic  environment  on 
optimally  controlled  phase  qubit.  Here  we  compare  possible  choices  of  microwave  pulses 
for  NOT  gate  operation.  Specifically,  we  consider  three  options:  single-quadrature  pulses 
and  two-quadrature  microwave  (control)  pulses  with  both  variable  and  fixed  frequency. 
Two-quadrature  pulses  led  to  significant  suppression  of  the  gate  error  compared  to  single¬ 
quadrature  pulses.  However,  the  presence  of  dissipative  environment  increased  the  gate 
error  even  for  two-quadrature  pulses  significantly  above  the  required  threshold  for  fault  tol¬ 
erant  quantum  computation  in  currently  available  phase  qubits.  We  further  investigate  how 
the  environmental  coupling  strength  affects  the  gate  error  and  found  that  an  improvement 
of  the  qubit  relaxation  time  is  crucial  for  effectiveness  of  the  DRAG  pulses.  We  determine 
that  two-quadrature  pulses  with  fixed  driving  frequency  suppress  the  gate  error  below  the 
required  threshold  for  a  reasonable  gate  time  of  7  ns,  but  for  qubits  with  the  relaxation 
time  ten  times  longer  than  the  currently  observed  relaxation  time. 

In  chapter  [4j  we  present  a  model  to  determine  the  quantum  efficiency  of  a  microwave 
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photon  detector  based  on  a  current-biased  JJ.  We  demonstrate  that  the  efficiency  to  detect 
single  photon  loaded  in  a  cavity  is  determined  by  relaxation  and  tunneling  rates.  We  also 
determine  that  the  bandwidth  of  the  detector  is  characterized  by  the  sum  of  the  vacuum 
Rabi  frequency  and  the  broadening  of  the  first  excited  state  of  the  JJ  due  to  tunneling  and 
relaxation  processes.  Our  results  indicate  that  for  currently  used  JJ  photon  detectors,  the 
quantum  efficiency  of  around  80%  can  be  achieved  even  for  modest  junction  parameters. 
In  addition,  we  present  alternative  schemes  for  photon  detection  and  present  a  systematic 
method  to  compute  the  power  absorbed  by  the  detector  in  such  schemes. 

In  chapter  [5j  we  discuss  the  effect  of  decoherence  on  the  Berry  curvature  and  present 
results  which  demonstrate  that  the  curvature  is  immune  to  decoherence  in  non-adiabatic 
measurement  protocols.  Surprisingly,  decoherence  forces  magnetization  to  respond  linearly 
to  the  drive  velocity  in  case  of  the  sinusoidal  drive,  thereby  assisting  in  the  measurement 
of  the  Berry  curvature.  Furthermore,  we  analyze  the  prospect  of  measuring  the  Berry 
curvature  in  superconducting  qubits.  Here  we  investigate  the  effect  of  the  third  level  in 
phase  or  transmon  qubits  on  the  measurement  of  the  Berry  curvature. 

In  the  remaining  chapters,  we  discuss  the  origin  and  effect  of  Johnson  noise  in  semicon¬ 
ducting  qubits.  In  chapter  [6j  we  identify  evanescent  wave  Johnson  noise  as  an  important 
source  of  decoherence  in  both  spin  and  charge  qubits.  Here  we  carry  out  analysis  within 
the  dipole  approximation  of  the  qubit-field  interaction.  Finally,  in  chapter  [7j  we  continue 
to  discuss  the  effect  of  evanescent  wave  Johnson  noise  on  extended  qubit  geometries.  Here 
we  go  beyond  the  dipole  approximation  and  present  result  for  relaxation  rates  for  both  spin 
and  charge  qubits.  We  also  consider  the  finite  thickness  nature  of  metallic  gates  and  present 
enhancement  of  electric  evanescent-wave  Johnson  noise  compared  to  field  fluctuations  from 
extended  gate  geometry.  These  results  serve  as  an  important  constraint  on  future  device 


architectures. 
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Appendix  A 


Supplemental  information  for 


chapter  p 


A.l  Hamiltonian  of  a  flux-biased  phase  qubit 

In  this  appendix,  we  detail  the  derivation  of  the  Hamiltonian  of  a  flux-biased  phase  qubit 


used  in  Eq.  3.2  The  phase  qubit  consists  of  a  superconducting  loop  interrupted  by  a 


Josephson  junction.  The  Josephson  junction  consists  of  an  insulating  barrier  sandwiched 
between  two  superconducting  layers.  In  1962,  B.D.  Josephson  predicted  that  super  current 
would  tunnel  through  the  insulating  barrier  even  in  the  absence  of  applied  voltage  [14DJ. 
He  showed  that  magnitude  of  super  current  Is  depends  only  on  the  phase  difference  5 
(Josephson  phase)  between  the  order  parameters  of  two  superconducting  electrodes.  This 
is  called  the  first  Josephson  equation  which  reads 


Is  =  I0  sin  5  , 


(A.l) 


where  Iq  is  the  critical  current  of  the  junction.  The  second  Josephson  equation  relates  time 
evolution  of  the  Josephson  phase  with  the  voltage  drop  V  across  the  junction,  that  is, 


dd  V 

—  =  2e— 
dt  h 


(A.2) 
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where  h  and  e  are  Planck’s  constant  and  electronic  charge,  respectively.  The  dynamics  of 
the  phase  qubit  is  best  understood  within  a  phenomenological  resistively  and  capacitively 
shunted  junction  model  (RCSJ)  of  a  Josephson  junction  j  1 41] .  In  a  flux-biased  phase  qubit, 
biasing  is  done  by  coupling  external  flux  to  the  superconducting  loop,  instead  of  directly 
injecting  the  bias  current.  This  technique  helps  prevent,  to  some  extent,  noise  in  the  bias 
circuitry  to  couple  to  the  qubit. 

Applying  Kirchoff’s  law  to  the  RCSJ  circuit  model  along  with  Josephson  equations 
and  phase  quantization  condition  in  a  superconductor,  we  obtain  the  following  equation  of 
motion, 


$0\2  d2S 


<f>r 


C  ^  -770+  77W7+  7^  WV 


1  d5  /+o\  d 


2n  J  dt 2  \ 2vr y  R  dt  \2ir  J  d5 


$0 

47tL 


6- 


27r$6y 

+0  J 


—  /q  cos  6 


=  0 


(A.3) 


where  C,  R  and  L  are  junction  capacitance,  shunt  resistance  and  loop  inductance,  respec¬ 
tively.  The  parameters  <!>;,  and  4>o  =  are  external  bias  flux  applied  to  the  phase  qubit 
and  superconducting  flux  quantum,  respectively.  Eq.  A. 3| represents  the  equation  of  motion 
of  a  fictitious  particle  with  mass  C  under  the  potential  energy  given  by 


- I 


+0 

47tL 


5- 


27r+bV 
+0  J 


—  Jo  cos  5 


(A.4) 


The  parameter  can  be  tuned  into  appropriately  to  change  the  potential  energy  landscape 
of  a  flux-biased  phase  qubit.  For  instance,  when  +5  =  4>0)  the  potential  energy  takes  the 
shape  of  a  symmetric  double  well.  Typically,  for  the  qubit  operation,  this  parameter  is 
tuned  in  between  0.95  —  0.98  of  4>c,  where  4>c  is  a  critical  flux  of  the  qubit  for  which  the 
second  minimum  in  U(5)  disappears.  With  this  expression  for  the  potential  energy,  the 
Hamiltonian  of  a  flux-biased  phase  qubit  is  given  by: 


TT  _  £ 0 

q~  2 C  +  2ir 


+0 

47tL 


6- 


2i r3>by 


—  /q  cos  5 


(A.5) 


+0  7 

Here  Q  corresponds  to  charge  across  the  Josephson  junction.  Until  now,  charge  and  phase 
of  the  JJ  were  treated  as  classical  variables.  In  order  to  study  the  quantum  properties 
of  the  phase  qubit,  these  variables  must  be  promoted  to  operators.  Hence,  for  a  phase 
qubit,  Q  and  6  are  treated  as  a  conjugate  momentum  and  position  operators,  respectively 
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that  satisfy  the  commutation  relationship  Q,5  =  i2e.  In  S— representation,  the  operator 
5^-5  and  Q  — >  —2ie-j^.  Therefore,  in  this  representation,  time-independent  Schrodinger’s 


equation  takes  the  following  form: 


— 2e2  d 2V>  i  $0 
C  +  2n 


2vrTfe\2 


5 - - —  -  I0  cos  5 

4irL  V  4>0  ) 


i/>(8)  =  Eip(S) 


(A.6) 


In  the  following  section,  we  will  outline  a  numerical  method  to  solve  this  time-independent 
Schrodinger’s  equation. 


A. 2  Numerical  solution  of  the  Schrodinger’s  equation 

Here  we  will  detail  the  algorithm  used  to  obtain  the  numerical  solution  of  the  time- 
independent  Schrodinger’s  equation.  The  Schrodinger’s  equation  is  discretized  over  a  one 
dimensional  spatial  interval  [<5o,  <5/v]  of  IV +2  points  such  that  5j  =  So+je,  where  e  is  the  step 
size,  and  j  E  [0,  N] .  The  boundary  conditions  on  wave  functions  over  this  interval  are  such 
that  V’nfA)  =  iPii(8n)  =  0.  On  the  discretized  interval,  time-independent  Schrodinger’s 
equation  of  the  qubit  system  can  be  written  as 


Atpn{6j)  +  U(5j)^n{5j)  =  Enijjn{8j)  ,  where 


U(Sj)  =  A1 


27r$by 
$0  ) 


—  /5l  cos (Sj 


(A.7) 

(A.8) 


Here,  A  =  -1,  En  =  -En/A0:  A0  =  -2 e2/C,  A1  =  - 4A® and  We  have 

rescaled  the  potential  and  kinetic  energy  terms  for  numerical  stability.  The  derivative  in 
Eq.  A.7  is  approximated  using  the  finite  differences  scheme  as  follows: 

-ipn(5j+2)  +  16ipn(6j+i)  -  30 ipn(5j)  +  16tpn(Sj-i)  -  2) 


12e2 


(A.9) 


The  error  in  this  approximation  scheme  is  of  the  order  of  e4.  With  this  approximation  for 


the  derivatives,  evaluating  Eq.  A.7  at  each  discretized  point,  we  obtain  a  system  of  N  —  1 
linear  equations,  which  can  be  written  in  a  matrix  form.  The  eigen  functions  and  eigenvalues 
of  this  matrix  are  the  desired  energy  eigenvalues  and  wave  functions  of  the  Schrodinger’s 
equation. 


105 


Appendix  B 

Supplemental  information  for 
chapter 


B.l  Derivation  of  Green’s  tensor  for  half  space 

In  this  appendix,  we  will  derive  Green’s  function  for  the  metallic  half  space.  We  begin 
by  calculating  Green’s  function  for  parallel  plate  geometry  in  three  spatial  dimensions. 
Our  starting  point  for  the  calculation  of  Green’s  function  in  the  free  space  between  two 
conducting  half  spaces  of  relative  permittivities  t\  (r,ui)  and  £2 (r,  w)  in  three  dimensions 
will  be: 

(~5ij  (v2  +  Uj2e(^Uj))  +  didj)  Dlk(f,f ')  =  53(r  -  r')5jk  (B.l) 

The  dependence  of  Dik(r,r')  on  u  is  suppressed  to  simplify  notation.  We  take  the  bound¬ 
aries  of  the  two  conducting  half  spaces  to  be  located  at  z  =  0  and  z  =  l,  with  vacuum  for 
0  <  z  <  l.  Because  the  geometry  is  translationally  invariant  in  the  x  and  y  directions,  we 
expand  Dik  in  a  Fourier  series  for  x  —  x'  and  y  —  y': 


(B.2) 


Dij(k,z,z’)  =  I  e~i%r*WD(r,r')d2\ 

0  =  J  z,  z')d2k 


(B.3) 
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Here  In  =  rji  —  r'li  =  (x  —  x',y  —  y').  After  B.3  is  plugged  into  B.l,  the  resulting  set  of 
differential  equations  will  significantly  decouple  if  we  choose  k  to  lie  along  the  x-axis,  i.e. , 
choose  k  =  (k,  0).  Once  this  is  done,  B.l  reduces  to 


q2  ur‘e\  ~  d 

7^2  +  )  Dxj(k,z,z')  -  ik—Dzj(k,  z,  z')  =  5xj5(z  -  z') 


7  2  u  e 

-k2  + 


dz' 

d 


Dzj(k,  z,  z')  -  ik—Dxj(k,  z,  z')  =  Szj5(z  -  z') 
^—k  +  ^2  +  ^  Dyj{k,  Z,Z  )  =  Syjd(z  —  Z  ) 

Where  j  =  x,  y,  z.  Combining  these  equations  we  arrive  at 
(Jp2  +  ^  Dyy{k,  z,  z)  =  5(z  -  z') 

/  f)2  \  n2  2 

{^2  +  Q2)  Dxx(k,  Z,  z)  =  -  A 

^  zA/  3  ~ 

Dzx{k">  Z”>  Z  )  =  ~^2  %•)  z  ) 

~  ik  3  ~  1 

Dzz(k,z,z')  =  -^—Dxz(k,z,z')  +  ~^6(z  -  z') 
qz  oz  qz 


We  solve  first  for  Dyy]  its  boundary  value  problem 


is 


Dyy(k,  Z,Z)  -  < 


Ae  4912  ,  z  <  0 

Cie-i92  +  C2eiqz  +  ,  0  <z<l 

Beiq2(z-l)  ^  l  <  z 


(B.4) 

(B.5) 

(B.6) 


(B.7) 

(B.8) 

(B.9) 

(B.10) 


(B.ll) 


2  _  2/2  7  2  2  _  2/2  z2  2_  2/2  z2 

q  —  uj  /c  —  k  ,  ^  =  eio;  /c  —  k  ,  g2  =  e2^  /c  —  k 


(B.12) 


Since  and  dDyy/dz  must  be  continuous  at  the  interfaces  at  z  =  0,/,  we  can  now 
determine  the  constants  C\  and  C2  to  be 


Ci  = 


2qA 


0iqz 


Q  +  Ql  ,-iqz' 


q-qi 


rjn  —  *  ( p-iqz’  I  9  iq(z’-2l) 

2  2qA\  +q~q2 

A  = ,  _  {q  +  qi){q  +  q2)  -2iqi 
-  (q-qi)(q-Q2) 


(B.13) 

(B.14) 

(B.15) 
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We  plug  these  back  into  Dyv  and  obtain 


Dyy(k,  Z,  z')  = 


yy 

—ci  ( v 


+  y/ei  —  v? 


2tovA  \v  —  C\  —  u2 


e-iuv(z+z')/c  _|_  v  +  V£2  ~  u~  riu>v(z+z' -21)/ c 


—  Ve2  — 


u- 


+  (2cos(o;u[z  -  z']/c)/ A  - 

Zuy  c/  \ 

where  we  have  introduced  the  definitions  q  =  to\J  1  —  u2/c  =  tov/c,  u  =  ck/u,  and 

\/l  —  u2  ,  0  <  u  <  1 


(B.16) 


v  =  < 


i\/u2  —  1  ,  u  >  1 


(B.17) 


A  similar  calculation  will  yield  -Dxa;,  although  the  boundary  conditions  for  this  problem  are 
and  e/q2(dDxx/dz)  continuous  across  the  boundaries  at  z  =  0,1.  The  result  is 


Dxx  (k,  z,  z  )  — 


CIV  I  €\V 


2toA 
civ 


+  Vei  -  u2  r_iuJV(z+z')/c  £2V_ 


+  Ve  2  — 


l  \eiu  —  —  u2  C2U  —  y/e2  —  u 


u  Ciujy(z+z'-2l)/c 


+  ^  ( 2  cos(u;u[z  -  ^]/c)/Ai  -  e^-2'l/c 
zee?  V 

A  =i_  (ew  +  Ve^~^)(e2y  +  Ve2  -  u2)  p_2w/c 
(eju  -  \/ei  —  u2)(e2V  -  V«2  -  u2) 


From  Eqs. 


B.9 


and 


B.10 


Dzz(k,  z,  z')  is  determined  to  be 
k2  ~ 

Dzz(k,  z,z  )  —  ^)Dxxi_  2  j  -2) 


(B.18) 

(B.19) 

(B.20) 


Next,  we  back  Fourier  transform  to  obtain  these  quantities  in  the  real  space.  Eqs.  |BT6l 


and 


B.18 


were  obtained  assuming  that  k  lies  entirely  in  the  x  direction.  Because  of  spatial 


invariance  in  the  x  —  y  plane,  B.18  would  be  equal  to  Dyy(k,  z)  if  k  were  chosen  to  point  in 
the  y-direction.  Eqs.  |B.16  and  B.18  are  then  the  contributions  into  the  kx  and  ky  Fourier 
components  of  Dyy(f,r',to),  respectively.  We  can  therefore  construct 

Dyy{f\\  ,z,to)  =  Dxx(r\\  ,z,to)  =  j  (A/2/(fc’  z> u)  sin2  0  +  Dxx(k,  z,  to)  cos2  d^j  d2k 

(B.21) 

The  exponential  term  vanishes  in  the  limit  r—tr1,  and  the  angular  integration  simply  gives 
a  factor  of  ir.  The  trigonometric  factors  are  squared  because  the  Green’s  function  is  related 
to  the  expectation  value  of  the  square  of  the  electric  field. 


Dyy{fu ,  z,  to)  =  Dxx(f\\ ,  z,  to)  = 


47T 


Dyy{k ,  z,  to)  +  Dxx(k ,  z,  to) )  kdk 


(B.22) 
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Meanwhile  the  integral  expression  for  Dzz(rn,  z,u)  may  be  trivially  derived  and  is  equal  to 

1  k 3  ~ 

Dzz(f\\,  z,uj)  =  —  /  —Dxx(k,z,z')dk  (B.23) 

^7r  Jo  Q 


The  case  of  a  half  space  can  be  obtained  from  B.22  by  taking  the  limit  l  — >  oo  in  the 
evanescent  range  of  q  (when  it  is  imaginary): 

72 


~  CX  l  V 

£&  D»(t’ Z)  =  2^~v  U  + 

eiv  -  \J £]  - 


V^l  ^liuvz/c  _|_  ^ 


lim  Dxx(k,z)  | k>U]/c  = 


—civ 


l — ^OO 


2  u  +  \fe\  — 


^_^2iuivz/c  _  ^ 


(B.24) 

(B.25) 


This  leads  to  the  following  expression  for  D 


hs 


DivS)(n,z,u;)  =  D^xs\rn,z,uj)  = 


xx, yy 

00  udu 


yy 


87 rc  J o  v 


(e2i“vz/c  ( rs(u )  +  ( u 2  -  1  )rp(u))  +2  -  u2)  , 

(B.26) 


where  we  have  introduced 

rP  («)  = 


ev  —  V e  —  u2 


and  rs (it)  = 


u  —  \J  e  —  u2 


ev  +  \J  e  —  u2  v  +  y/  e  —  u2 


(B.27) 


rp  is  the  reflection  coefficient  for  the  part  of  the  electric  field  perpendicular  to  the  plane  of 
incidence  and  rs  is  the  reflection  coefficient  for  the  part  of  the  electric  field  parallel  to  the 
plane  of  incidence.  The  imaginary  part  of  the  response  function  (near-field  only)  is  given 
by 

Im  Xxx(r,  r*,  u)  =  Im  D^xs)  (rj( ,  z,  u) 
e0c2  11 


LV 


87TC0  C3 


Re 


J°°  ^  e2  ^/cjrj(u)  +  (u2  _  i)rp(u))  (B.28) 
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Appendix  C 


Supplemental  information  for 
chapter  [7] 


C.l  Derivation  of  generalized  Fermi’s  golden  rule 


In  this  appendix,  we  will  derive  the  generalized  Fermi’s  golden  rule  used  in  Eq.|7.5[  Consider 
a  system-bath  interaction  Hamiltonian 


Hsb  =  ~ 


2m 


^A(r)  ■  p  +  p  ■  A(r)^ 


(C.l) 


Since  [p,  A(r)]  =  —iKV  ■  A[r),  the  interaction  Hamiltonian  can  be  written  as:  Hsb  = 
—  ^ rA(f)  ■  p  +  •  A{f).  The  relaxation  rate  for  |— )  — >  |+)  transition  at  zero  tem¬ 

perature  can  be  calculated  using  the  Fermi’s  golden  rule.  To  calculate  the  total  relaxation 
rate,  we  have  to  sum  over  all  initial  reservoir  states  weighted  with  the  probability  to  find 
these  states  and  over  all  final  reservoir  states,  namely  : 


r ^  [(R' ,tl)+\Ai{r)pi  -  ih/2diAi(r)\i/j-,  R)* (R' ,ilj+\Aj(r')pj 
ij  R'R 

-  ih/2djAj(r')\ilj-,R)  +  i^  j  P(R)  <5(e_  +  eR  -  e+  -  eR>) , 


(C.2) 


where  e±  are  energies  of  |±)  states,  |f?)(|f?/))  are  reservoir  initial  (final)  states  respectively, 
and  P(R)  is  the  probability  to  find  the  reservoir  in  initial  state  R.  Next,  rewrite  the  delta 
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function  in  terms  of  its  Fourier  transform 


5{e- +  £r  — e+ —  £} j/)  =  — :  [  d(t  —  t')  elujlyt  4')  where  w  =  (e_  —  £+)/h, 

2irn  J 

and  introduce  time-dependent  Heisenberg  operator  A{f ,  t)  to  obtain 


1 


.  =  - 


2  Lm*2/i2 


/  dr'tp+(r)(ih)di,ip*_(r)  /  di* 'ip+(r')  (—ih)dj'ip-(r 


x^P(i2)  /(i2,|Ai(r,t)|i2)*(i2,|AJ-(r^,0|i2>e<w(t“t')d(i-0  + 

JD/  D 

?/  drtp+{r)  {ih)di^*_(r)  J  dr*  ^+(r/)V;- 


„2 


x  (-ih)djJ2  p(R)  [ (p'\Ai(r,  t)\R)*(R'\Aj(rf,  t')\R)  d(t  —  t')  + 

^  Y.  J  dr  ip+{r)^*_{r)  J  d?  ^+(^)  iMO 

ij 

xdidj^PiR)  j  (R'\Ai(r,  t)\R)*  (R,\Aj(r' ,  t')\R)  d(t  —  t')  +  i  -H-  j 

R'R  ' 

This  expression  can  be  simplified  further  by  using  the  completeness  relation  for  the  reservoir 
final  states  and  the  following  definition  of  the  equilibrium  cross  correlation  function 


({Ai(f,t)  ,Aj(r,t')})  =  E  P(R)  (i?| { Ai (r,  t) ,  Aj {f ,  t')}\R) ,  and  its  Fourier  transform 

R 

J  e*WT-  ({Hj  (r,  r)  ,  Aj  (r^,  0) })  dr  =  Sij  (f,f',Lo')  ,  where  {  ,  }  is  anticommutator. 

Here  Sij  is  the  spectral  density  of  the  vector  potential.  This  simplifies  the  above  expression 
for  the  relaxation  rate  to 

f  dr  i>+{r)  f  dr'  (—ih)dji/j-(r')  Sij  (f,r* ,ui)  — 


m 


'2H2 


*7 


m 


e2h2 
4m*2  h2 


j_j  yy  J  dr ip+(f)  {ih)dj^*_{r)  J  d?  V’+(^)V’-(^)  {~iK)djSij  (f,? ,w)  - 

ij 

yy  J  drip+(r)ij)*L(r)  J  df1  ^(r')'ip-(r)  didjSij  ,oj)  .  (C.3) 

ij 


Ill 


The  above  expression  can  also  be  written  compactly  as 


r-^+  h2 


^  y,  J  dr  J  dr'  M*l{r)  M-}.(r')Sij(r,r'  ,u) ,  where  Mr(r)  is  defined  as 
ij 


Mr(r)  =  -e-t:'ip*+(r)P'ip-{r)  ~  ^^(r)  ^-(r)  V  . 
m  2m 


C.2  Derivation  of  Green’s  tensor  for  a  thin  film 


Here  we  present  the  calculation  for  the  retarded  photon  Green’s  tensor  outside  of  a  thin 


conducting  sheet  of  permittivity  e.  Green’s  function  will  satisfy  Eq.  B.l  Here  f'  is  simply 


a  parameter  for  the  purposes  of  solving  this  set  of  equations,  and  we  take  it  to  lie  in  the 
vacuum  outside  the  conducting  sheet.  We  will  suppress  the  dependence  of  Dik(r,r',oj)  on 
r'  and  u  to  simplify  the  notation.  The  geometry  of  the  problem  is  contained  entirely  in  the 
permittivity  function  e  (r,u).  We  take  the  boundaries  of  the  conducting  sheet  to  be  located 
at  z  =  —a  and  z  =  0,  with  vacuum  outside.  Because  the  geometry  is  still  translationally 


invariant  in  the  x  and  y  directions,  we  employ  the  same  Fourier  expansion  of  Eq.  7.13 


as  m 


Section  [7.2|  Solving  Eq.  B.l  for  a  problem  with  planar  symmetry  is  greatly  simplified  by 
separately  considering  the  Fourier  components  of  Eq.  |7.13  that  are  polarized  in  the  x  and 
y  directions.  Dyy{r)  may  then  be  reconstructed  as 

Dyy(r)  =  J  J^y2elk  H  ^Dyy^Ak,  z)  cos2  9  +  Dyytky(k,z)  sin2  0^j  ,  (C.4) 

where  Dyy^x  =  Dyy  when  ky  =  0,  and  Dyy  ^y  =  Dyy  when  kx  =  0.  The  boundary  value 
problem  for  Dyy  kx  (k,  z)  then  becomes 


Dyy:kx  (k,  z)  —  < 


Ae~iqz 
Bie~iqiz  +  B2eiqiz 
Ceiqz  +  5 \-qeiq\z~z'\ 


z  <  —a 
—a  <  z  <  0 
z  >  0 


(C.5) 


J  ~  — 


Our  interest  lies  in  the  behavior  of  the  fields  for  z  >  0,  so  we  need  only  to  calculate  C. 
Enforcing  that  Dyy  ^x  and  dDyy^x/dz  are  continuous  across  the  boundaries  results  in 


C  =  —rs(k,u,ci)eiqz 
2  q 


(C.6) 
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where 


rs(k,  uj,  a)  = 


(g2  -  qf)  sin(gia) 


(9?  +  S2)  sin(gia)  +  2iqq\  cos(gia) 


^iq\gQ  9l  _  ^—iqia  9  ~f~  1/1 

g  +  gi  g  -  gi 


=  2isingia  e 


-l 


(C.7) 


so  that 


Dyy>kx(k,z )  =  ^  (rs(fc,u;,a)ei(^+z')  +  ewk-*'I 


(C.8) 


The  term  proportional  to  exp(ig|z  —  z'l)  is  the  free  photon  contribution  to  the  power  spec¬ 
trum.  It  will  have  an  imaginary  component  and  thus  contribute  to  relaxation  only  in  the 
radiative  regime,  k  <  u>/c.  Within  a  skin  depth  of  separation  from  the  metal,  evanescent 
waves  are  orders  of  magnitude  larger  in  field  strength  than  these  free  photons.  They  may 
be  safely  ignored  in  this  context.  A  similar  calculation  yields  the  result  for  Dyy)ky: 


9 

ic  q 


Dyy)ky(k,z)  =  -7^2  [rp(k,u,a)e 


Jq{z+z')  _  eiq\z-z'\ 


(C.9) 


where 


rp(k,uj,a)  = 


(e2g2  -  g2)  sin(gia) 


fe2  +  eV)  sin(gia)  +  2iqq\e  cos(gia 
=  2zsingia  (  e 


«gi«M  _  c-iq,a^Q  +  Q  1 


-1 


(C.10) 


eg  +  gi  eg  -  gi 

A  Taylor  expansion  of  Eqs.  C.7  and|C.10|for  large  a  in  the  evanescent  range  of  wavevectors, 
i.e. ,  a  Taylor  expansion  in  powers  of  exp(— 2|gi|a),  gives  a  monotonically  increasing  function 
of  film  thickness,  a.  However,  a  more  careful  treatment  reveals  that  this  is  an  error.  The 


naive  expansions  of  C.7  and  C.10  for  large  a  neglect  an  enhancement  of  the  field  spectrum 
that  occurs  for  small  k.  In  fact,  EWJN  is  enhanced  as  the  thickness  is  decreased  for  any 
good  conductor.  Specifically,  the  enhancement  is  preserved  for  a  particular  spatial  Fourier 
component  of  the  Green’s  function  as  long  as  |  <  1.  EWJN  will  eventually  vanish  as 

a  — >  0,  but  this  does  not  occur  until  an  unphysically  small  thickness  is  reached,  on  the 
order  of  10~14  m  for  copper  at  T  =  0A'  which  is  well  below  the  applicability  of  the  local 
permittivity  model. 


113 


Bibliography 


[1]  R.  P.  Feynman,  “Simulating  physics  with  computers,”  Int.  J.  Theor.  Phys.,  vol.  21, 
no.  6/7,  p.  467,  1982. 

[2]  D.  Deutsch,  “Quantum  theory,  the  church-turing  principle  and  the  universal  quantum 
computer,”  Proc.  R.  Soc.  Lond.  A,  vol.  400,  pp.  97-117,  1985. 

[3]  P.  W.  Shor,  “Polynomial-time  algorithms  for  prime  factorization  and  discrete  loga¬ 
rithms  on  a  quantum  computer,”  SIAM  J.  Compute  vol.  26,  p.  1484,  1997. 

[4]  L.  K.  Grover,  “Quantum  mechanics  helps  in  searching  for  a  needle  in  a  haystack,” 
Phys.  Rev.  Lett.,  vol.  79,  p.  325,  1997. 

[5]  M.  A.  Nielsen  and  I.  L.  Chuang,  Quantum  Computation  and  Quantum  Information. 
Cambridge  University  Press,  2000. 

[6]  J.  I.  Cirac  and  P.  Zoller,  “Quantum  computations  with  cold  trapped  ions,”  Phys.  Rev. 
Lett.,  vol.  74,  p.  4091,  1995. 

[7]  C.  Monroe,  D.  M.  Meekhof,  B.  E.  King,  W.  M.  Itano,  and  D.  J.  Wineland,  “Demon¬ 
stration  of  a  fundamental  quantum  logic  gate,”  Phys.  Rev.  Lett.,  vol.  75,  p.  4714, 
1995. 

[8]  H.-J.  Briegel,  T.  Calarco,  D.  Jaksch,  J.  I.  Cirac,  and  P.  Zoller,  “Quantum  computing 
with  neutral  atoms,”  Journal  of  Modern  Optics,  vol.  47,  p.  415,  2000. 

[9]  G.  K.  Brennen,  C.  M.  Caves,  P.  S.  Jessen,  and  I.  H.  Deutsch,  “Quantum  logic  gates 
in  optical  lattices,”  Phys.  Rev.  Lett.,  vol.  82,  p.  1060,  1999. 

[10]  E.  Knill,  R.  Laflamme,  and  G.  J.  Milburn,  “A  scheme  for  efficient  quantum  compu¬ 
tation  with  linear  optics,”  Nature,  vol.  409,  p.  46,  2001. 

[11]  B.  E.  Kane,  “A  silicon-based  nuclear  spin  quantum  computer,”  Nature,  vol.  393, 
p.  133,  2001. 


114 


[12]  A.  P.  Nizovtsev,  S.  Y.  Kilin,  F.  Jelezko,  T.  Gaebal,  I.  Popa,  A.  Gruber,  and 

J.  Wrachtrup,  “A  quantum  computer  based  on  nv  centers  in  diamond:  Optically  de¬ 
tected  nutations  of  single  electron  and  nuclear  spins,”  Opt.  Spectrosc .,  vol.  99,  p.  233, 
2005. 

[13]  D.  Loss  and  P.  D.  DiVincenzo,  “Quantum  computation  with  quantum  dots,”  Phys. 
Rev.  A,  vol.  57,  p.  120,  1998. 

[14]  J.  Clarke  and  F.  K.  Wilhelm,  “Superconducting  quantum  bits,”  Nature,  vol.  453, 
p.  1031,  2008. 

[15]  L.  Chirolli  and  G.  Burkard,  “Decoherence  in  solid-state  qubits,”  Advances  in  Physics , 
vol.  57,  p.  225,  2008. 

[16]  J.  M.  Martinis,  K.  B.  Cooper,  R.  McDermott,  M.  Steffen,  M.  Ansmann,  K.  D.  Osborn, 

K.  Cicak,  S.  Oh,  D.  P.  Pappas,  R.  W.  Simmonds,  and  C.  C.  Yu,  “Decoherence  in 
josephson  qubits  from  dielectric  loss,”  Phys.  Rev.  Lett.,  vol.  95,  p.  210503,  2005. 

[17]  D.  A.  Lidar,  I.  L.  Chuang,  and  K.  B.  Whaley,  “Decoherence-free  subspaces  for  quan¬ 
tum  computation,”  Phys.  Rev.  Lett.,  vol.  81,  p.  2594,  1998. 

[18]  S.  D.  Sarma,  M.  Freedman,  and  C.  Nayak,  “Topologically-protected  qubits  from  a  pos¬ 
sible  non-abelian  fractional  quantum  hall  state,”  Phys.  Rev.  Lett.,  vol.  94,  no.  166802, 
2005. 

[19]  A.  Stern  and  N.  H.  Lindner,  “Topological  quantum  computation — from  basic  concepts 
to  first  experiments,”  Science,  vol.  339,  p.  1179,  2013. 

[20]  E.  Knill  and  R.  Laflamme,  “Theory  of  quantum  error-correcting  codes,”  Phys.  Rev. 
A,  vol.  55,  p.  900,  1997. 

[21]  L.  Viola,  E.  Knill,  and  S.  Lloyd,  “Dynamical  decoupling  of  open  quantum  systems,” 
Phys.  Rev.  Lett.,  vol.  82,  no.  2417,  1999. 

[22]  D.  Gottesman,  “Theory  of  fault-tolerant  quantum  computation,”  Phys.  Rev.  A, 
vol.  57,  p.  127,  1998. 

[23]  H.  K.  Lo,  S.  Popescu,  and  T.  P.  Spiller,  eds.,  Introduction  to  Quantum  Computation 
and  Information.  World  Scientific  Publishing  Co.  Pte.  Ltd,  1998. 

[24]  Y.  Liu,  L.  F.  Wei,  and  F.  Nori,  “Generation  of  nonclassical  photon  states  using  a 
superconducting  qubit  in  a  microcavity,”  Europhys.  Lett.,  vol.  67,  p.  941,  2004. 


115 


[25]  M.  Hofheinz,  E.  Weig,  M.  Ansmann,  R.  C.  Bialczak,  E.  Lucero,  M.  Neeley,  A.  D. 
O’Connell,  H.  Wang,  J.  M.  Martinis,  and  A.  N.  Cleland,  “Generation  of  fock  states 
in  a  superconducting  quantum  circuit,”  Nature ,  vol.  454,  p.  310,  2008. 

[26]  J.  Majer,  J.  M.  Chow,  J.  M.  Gambetta,  J.  Koch,  B.  R.  Johnson,  J.  A.  Schreier, 
L.  Frunzio,  D.  I.  Schuster,  A.  A.  Houck,  A.  Wallraff,  A.  Blais,  M.  H.  Devoret,  S.  M. 
Girvin,  and  R.  J.  Schoelkopf,  “Coupling  superconducting  qubits  via  a  cavity  bus,” 
Nature ,  vol.  449,  p.  443,  2007. 

[27]  M.  D.  Eisaman,  J.  Fan,  A.  Migdall,  and  S.  V.  Polyakov,  “Single-photon  sources  and 
detectors,”  Rev.  Sci.  Instrum.,  vol.  82,  p.  071101,  2011. 

[28]  A.  Ekert,  M.  Ericsson,  P.  Hayden,  H.  Inamori,  J.  A.  Jones,  D.  K.  L.  Oi,  and  V.  Vedral, 
“Geometric  quantum  computation,”  J.  Mod.  Opt.,  vol.  47,  p.  2501,  2000. 

[29]  P.  J.  Leek,  J.  M.  Fink,  A.  Blais,  R.  Bianchetti,  M.  Goppl,  J.  M.  Gambetta,  D.  I. 
Schuster,  L.  Frunzio,  R.  J.  Schoelkopf,  and  A.  Wallraff,  “Observation  of  berry’s  phase 
in  a  solid-state  qubit,”  Science,  vol.  318,  p.  1889,  2007. 

[30]  V.  Gritsev  and  A.  Polkovnikov,  “Dynamical  quantum  hall  effect  in  the  parameter 
space,”  PNAS,  vol.  109,  p.  6457,  2012. 

[31]  A.  Poudel  and  M.  G.  Vavilov,  “Effect  of  an  ohmic  environment  on  an  optimally 
controlled  flux-biased  phase  qubit,”  Phys.  Rev.  B,  vol.  82,  p.  144528,  2010. 

[32]  A.  Poudel,  R.  McDermott,  and  M.  G.  Vavilov,  “Quantum  efficiency  of  a  microwave 
photon  detector  based  on  a  current-biased  josephson  junction,”  Phys.  Rev.  B,  vol.  86, 
p.  174506,  2012. 

[33]  L.  S.  Langsjoen,  A.  Poudel,  M.  G.  Vavilov,  and  R.  Joynt,  “Qubit  relaxation  from 
evanescent-wave  johnson  noise,”  Phys.  Rev.  A,  vol.  86,  p.  010301(R),  2012. 

[34]  A.  Poudel,  L.  S.  Langsjoen,  M.  G.  Vavilov,  and  R.  Joynt,  “Relaxation  in  quantum 
dot  due  to  evanescent-wave  johnson  noise  from  a  metallic  backgate,”  Phys.  Rev.  B, 
vol.  87,  p.  045301,  2013. 

[35]  M.  H.  Devoret  and  J.  M.  Martinis,  “Implementing  qubits  with  superconducting  inte¬ 
grated  circuits,”  Quantum  Inf.  Process.,  vol.  3,  p.  163,  2004. 

[36]  M.  Constantin  and  C.  C.  Yu,  “Microscopic  model  of  critical  current  noise  in  josephson 
junctions,”  Phys.  Rev.  Lett.,  vol.  99,  p.  207001,  2007. 


116 


[37]  R.  Sousa,  K.  B.  Whaley,  T.  Hecht,  J.  Delft,  and  F.  K.  Wilhelm,  “Microscopic  model 
of  critical  current  noise  in  josephson-j unction  qubits:  Subgap  resonances  and  andreev 
bound  states,”  Phys.  Rev.  B,  vol.  80,  p.  094515,  2009. 

[38]  E.  Paladino,  L.  Faoro,  G.  Falci,  and  R.  Fazio,  “Decoherence  and  1/f  noise  in  josephson 
qubits,”  Phys.  Rev.  Lett.,  vol.  88,  p.  228304,  2002. 

[39]  L.  Faoro,  J.  Bergli,  B.  L.  Altshuler,  and  Y.  M.  Galperin,  “Models  of  environment  and 
tl  relaxation  in  josephson  charge  qubits,”  Phys.  Rev.  Lett.,  vol.  95,  p.  046805,  2005. 

[40]  R.  W.  Simmonds,  K.  M.  Lang,  D.  A.  Hite,  S.  Nam,  D.  P.  Pappas,  and  J.  M.  Martinis, 
“Decoherence  in  josephson  phase  qubits  from  junction  resonators,”  Phys.  Rev.  Lett., 
vol.  93,  p.  077003,  2004. 

[41]  S.  Sendelbach,  D.  Hover,  A.  Kittel,  M.  Muck,  J.  M.  Martinis,  and  R.  McDer¬ 
mott,  “Magnetism  in  squids  at  millikelvin  temperatures,”  Phys.  Rev.  Lett.,  vol.  100, 
p.  227006,  2008. 

[42]  L.  Faoro  and  L.  B.  Ioffe,  “Microscopic  origin  of  low-frequency  flux  noise  in  josephson 
circuits,”  Phys.  Rev.  Lett.,  vol.  100,  p.  227005,  2008. 

[43]  T.  Fujisawa,  T.  Hayashi,  and  S.  Sasaki,  “Time-dependent  single-electron  transport 
through  quantum  dots,”  Rep.  Prog.  Phys.,  vol.  69,  p.  759,  2006. 

[44]  T.  Fujisawa,  T.  Hayashi,  H.  D.  Cheong,  Y.  H.  Jeong,  and  Y.  Hirayama,  “Rotation 
and  phase-shift  operations  for  a  charge  qubit  in  a  double  quantum  dot,”  Physica  E, 
vol.  21,  p.  1046,  2004. 

[45]  X.  Hu,  B.  Koiller,  and  S.  D.  Sarma,  “Charge  qubits  in  semiconductor  quantum 
computer  architecture:  Tunnel  coupling  and  decoherence,”  Phys.  Rev.  B,  vol.  71, 
p.  235332,  2005. 

[46]  S.  W.  Jung,  T.  Fujisawa,  Y.  Hirayama,  and  Y.  H.  Jeong,  “Background  charge  fluctu¬ 
ation  in  a  gaas  quantum  dot  device,”  Appl.  Phys.  Lett.,  vol.  85,  p.  768,  2004. 

[47]  Y.  A.  Pashkin,  O.  Astafiev,  T.  Yamamoto,  Y.  Nakamura,  and  J.  S.  Tsai,  “Josephson 
charge  qubits:  a  brief  review,”  Quantum  Inf.  Process.,  vol.  8,  p.  55,  2009. 

[48]  Z.  Shi,  C.  B.  Simmons,  D.  R.  Ward,  J.  R.  Prance,  R.  T.  Mohr,  T.  S.  Koh,  J.  K.  Gam¬ 
ble,  X.  Wu,  D.  E.  Savage,  M.  G.  Lagally,  M.  Friesen,  S.  N.  Coppersmith,  and  M.  A. 
Eriksson,  “Coherent  quantum  oscillations  in  a  silicon  charge  qubit,”  arXiv:1208.0519, 
2012. 


117 


[49]  I.  V.  Yurkevich,  J.  Baldwin,  I.  V.  Lerner,  and  B.  L.  Altshuler,  “Decoherence  of  charge 
qubit  coupled  to  interacting  background  charges,”  Phys.  Rev.  B ,  vol.  81,  p.  121305, 
2010. 

[50]  S.  Vorojtsov,  E.  R.  Mucciolo,  and  H.  U.  Baranger,  “Phonon  decoherence  of  a  double 
quantum  dot  charge  qubit,”  Phys.  Rev.  B,  vol.  71,  p.  205322,  2005. 

[51]  D.  C.  B.  Valente,  E.  R.  Mucciolo,  and  F.  K.  Wilhelm,  “Decoherence  by  electro¬ 
magnetic  fluctuations  in  double-quantum-dot  charge  qubits,”  Phys.  Rev.  B ,  vol.  82, 
p.  125302,  2010. 

[52]  D.  P.  DiVincenzo,  D.  Bacon,  J.  Kempe,  G.  Burkard,  and  K.  B.  Whaley,  “Universal 
quantum  computation  with  the  exchange  interaction,”  Nature ,  vol.  408,  p.  339,  2000. 

[53]  M.  Friesen,  C.  Tahan,  R.  Joynt,  and  M.  A.  Eriksson,  “Spin  readout  and  initialization 
in  a  semiconductor  quantum  dot,”  Phys.  Rev.  Lett.,  vol.  92,  p.  037901,  2004. 

[54]  C.  Tahan,  M.  Friesen,  and  R.  Joynt,  “Decoherence  of  electron  spin  qubits  in  si-based 
quantum  computers,”  Phys.  Rev.  B,  vol.  66,  p.  035314,  2002. 

[55]  R.  Hanson,  L.  P.  Kouwenhoven,  J.  R.  Petta,  S.  Tarucha,  and  L.  M.  K.  Vandersypen, 
“Spins  in  few-electron  quantum  dots,”  Rev.  Mod.  Phys.,  vol.  79,  p.  1217,  2007. 

[56]  A.  V.  Khaetskii,  L.  Loss,  and  L.  Glazman,  “Electron  spin  decoherence  in  quantum 
dots  due  to  interaction  with  nuclei,”  Phys.  Rev.  Lett.,  vol.  88,  p.  186802,  2002. 

[57]  A.  V.  Khaetskii  and  Y.  V.  Nazarov,  “Spin  relaxation  in  semiconductor  quantum  dots,” 
Phys.  Rev.  B,  vol.  61,  p.  12639,  2000. 

[58]  S.  Erlingsson,  Y.  V.  Nazarov,  and  V.  I.  Fal’ko,  “Nucleus-mediated  spin-flip  transitions 
in  gaas  quantum  dots,”  Phys.  Rev.  B,  vol.  64,  p.  195306,  2001. 

[59]  S.  I.  Erlingsson  and  Y.  V.  Nazarov,  “Hyperfine-mediated  transitions  between  a  zeeman 
split  doublet  in  gaas  quantum  dots:  The  role  of  the  internal  field,”  Phys.  Rev.  B, 
vol.  66,  p.  155327,  2002. 

[60]  H.-P.  Breuer  and  F.  Petruccione,  The  theory  of  open  quantum  systems.  Oxford  Uni¬ 
versity  Press,  2002. 

[61]  A.  O.  Caldeira  and  A.  J.  Leggett,  “Quantum  tunnelling  in  a  dissipative  system,”  Ann. 
Phys.,  vol.  149,  p.  374,  1983. 


118 


[62]  J.  P.  Paz  and  W.  H.  Zurek,  “Course  8:  Environment-induced  decoherence  and  the 
transition  from  quantum  to  classical,”  in  Coherent  Atomic  Matter  Waves  (R.  Kaiser, 
C.  Westbrook,  and  F.  David,  eds.),  2001. 

[63]  U.  Weiss,  Quantum  Dissipative  Systems.  World  Scientific  Publishing  Co.  Pte.  Ltd, 
second  ed.,  1999. 

[64]  D.  Vion,  A.  Aassime,  A.  Cottet,  P.  Joyez,  H.  Pothier,  C.  Urbina,  D.  Esteve,  and 
M.  H.  Devoret,  “Manipulating  the  quantum  state  of  an  electrical  circuit,”  Science , 
vol.  296,  p.  886,  2002. 

[65]  J.  M.  Martinis,  S.  Nam,  J.  Aumentado,  and  C.  Urbina,  “Rabi  oscillations  in  a  large 
josephson-j unction  qubit,”  Phys.  Rev.  Lett.,  vol.  89,  p.  117901,  2002. 

[66]  I.  Chiorescu,  Y.  Nakamura,  C.  J.  P.  M.  Harmans,  and  J.  E.  Mooij,  “Coherent  quantum 
dynamics  of  a  superconducting  flux  qubit,”  Science,  vol.  299,  p.  1869,  2003. 

[67]  F.  Motzoi,  J.  M.  Gambetta,  P.  Rebentrost,  and  F.  K.  Wilhelm,  “Simple  pulses 
for  elimination  of  leakage  in  weakly  nonlinear  qubits,”  Phys.  Rev.  Lett.,  vol.  103, 
p.  110501,  2009. 

[68]  E.  Knill,  “Quantum  computing  with  realistically  noisy  device,”  Nature,  vol.  434,  p.  44, 
2005. 

[69]  M.  Steffen,  J.  M.  Martinis,  and  I.  L.  Chuang,  “Accurate  control  of  josephson  phase 
qubits,”  Phys.  Rev.  B ,  vol.  68,  p.  224518,  2003. 

[70]  P.  Rebentrost,  I.  Serban,  T.  Schulte-Herbriiggen,  and  F.  K.  Wilhelm,  “Optimal  control 
of  a  qubit  coupled  to  a  non-markovian  environment,”  Phys.  Rev.  Lett.,  vol.  102, 
p.  090401,  2009. 

[71]  P.  Rebentrost  and  F.  K.  Wilhelm,  “Optimal  control  of  a  leaking  qubit,”  Phys.  Rev. 
B,  vol.  79,  p.  060507,  2009. 

[72]  G.  Lindblad,  “On  the  generators  of  quantum  dynamical  semigroups,”  Comm.  Math. 
Phys.,  vol.  48,  p.  119,  1976. 

[73]  D.  J.  Van  Harlingen,  T.  L.  Robertson,  B.  L.  T.  Plourde,  P.  A.  Reichardt,  T.  A.  Crane, 
and  J.  Clarke,  “Decoherence  in  josephson-junction  qubits  due  to  critical-current  fluc¬ 
tuations,”  Phys.  Rev.  B,  vol.  70,  p.  064517,  2004. 


119 


[74]  J.  Claudon,  A.  Fay,  L.  P.  Levy,  and  O.  Buisson,  “Decoherence  processes  in  a  current 
biased  dc  squid,”  Phys.  Rev.  B,  vol.  73,  p.  180502,  2006. 

[75]  R.  C.  Bialczak,  R.  McDermott,  M.  Ansmann,  M.  Hofheinz,  N.  Katz,  E.  Lucero, 
M.  Neeley,  A.  D.  O’Connell,  H.  Wang,  A.  N.  Cleland,  and  J.  M.  Martinis,  “1//  flux 
noise  in  josephson  phase  qubits,”  Phys.  Rev.  Lett.,  vol.  99,  p.  187006,  2007. 

[76]  J.  M.  Martinis,  “Superconducting  phase  qubits,”  Quantum  Inf.  Process.,  vol.  8,  p.  81, 
2009. 

[77]  J.  Lisenfeld,  A.  Lukashenko,  M.  Ansmann,  J.  M.  Martinis,  and  A.  V.  Ustinov,  “Tem¬ 
perature  dependence  of  coherent  oscillations  in  josephson  phase  qubits,”  Phys.  Rev. 
Lett.,  vol.  99,  p.  170504,  2007. 

[78]  M.  D.  Bowdrey,  D.  K.  L.  Oi,  A.  J.  Short,  K.  Banaszek,  and  J.  A.  Jones,  “Fidelity  of 
single  qubit  maps,”  Phys.  Lett.  A,  vol.  294,  p.  258,  2002. 

[79]  J.  M.  Chow,  L.  DiCarlo,  J.  M.  Gambetta,  F.  Motzoi,  L.  Frunzio,  S.  M.  Girvin,  and 

R.  J.  Schoelkopf,  “Implementing  optimal  control  pulse  shaping  for  improved  single¬ 
qubit  gates,”  Phys.  Rev.  A,  vol.  82,  p.  040305,  2010. 

[80]  D.  Walls  and  G.  Milburn,  Quantum  Optics.  Springer- Verlang,  Berlin,  1994. 

[81]  P.  Kok,  W.  J.  Munro,  K.  Nernoto,  T.  C.  Ralph,  J.  P.  Dowling,  and  G.  J.  Milburn, 
“Linear  optical  quantum  computing  with  photonic  qubits,”  Rev.  Mod.  Phys.,  vol.  79, 
p.  135,  2007. 

[82]  R.  H.  Hadfield,  “Single-photon  detectors  for  optical  quantum  information  applica¬ 
tions,”  Nature  Photonics,  vol.  3,  p.  696,  2009. 

[83]  A.  Blais,  R.-S.  Huang,  A.  Wallraff,  S.  M.  Girvin,  and  R.  J.  Schoelkopf,  “Cavity 
quantum  electrodynamics  for  superconducting  electrical  circuits:  An  architecture  for 
quantum  computation”,”  Phys.  Rev.  A,  vol.  69,  p.  062320,  2004. 

[84]  A.  Wallraff,  D.  I.  Schuster,  A.  Blais,  L.  Frunzio,  R.  S.  Huang,  J.  Majer,  S.  Kumar, 

S.  M.  Girvin,  and  R.  J.  Schoelkopf,  “Strong  coupling  of  a  single  photon  to  a  super¬ 
conducting  qubit  using  circuit  quantum  electrodynamics,”  Nature,  vol.  431,  p.  162, 
2004. 

[85]  J.  Q.  You  and  F.  Nori,  “Atomic  physics  and  quantum  optics  using  superconducting 
circuits,”  Nature,  vol.  474,  p.  589,  2011. 


120 


[86]  J.  Q.  You  and  F.  Nori,  “Quantum  information  processing  with  superconducting  qubits 
in  a  microwave  field,”  Phys.  Rev.  B,  vol.  68,  p.  064509,  2003. 

[87]  I.  Buluta,  S.  Ashhab,  and  F.  Nori,  “Natural  and  artifical  atoms  for  quantum  compu¬ 
tation,”  Rep.  Prog.  Phys.,  vol.  74,  p.  104401,  2011. 

[88]  A.  A.  Houck,  D.  I.  Schuster,  J.  M.  Gambetta,  J.  A.  Schreier,  B.  R.  Johnson,  J.  M. 
Chow,  L.  Frunzio,  J.  Majer,  M.  H.  Devoret,  S.  M.  Girvin,  and  R.  J.  Schoelkopf, 
“Generating  single  microwave  photons  in  a  circuit,”  Nature,  vol.  449,  p.  328,  2007. 

[89]  I.  Chiorescu,  P.  Bertet,  K.  Semba,  Y.  Nakamura,  C.  J.  P.  M.  Harmans,  and  J.  E. 
Mooij,  “Coherent  dynamics  of  a  flux  qubit  coupled  to  a  harmonic  oscillator,”  Nature, 
vol.  431,  p.  159,  2004. 

[90]  G.  Romero,  J.  J.  Garcia-Ripoll,  and  E.  Solano,  “Microwave  photon  detector  in  circuit 
qed,”  Phys.  Rev.  Lett.,  vol.  102,  p.  173602,  2009. 

[91]  Y.-F.  Chen,  D.  Hover,  S.  Sendelbach,  L.  Maurer,  S.  T.  Merkel,  E.  J.  Pritchett,  F.  K. 
Wilhelm,  and  R.  McDermott,  “Microwave  photon  counter  based  on  josephson  junc¬ 
tions,”  Phys.  Rev.  Lett.,  vol.  107,  p.  217401,  2011. 

[92]  N.  Gisin,  G.  Ribordy,  W.  Tittel,  and  H.  Zbinden,  “Quantum  cryptography,”  Rev. 
Mod.  Phys.,  vol.  74,  p.  145,  2002. 

[93]  B.  Osberg,  “Quantum  microwave  photodetection  using  superconducting  josephson 
circuits,”  University  of  Waterloo  (Msc.  Thesis),  http://hdl.handle.net/10012/4689,, 
2009. 

[94]  B.  Peropadre,  G.  Romero,  G.  Johansson,  C.  Wilson,  E.  Solano,  and  J.  J.  Garcia- 
Ripoll,  “Approaching  perfect  microwave  photodetection  in  circuit  qed,”  Phys.  Rev. 

A,  vol.  84,  no.  6,  p.  063834,  2011. 

[95]  J.  Q.  You,  Y.-X.  Liu,  C.  P.  Sun,  and  F.  Nori,  “Persistent  single-photon  production 
by  tunable  on-chip  micromaser  with  a  superconducting  quantum  circuit,”  Phys.  Rev. 

B,  vol.  75,  p.  104516,  2007. 

[96]  M.  H.  Devoret,  J.  M.  Martinis,  D.  Esteve,  and  J.  Clarke,  “Resonant  activation 
from  the  zero- voltage  state  of  a  current-biased  josephson  junction,”  Phys.  Rev.  Lett., 
vol.  53,  p.  1260,  1984. 


[97]  E.  Jaynes  and  F.  Cummings,  “Comparison  of  quantum  and  semiclassical  radiation 
theories  with  application  to  the  beam  maser,”  Proceedings  of  the  IEEE,  vol.  51,  p.  89, 
1963. 


121 


[98]  R.  R.  Puri  and  G.  S.  Agarwal,  “Collapse  and  revival  phenomena  in  the  jaynes- 
cummings  model  with  cavity  damping,”  Phys.  Rev.  A,  vol.  33,  p.  3610,  1986. 

[99]  S.  Gurvitz  and  Y.  S.  Prager,  “Microscopic  derivation  of  rate  equations  for  quantum 
transport,”  Phys.  Rev.  B ,  vol.  53,  p.  15932,  1996. 

[100]  S.  Gurvitz,  “Rate  equations  for  quantum  transport  in  multidot  systems,”  Phys.  Rev. 
B,  vol.  57,  p.  6602,  1998. 

[101]  J.  Ping,  X.-Q.  Li,  and  S.  Gurvitz,  “Quantum  coherence  and  entanglement  induced 
by  the  continuum  between  distant  localized  states,”  Phys.  Rev.  A,  vol.  83,  p.  042112, 
2011. 

[102]  J.  Ping,  Y.  Ye,  X.-Q.  Li,  Y.  Yan,  and  S.  Gurvitz,  “Undetectable  quantum  transfer 
through  a  continuum,”  arXiv:1106.2901v2,  06  2011. 

[103]  A.  A.  Clerk,  M.  H.  Devoret,  S.  M.  Girvin,  F.  Marquardt,  and  R.  J.  Schoelkopf, 
“Introduction  to  quantum  noise,  measurement,  and  amplification,”  Rev.  Mod.  Phys., 
vol.  82,  p.  1155,  2010. 

[104]  J.  R.  Johansson,  G.  Johansson,  C.  M.  Wilson,  and  F.  Nori,  “Dynamical  casimir  effect 
in  superconducting  microwave  circuits,”  Phys.  Rev.  A,  vol.  82,  p.  052509,  2010. 

[105]  T.  Brandes,  “Waiting  times  and  noise  in  single  particle  transport,”  Ann.  Phys.,  vol.  17, 
p.  477,  2008. 


[106]  M.  V.  Berry,  “Quantal  phase  factors  accompanying  adiabatic  changes,”  Proc.  R.  Soc. 
Lond.  A,  vol.  392,  no.  1802,  p.  45,  1984. 


[107]  L.-M.  Duan,  J.  I.  Cirac,  and  P.  Zoller,  “Geometric  manipulation  of  trapped  ions  for 
quantum  computation,”  Science,  vol.  292,  p.  1695,  2001. 


[108]  J.  P.  Provost  and  G.  Valle,  “Riemannian  structure  on  manifolds  of  quantum  states,” 
Comm.  Math.  Phys.,  vol.  76,  no.  3,  p.  289,  1980. 


[109]  S.  Berger,  M.  Pechal,  J.  A.  A.  Abdumalikov,  C.  Eichler,  L.  Steffen,  A.  Fedorov, 
A.  Wallraff,  and  S.  Filipp,  “Geometric  phases  in  superconducting  qubits  beyond  the 
two-level  approximation,”  Phys.  Rev.  B,  vol.  85,  p.  220502(R),  2012. 


[110]  S.  Berger,  M.  Pechal,  J.  A.  A.  Abdumalikov,  C.  Eichler,  L.  Steffen,  A.  Fedorov, 
A.  Wallraff,  and  S.  Filipp,  “Exploring  the  effect  of  noise  on  geometric  phases  using 
superconducting  qubits,”  arXiv:1302.3305,  2013. 


122 


[111]  J.  E.  Avron,  M.  Fraas,  G.  M.  Graf,  and  O.  Kenneth,  “Quantum  response  of  dephasing 
open  systems,”  New  J.  Phys.,  vol.  13,  p.  053042,  2011. 

[112]  R.  S.  Whitney,  Y.  Makhlin,  A.  Shnirman,  and  Y.  Gefen,  “Geometric  nature  of  the 
environment-induced  berry  phase  and  geometric  dephasing,”  Phys.  Rev.  Lett.,  vol.  94, 
p.  070407,  2005. 

[113]  T.  Hayashi,  T.  Fujisawa,  H.  D.  Cheong,  Y.  H.  Jeong,  and  Y.  Hirayama,  “Coherent 
manipulation  of  electronic  states  in  a  double  quantum  dot,”  Phys.  Rev.  Lett.,  vol.  91, 
p.  226804,  2003. 

[114]  J.  Petta,  A.  C.  Johnson,  C.  M.  Marcus,  M.  P.  Hanson,  and  A.  C.  Gossard,  “Manipula¬ 
tion  of  a  single  charge  in  a  double  quantum  dot,”  Phys.  Rev.  Lett.,  vol.  93,  p.  186802, 

2004. 

[115]  J.  Gorman,  D.  G.  Hasko,  and  D.  A.  Williams,  “Charge-qubit  operation  of  an  isolated 
double  quantum  dot,”  Phys.  Rev.  Lett.,  vol.  95,  p.  090502,  2005. 

[116]  Y.  Hu,  H.  O.  H.  Churchill,  D.  J.  Reilly,  J.  Xiang,  C.  M.  Lieber,  and  C.  M.  Marcus, 
“A  ge/si  heterostructure  nanowire-based  double  quantum  dot  with  integrated  charge 
sensor,”  Nature  Naotechnology,  vol.  2,  p.  622,  2007. 

[117]  N.  Mason,  M.  J.  Biercuk,  and  C.  M.  Marcus,  “Local  gate  control  of  a  carbon  nanotu- 
ube  double  quantum  dot,”  Science,  vol.  303,  p.  655,  2004. 

[118]  J.  Taylor,  J.  Petta,  A.  C.  Johnson,  A.  Yacoby,  C.  M.  Marcus,  and  M.  D.  Lukin, 
“Relaxation,  dephasing,  and  quantum  control  of  electron  spins  in  double  quantum 
dots,”  Phys.  Rev.  B,  vol.  76,  p.  035315,  2007. 

[119]  V.  N.  Stavrou  and  X.  Hu,  “Charge  decoherence  in  laterally  coupled  quantum  dots 
due  to  electron-phonon  interactions,”  Phys.  Rev.  B,  vol.  72,  p.  075362,  2005. 

[120]  M.  Thowart,  J.  Eckel,  and  E.  R.  Mucciolo,  “Non-markovian  dynamics  of  double  quan¬ 
tum  dot  charge  qubits  due  to  acoustic  phonons,”  Phys.  Rev.  B,  vol.  72,  p.  235320, 

2005. 

[121]  E.  M.  Lifshitz,  “The  theory  of  molecular  attractive  forces  between  solids,”  Sov.  Phys. 
JETP ,  vol.  2,  p.  73,  1956. 

[122]  G.  Agarwal,  “Quantum  electrodynamics  in  the  presence  of  dielectrics  and  conductors, 
iii.  relations  among  one-photon  transition  probabilities  in  stationary  and  nonstation¬ 
ary  fields,  density  of  states,  the  field-correlation  functions,  and  surface-dependent 
response  functions,”  Phys.  Rev.  A,  vol.  11,  p.  253,  1975. 


123 


[123]  S.  M.  Rytov,  Theory  of  Electrical  Fluctuation  and  Thermal  Radiation.  AN  SSSR 
Press,  1953. 

[124]  J.  B.  Pendry,  “Radiative  exchange  of  heat  between  nanostructures,”  J.  Rhys.:  Con- 
dens.  Matter ,  vol.  11,  p.  6621,  1999. 


[125]  A.  I.  Volokitin  and  B.  N.  J.  Persson,  “Near-held  radiative  heat  transfer  and  noncontact 
friction,”  Rev.  Mod.  Phys.,  vol.  79,  p.  1291,  2007. 


[126]  K.  Joulain,  J.-P.  Mulet,  F.  Marquier,  R.  Carminati,  and  J.-J.  Greffet,  “Surface  elec¬ 
tromagnetic  waves  thermally  excited:  Radiative  heat  transfer,  coherence  properties 
and  casimir  forces  revisited  in  the  near  held,”  Surf.  Sci.  Rep .,  vol.  57,  p.  59,  2005. 

[127]  C.  Henkel,  S.  Potting,  and  M.  Wilkens,  “Loss  and  heating  of  particles  in  small  and 
noisy  traps,”  Appl.  Phys.  B,  vol.  69,  p.  379,  1999. 


[128]  D.  M.  Harber,  J.  M.  McGuirk,  J.  M.  Obrecht,  and  E.  A.  Cornell,  “Thermally  induced 
losses  in  ultra-cold  atoms  magnetically  trapped  near  room-temperature  surfaces,”  J. 
Low  Temp.  Phys.,  vol.  133,  p.  229,  2003. 


[129]  F.  Marquardt  and  V.  A.  Abalmassov,  “Spin  relaxation  in  a  quantum  dot  due  to 
nyquist  noise,”  Phys.  Rev.  B,  vol.  71,  p.  165325,  2005. 

[130]  G.  Ford  and  W.  Weber,  “Electromagnetic  interactions  of  molecules  with  metal  sur¬ 
faces,”  Phys.  Rep.,  vol.  113,  p.  195,  1984. 


[131]  C.  Henkel  and  K.  Joulain,  “Electromagnetic  held  correlations  near  a  surface  with  a 
nonlocal  optical  response,”  Appl.  Phys.  B,  vol.  84,  p.  61,  2006. 


[132]  H.  B.  Callen  and  T.  A.  Welton,  “Irreversibility  and  generalized  noise,”  Phys.  Rev., 
vol.  83,  p.  34,  1951. 


[133]  E.  M.  Lifshitz  and  L.  P.  Pitaevskii,  Statistical  Physics,  Part  2,  vol.  9  of  Course  in 
Theoretical  Physics.  Pergamon,  1980. 


[134]  M.  Xiao,  M.  G.  House,  and  H.  W.  Jiang,  “Measurement  of  the  spin  relaxation  time 
of  single  electrons  in  a  silicon  metal-oxide-semiconductor-based  quantum  dot,”  Phys. 
Rev.  Lett.,  vol.  104,  p.  096801,  2010. 


[135]  J.  M.  Elzerman,  R.  Hanson,  L.  H.  Willems  van  Beveren,  B.  Witkamp,  L.  M.  K. 
Vandersypen,  and  L.  P.  Kouwenhoven,  “Single-shot  read-out  of  an  individual  electron 
spin  in  a  quantum  dot,”  Nature,  vol.  430,  p.  431,  2004. 


124 


[136]  M.  G.  Borselli,  K.  Eng,  E.  T.  Croke,  B.  M.  Maune,  B.  Huang,  R.  S.  Ross,  A.  A.  Kiselev, 
P.  W.  Deelman,  I.  Alvarado- Rodriguez,  A.  E.  Schmitz,  M.  Sokolich,  K.  S.  Holabird, 
T.  M.  Hazard,  M.  F.  Gyure,  and  A.  T.  Hunter,  “Pauli  spin  blockade  in  undoped  si/sige 
two-electron  double  quantum  dots,”  Appl.  Phys.  Lett .,  vol.  99,  p.  063109,  2011. 

[137]  K.  D.  Petersson,  J.  R.  Petta,  H.  Lu,  and  A.  C.  Gossard,  “Quantum  coherence  in  a 
one-electron  semiconductor  charge  qubit,”  Phys.  Rev.  Lett.,  vol.  105,  p.  246804,  2010. 

[138]  S.  Arnasha,  K.  MacLean,  I.  P.  Radu,  D.  M.  Zumbiihl,  M.  A.  Kastner,  M.  P.  Hanson, 
and  A.  C.  Gossard,  “Electrical  control  of  spin  relaxation  in  a  quantum  dot,”  Phys. 
Rev.  Lett.,  vol.  100,  p.  046803,  2008. 

[139]  A.  V.  Khaetskii  and  Y.  V.  Nazarov,  “Spin-flip  transitions  between  zeernan  sublevels 
in  semiconductor  quantum  dots,”  Phys.  Rev.  B.,  vol.  64,  p.  125316,  2001. 

[140]  B.  D.  Josephson,  “Possible  new  effects  in  superconductive  tunnelling,”  Phys.  Lett., 
vol.  1,  p.  251,  1962. 

[141]  M.  Tinkham,  Introduction  to  Superconductivity.  McGraw-Hill,  Inc.,  1980. 


